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ABSTRACT 


"^Rock  and  Joint  properties  were  reviewed  as  a  basis  for 
developing  the  analyses.  The  review  discussed:  failure  criteria 
for  the  rock;  and  Joint  properties.  Improved  methods  of  sampling 
and  testing  joints  were  explored  in  detail.  Pile  Driver  test 
sections  were  used  as  a  basis  for  analysis. 

Two  approaches  of  analysis  of  the  rock  to  predict  performance 
of  rock  bolted  and  unlined  test  sections  were  used.  The  first 
method  called  finite  element  was  used  on  one  section  only.  This 
method's  complexity  prevented  its  application  to  other  sections. 
Discontinuities,  heterogeneities,  and  local  failures  can  be  taken 
into  account  by  finite  elements,  but  the  representation  is  only 
two  dimensional.  Joint  influence  analysis  ,•  the  second  method, 
was  applied  to  many  sections.  In  comparison  to  the  finite  element 
method,  joint  influence  analysis  is  relatively  simple  to  use. 
Joints  and  rock  bolts  are  assumed  not  to  modify  the  stress  dis¬ 
tribution  in  the  joint  influence  analysis;  however,  the  method 
represents  a  three  dimensional  approach. 

Comparisons  were  made  to  post-shot  cross  sections  of  the  Pile 
Driver  openings  to  evaluate  the  methods  of  analysis ^  The  finite 


element  method  is  more  powerful;  the  joint  influence  analysis  is 


easier  to  use. 


PREFACE 


This  investigation  was  authorized  by  the  Chief  of  Engineers 
(ENGMC-EJ4)  and  was  performed  in  PST  1968  under  Contract  No.  MCA^5~ 
67-C-0015,  Mod.  P001,  between  the  Cmaha  District,  Corps  of  Engineers 
and  Dr.  Richard  E.  Goodman,  Berkeley,  California. 

ftiis  report  was  prepared  under  the  supervision  of  Dr.  R.  E. 
Goodman,  Principal  Investigator.  Personnel  contributing  to  the 
report  were  Michael  Cleary,  Ashraf  Mahtab,  and  Fred  Sage. 

During  the  work  period  covered  by  this  report.  Colonel 
William  H.  McKenzie  III  was  District  Engineer;-  Charles  L.  Hipp 
was  Chief,  Engineering  Division;  C.  J.  Distefano  was  technical 
monitor  for  the  Omaha  District  under  the  general  supervision 
of  D.  C.  Aardappel,  Acting  Chief,  Protective  Structures  Branch. 

D.  G.  Heltmann  participated  in  the  monitoring  work. 


TABLE  OP  CONTENTS 


ABSTRACT  i 

PREFACE  ii 

1.  INTRODUCTION  1 

2.  MATHEMATICAL  MODELING  OF  TUNNELS  IN  JOINTED  ROCK  5 

2.1  GENERAL-FINITE  ELEMENT  METHOD  5 

2.2  GEOLOGICAL  DATA  5 

2.3  ROCK  PROPERTIES  6 

(1)  Pcwer  Law  Failure  Criterion  7 

(2)  Failure  Law  When  o  t  if  Negative  12 

2.4  JOINT  PROPERTIES  32 

(1)  Importance  of  Joints  12 

(2)  Methods  of  Sampling  Joints  14 

(3)  Methods  of  Testing  Joints  18 

(4)  Results  of  Shear  Tests  22 

2.5  APPLICATION  TO  PIIEDKIVER  -  DL  DRIFT  24 

(1)  Case  1  —  Results  25 

(2)  Case  2  —  Results  26 

(3)  Comparison  with  Actual  Damage  26 

(4)  Conclusion  26 

3.  ANALYSIS  OF  CR  AND  DL  DRIFTS  BY  JOINT  INFLUENCE 

DIAGRAM  METHOD  58 

3.1  DESCRIPTION  OF  PROCEDURES  58 

3.2  CR  EAST  -  WEST  DRIFT  —  RESULTS  OF  ANALYSIS  60 

(1)  Station  0+60  6l 

(2)  Station  0  +  70  6l 

(3)  Station  0  +  80  62 

(4)  station  0  +  90  to  1  +  20  63 

(5)  Station  1  +  30,  1  +  40  64 

(6)  station  1  +  50  64 

(7)  Station  1  +  60  64 

(8)  station  1  4  70  65 


iii 


3*3  CR  NORTH  MOTT  —  RESULTS  OP  ANALYSIS  66 

(1)  Station  1  +  60  66 

(2)  Station  1+70  67 

(3)  Station  1  +  80  67 

(If)  Station  1+90  67 

(5)  Station  2  +  00  to  2  +  36  68 

3.4  DL  1  AND  2  —  RESULTS  OF  ANALYSIS  68 

3.5  SUMMARY  OP  RESULTS  70 

k.  CONCLUSIONS  AND  RECOMMENDATIONS  167 

PERSONNEL  AND  ACKNOWLEDGEMENTS  168 

APPENDIX  1  -  FINITE  ELEMENT  ANALYSIS  COMPUTER  PROGRAM  169 
APPENDIX  2  -  JOINT  INFLUENCE  ANALYSIS  COMPUTER  PROGRAM  183 


TABLES 

Page 

2.1 

Suraary  of  Rock  Properties 

28 

2.2 

Results  of  Laboratory  Direct  Shear  Tests  on  Rock 
Specimens  with  Joints 

29 

2.3 

Summary  of  Results  from  In  Situ  Block  Shear  Test 
on  Weakness  Planes 

32 

A 

2,4 

Typical  Load  Deformation  Relations  for  Different 

Types  of  seams  and  Fractures 

33 

:  - 

2.5 

Rock  and  Joint  Properties  Used  in  Study  of  Section 
at  0  +  70  in  DL  Drift 

34 

3.1 

Applied  Pressures  to  Develop  Stress  Distributions 

Used  in  Joint  Analysis 

71 

3.2 

CR(EW)-  Drift-Lining  Type  and  Tunnel  Dimensions 

72 

3.3 

Drift  CR  -  Lateral  -  Joint  Set  Strikes  and  Mp3  at 
Each  Station 

73 

3.4 

Brj.ft  CR  -  Lateral  -  Joint  Orientation  Parameters 

74 

3.5 

CR  North  Drift  -  Lining  Types  and  Tunnel  Dimensions 

76 

3.6 

CR  North  Drift  -  Joint  Set  Strikes  and  Dips  at  Each 
Station 

77 

3.7 

CR  North  Drift  -  Joint  Orientation  Itaameters 

78 

3-8 

DL1  and  DL2  Joint  Set  Orientation 

79 

3-9 

Predicted  and  Actual  Tunnel  Damage  —  CR(EW) 

80 

3*10  Predicted  and  Actual  Tunnel  Damage  —  CR(North) 

82 

3.11  Predicted  and  Actual  Tunnel  Damage  —  DL1  and  DI2 

83 

FIGURES 

1.1 

Perspective  of  Drift  Complex 

4 

2.1 

Pre-shot  Geologic  Log,  DLL  and  DL2 

35 

2.2 

Method  of  Projecting  Geological  Data 

37 

2.3 

Geologic  Cross  Sections  Drift  DL1  -  2 

38 

4 

r 

2.4 

Comparison  of  Projected  and  Actual  End  Wall  as  a 

Test  of  the  Method  of  Projection 

39 

9 

2.5 

Octahedral  Stresses  at  Failure  for  Piledriver  Rock 

40 

> 

2.6 

Log  -  Log  Pit  of  Octahedral  Stresses  at  Failure 

4l 

<* 

2.7 

Finite  Element  Mesh  FD5  Used  in  Example  Below 

42 

2.8 

Progressive  Failure  of  Elements  as  P^  Increases 
from  0  -  32000  psi.  in  8  Increments;  P2/P1  =  0.25 

42 

2.9 

Artificial  Joint  specimen 

43 

V 


FIGUBT2S  -  continued 


2.10  Preparation  of  Direct  Shear  Specimen  Containing 

an.  Artificial  Extension  Joint  1*3 

2.11  large  Direct  Shear  specimens  of  a  Natural  Joint 

Obtained  by  Drilling  7  l/2"  Cores  Across  Joint  44 

2.12  Joint  Specimens  Obtained  by  Drilling  ICC  Cores 

Parallel  to  the  Joint  Plane  44 

2.13  large  Core  specimens  Containing  Joints  in  the 

Diametral  Plane  44 

2.14  large  Direct  Shear  Specimens  Obtained  by  Drilling 

Nine  Inch  Cores  Parallel  to  Selected  Joint  Surfaces  45 

2.15  Wire  Sawing  Technique  Used  by  Coyne  and  Bellier  to 

Extract  Block  Samples  Containing  Joints  46 

2.16  large  Direct  Shear  Machine  at  Imperial  College,  Royal 

School  of  Mines  46 

2.17  Stress  Distribution  in  a  Direct  Shear  Box  47 

2.18  Casagrande  Type  Direct  Shear  Machine  Used  by  Coyne 

and  Belller  46 

2. IS  Plane  Strain  Direct  shear  Machine  Qnployed  at  Dept. 

of  Geology,  Imperial  College  48 

2.20  Plane  Strain  Direct  Shear  Machine  at  University  of 

Illinois  48 

2.21  Krsmanovic  Shear  Machine,  Sarajevo,  Yugoslavia  49 

2.22  In  Situ  Shear  Test  Conducted  at  Vouglans  Dam  50 

2.23  In  Situ  Shear  Tests  Conducted  in  Czechoslovakia  51 

2.24  In  Situ  Shear  Tests  on  a  Clay  Parting  along  Bedding 
in  Shale  on  a  Natural  Slope  at  Khajure  Kach  Dam, 

Pakistan  51 

2.25  Typical  Shear  Stress-Deformation  Relationships  for 

Various  Weakness  Surfaces  52 

2.26  Finite  Element  Mesh  for  DL  Drift  0  +  70  53 

2.27a  DL  Drift  O+70  Displacements  and  Failed  Elements 

for  Case  1  54 

2.27b  DL  Drift  0+70  Displacements  and  Failed  Elements 

for  Case  1  55 

2.28  DL  Drift  0+70  Displacements  and  Failed  Elements 

for  Case  2  56 

2.29  Post  Shock  Cross  Section  0  +  70  -  DL  Drift  57 

3.1a  General  Plan,  CR  Drift  84 

3.1b  Plan  of  Structural  Geology  in  Vicinity  of  CR  and 

DL  Drifts  85 

vi 


FIGURES  -  continued 


3.1c  Geological  Logs  CR  East  West  Drift  Station  0+37 

to  1  +  12  86 

3. Id  Geological  Logs  CR  East  West  Station  1  +  12  to  1  +  76  87 

3*2  Joint  Influence  Diagrams  CR  Lateral  -  10,000  psi. 

Blast,  c  -  100  psi  88 

3.3  Joint  Influence  Diagrams  CR  Lateral  -  10,000  psi. 

Blast,  c  =  1000  psi.  89 

3.1m  Geologic  Cross  Section,  0  +  60,  16 '  Diara.,  16’  Rock 

Bolts  at  2*  go 

Influential  Joints,  0  +  60,  16'  Diara.,  lc'  Rock 

Bolts  at  2*  go 

3.l4c,d  Predicted  Post-Shock  Cross  Sections,  0  +  60,  16' 

Diara.,  16’  Rock  Bolts  at  2'  91 

3.5a  Geologic  Cross  Section,  0  +  70,  16'  Diara.,  16’ 

Rock  Bolts  at  2'  g2 

3«5t  'Influential  Joints,  0  +  70,  16'  Diara.,  16'  Rock 

Bolts  at  2’  gp 

3*5c,d  Predicted  Post-Shock  Cross  Sections,  0  +  70,  lo’ 

Diam.,  l6f  Rock  Bolts  at  2'  93 

3.6a  Geologic  Cross  Section,  0  +  80,  16'  Diam.,  l6‘ 

Rock  Bolts  at  2'  9k 

3.6b  Influential  Joints,  0  +  80,  16'  Diam.,  16’  Rock 

Bolts  at  2'  9k 

3«6c,d  Predicted  Post-Shock  Cress  Sections,  0  +  80,  16' 

Diam.,  16'  Rock  Bolts  at  2'  95 

3.7a  Geologic  Cross  Section,  0  +  90,  16’  Diara.,  16' 

Rock  Bolts  at  2'  96 

3.7b  Influential  Joints,  0  +  90,  16'  Diara.,  16'  Rock 

Bolts  at  2’  96 

3*7c,d  Predicted  Post-Shock  Cross  Sections,  0  +  90,  16' 

Diam.,  16'  Rock  Bolts  at  2'  97 

3.8a  Geologic  Cross  Section  1  +  00,  16 '  Diam.,  Unlined  98 

3.8b  Influential  Joints,  1+00,  16'  Diam,,  Unlined  98 

3.8c, d  Predicted  Post-Shock  Cross  Sections,  1  +  00,  16’ 

Diam.,  Unlined  99 

3.9a  Geologic  Cross  Section,  1  +  10,  16*  Diam.,  Unlined  100 

3.9b  Influential  Joints,  1  +  10,  16'  Diam.,  Unlined  100 

3*9c,d  Predicted  Post-Shock  Cross  Sections,  1  +  10,  16' 

Diam.,  Unlined  101 

3.10a  Geologic  Cross  Section,  1  +  20,  18'  Diam.,  Unlined  102 


FIGURES  -  continued 

3. 10b  Influential  Joints,  1  +  20,  18'  Diana.,  Unlined  102 

3.10c,d  Predicted  Post-Shock  Cross  Sections,  1  +  20,  18’ 

Diaa. ,  Unlined  103 

3.11a  Geologic  Cross  Section,  1  +  30,  13'  Diana.,  12* 

Rock  Bolts  at  2*  104 

3.11b  Influential  Joints,  1  +  30,  13'  Diaa.,  12*  Rock 

Bolts  at  2’  104 

3.11c,d  Predicted  Post-Shock  Cross  Sections,  1  +  30,  13* 

Diaa.,  12*  Rock  Bolts  at  2'  105 

3.32a  Geologic  Cross  Section,  1  +  40,  13’  Diaa.,  12’ 

Rock  Bolts  at  2'  106 

3.12b  Influential  Joints,  1  +  40,  13'  Diaa.,  12*  Rock 

Bolts  at  2’  106 

3.12e,d  Predicted  Post-Shock  Cross  Sections,  1  +  40,  13’ 

Diaa.,  12’  Rock  Bolts  at  2'  107 

3.13a  Geologic  Cross  Section,  1  +  60,  13*  Diaa.,  12* 

Rock  Anchors  at  2’  108 

k 

3.13b  Influential  Joints,  1+  60,  13’.Diaa.,  12’  Rock 

Anchors  at  2’  108 

3.13c,d  Predicted  Post-Shock  Cross  Sections,  1  +  60,  13’ 

Diaa.,  12'  Rock  Anchors  at  2'  109 

3.14a  Geologic  Cross  Sections,  1  +  70,  9'  Diaa.,  Unlined  110 

3.14b  Influential  Joints,  1  +  70,  9*  Diaa.,  Unlined  110 

3.l4c,d  Predicted  Post-Shock  Cross  Sections,  1  +  70,  9* 

Diaa.,  Unlined  111 

3.15  Sectional  Elevation  at  Centerline  and  Geology  -  CR 

Drift  112 

3.16  Cross  Sections  and  Geology  -  CR  Drift  113 

3.17a  Geological  Logs  CR  North  Drift  station  1  +  52  to 

1+98  123 

3.17b  Geological  Logs,  CR  North  Drift  Station  1  +  98  bo 

End  124 

3.1B  Joint*  Influence  Diagraas  CR  Drift  -  10,000  pei. 

Blast  125 

3.19a  Geologic  Cross  Section,  1  +  60,  CR  North,  13'  Diaa., 

12*  Rock  Bolts  at  2'  126 

3.19b  Influential  Joints,  1  +  60,  CR  North,  13'  Diaa., 

12’  Rock  Bolts  at  2'  126 

3d9c,d  Predicted  Poet-Shock  Cross  Sections,  1  +  60,  CR 

North,  13''  Diaa.,  12'  Rock  Bolts  at  2'  127 

3.20a  Geologic  Cross  Section,  1  +  70  CR  North,  13' 

Diaa. ,  Unlined  128 


4 


viii 


FIGURES  - 

continued 

Page 

3.20b 

Influential  Joints,  1  +  70  CR  North,  13*  Diam., 
Unlined 

128 

3«20c,d 

Predicted  Post-Shock  Cross  Sections,  1  +  70  CR 

North,  13’  Diam.,  Unlined 

129 

3.21a 

Geologic  Cross  Section,  1  +  80  CR  North,  18’  Diam., 
Unlined 

130 

3.21b 

Influential  Joints,  1  +  80  CR  North,  18’  Diam., 
Unlined 

130 

3.21c, d 

Predicted  Post-Shock  Cross  Sections,  1  +  80  CR 

North;  18*  Diam.,  Unlined 

131 

3.22a 

Geologic  Cross  Section,  1  +  90  CR  North,  18*  Diam., 
Unlined 

132 

3.22b 

Influential  Joints,  1  +  90  CR  North,  18'  Diam., 
Unlined 

132 

3-22c,d 

Predicted  Post-Shock  Cross  Sections,  1  +  90  CR 

North,  18’  Diam.,  Unlined 

133 

3.23a 

Geologic  Cross  Section,  2  +  00  CR  North,  lo'  Diam., 
16'  Rock  Bolts  at  2’ 

134 

3.23b 

Influential  Joints,  2  +  00  CR  North,  16'  Diam., 

16'  Rock  Bolts  at  2' 

134 

3-23c,d 

Predicted  Post-Shock  Cross  Sections,  2  +  00  CR 

North,  16'  Diam.,  16’  Rock  Bolts  at  2' 

135 

3.24a 

Geologic  Cross  Section,  2  +  10  CR  North,  16'  Diam., 
16'  Rock  Bolts  at  2' 

136 

3.24b 

Influential  Joints,  2  +  10  CR  North,  16'  Diam., 

16'  Rock  Bolts  at  2' 

136 

3-24c,d 

Predicted  Post-Shock  Cross  Sections,  2  +  10  CR 

North,  16'  Diam.,  16'  Rock  Bolts  at  2' 

137 

3.25a 

Geologic  Cross  Section,  2  +  2C  CR  North,  16'  Diam., 
lo'  Rock  Bolts  at  2' 

138 

3.25b 

Influential  Joints,  2  +  20  CR  North,  16*  Diam., 

16'  Rock  Bolts  at  2* 

138 

3.25c, d 

Predicted  Post-Shock  Cross  Sections,  2  +  20  CR 

North,  16'  Diam.,  16'  Rock  Bolts  at  2' 

139 

3.26a 

Geologic  Cross  Section,  2  +  30  CR  North,  16'  Diam., 
16*  Rock  Bolts  at  2' 

l4o 

3.26b 

Influential  Joints,  2  +  30  CR  North,  16'  Diam., 

16'  Rock  Bolts  at  2' 

l4o 

3.26c,d 

Predicted  Post-Shock  Cross  Sections,  2  +  30  CR  North 
16 '  Diam.,  16'  Rock  Bolts  at  2' 

’l4l 

3.27a 

Geologic  Cross  Section,  2  +  38  CR  North,  16*  Diam., 
16 '  Rock  Bolts  at  2’ 

142 

FIGURES  -  continued 

3.27b  Influential  Joints,  2  +  38  CR  North,  16'  Dies.,  16' 


Rock  Bolts  at  2*  142 

3. 27c, 4  Predicted  Post-Shock  Cross  Sections,  2  +  38  CR 

North,  16*  Diaa.,  16'  Rock  Bolts  at  2*  143 

3*28  Sectional  Elevation  at  Centerline  and  Geology  -  CR 

North  Drift  144 

3.29  Cross  Sections  and  Geology  -  CR  North  Drift  145 

3*30  Geological  Logs  DL  1,  2  Station  0  +  50  to  0  +  90  153 

3.31  Influence  Diagrams  for  DL  1  2  Drift,  10,000  pel. 

Blast  Pressure  154 

3.32a  Geologic  Cross  Section  -  Sta.  0  +  50  155 

3.32b  Influential  Joints  -  Sta.  0  +  50  155 

3.32c  Predicted  Po6t-Shock  Section  -  Sta.  0  +  50  156 

3* 33*  Geologic  Cross  Section  -  Sta.  0  +  60  157 

3.33b  Influential  Joints  -  Sta.  0  +  60  157 

3.33c  Predicted  Post-Shock  section  -  Sta.  0  +  60  158 

3. 34a  Geologic  Cross  Section  -  Sta.  0  +  70  159 

3*  3 to  Influential  Joints  -  Sta.  0  +  70  159 

3.3to  Predicted  Post-Shock  Section  -  Sta.  0  +  70  l60 

3* 35a  Geologic  Cross  Section  -  Sta.  0  +  80  l6l 

3.35b  Influential  Joints  -  Sta.  0  +  80  l6l 

3.35c  Predicted  Post  Shock  Section  -  Sta.  0  +  80  162 

3.36  Sectional  Elevation  at  Centerline  and  Geology  - 

DL  Drift  I63 

3*37  Cross  Sections  and  Geology  -  DL  Drift  164 


x 


EFFECT  OF  JOINTS  ON  THE  STRENGTH  OF  TUNNELS 


1.  INTRODUCTION 


This  is  a  report  of  an  investigation  performed  under  Contract  Number 
DACA45-67-C-0015  Mod.  P001,  between  the  Corps  of  Engineers,  Omaha  District, 
and  Dr.  Richard  E.  Goodman,  as  contra 'tor.  It  describes  results  of  re¬ 
search  on  an  application  of  analysis  techniques  to  the  Piledriver  Drift. 

The  work  was  performed  during  the  period,  July  1,  1967  to  June  30,  1968. 
This  work  was  specifically  tied  to  the  Piledriver  Project.  An  effort  was 
made  to  apply  previously  developed  analytical  techniques  to  evaluate  the 
influence  of  jointing  and  faulting  in  the  bedrock  on  the  performance  of 
Piledriver  test  drifts.  Frequent  reference  is  made  to  the  final  reports  of 
two  previous  contracts^-* 

In  the  first  contract,  the  principal  investigators  studied  the  effect 
of  geological  factors  on  the  behavior  of  several  experimental  lined  and 
unlined  sections  of  the  Piledriver  test.  Particular  attention  was  paid  to 
the  effect  of  joint  planes  on  the  strength  of  a  tunnel  under  static  load. 
Consideration  was  also  given  to  the  strengthening  effect  of  rock  bolts  and 
of  composite  steel-concrete  linings.  The  first  contract  was  a  specifically 
directed  investigation  of  short  duration.  In  the  second  contract,  several 
ideas  for  analysis  of  geological  factors  which  had  been  considered  in  the 
first  contract  were  given  further  attention  and  the  importance  of  weakness 
surfaces  on  the  mechanics  of  tunnels  was  underscored.  It  was  stated  that 
the  technology  of  description  and  prediction  of  joints  and  fractures  was 
sufficiently  far  advanced  to  permit  analysis  but  that  the  mechanics  of  the 


^Goodman,  R.  E.,  Geological  Factors  in  Design  of  Blast  Resistant 
Tunnels ,  U.S.  Army,  Corps  of  Engineers,  Omaha  District,  Tech.  Report  No.  2, 
September  1966. 

^Goodman,  R.  E.,  Analysis  of  Structures  in  Jointed  Rock,  U.  S.  Army, 
Corps  of  Engineers,  Omaha  District,  Tech.  Report  No.  3,  September  1967. 


2 


interaction  of  joint  blocks  around  an  opening  had  not  been  considered  in 
any  depth. 

Two  methods  of  analysis  were  proposed  and  developed  to  some  extent  in 
the  second  report.  The  first  was  a  three-dimensional  ubiquitous  joint 
analysis ,  referred  to  also  by  the  more  simple  title  of  joint  influence 
analysis .  It  involves  the  calculation  of  the  area  around  the  tunnel  within 
which  the  traction  across  a  joint  of  fixed  orientation  and  given  cohesion 
and  friction  is  sufficient  to  induce  sliding  along  or  opening  of  the  joint. 

The  second  method  of  analysis  discussed  in  the  second  contract  report 
was  a  joint  stiffness  analysis  using  the  method  of  finite  element  analysis. 

In  this  method,  an  attempt  is  made  to  compile  a  mathematical  model  of  a 
particular  cross  section  of  a  tunnel,  and  to  introduce  into  the  model, 
representations  for  the  deformability  and  strength  of  both  the  solid  rock 
blocks  and  the  joints  and  fractures  between  them. 

In  this  report,  the  two  above  methods  have  been  applied  to  analysis 
of  selected  portions  of  the  Piledriver  experiment.  The  Piledriver  experi¬ 
ment  includes  test  galleries  of  various  distances  from  a  blast  as  depicted 
in  Fig.  1.1.  In  this  study,  the  CR  and  DL  Drifts  shown  in  Fig.  1.1  were 
considered.  These  drifts  include  unlined  sections,  rock-bolted  sections, 
and  lined  sections  with  various  combinations  of  concrete,  steel  and  compres¬ 
sible  back-packing.  Detailed  descriptions  of  the  test  openings  are  given 
in  a  number  of  other  reports  and  are  summarized  in  Technical  Report  Number 

23. 

The  principal  investigator  studied  general  descriptive  information 
about  the  Piledriver  test,  including  geological  maps  and  logs  compiled 
before  the  experiment,  lie  made  two  trips  —  one  with  Dr.  Robert  Taylor 
and  one  alone  —  to  examine  the  failure  pattern  and  damage  patterns  of  drifts 
after  the  blast.  He  also  studied  post-blast  cross  sections  and  profiles, 
portions  of  the  Corps  of  Engineers  Report  describing  damage  in  the  drifts^, 
and  post-shock  geological  notes  made  by  the  geologist  on  the  pre-construction 
logs.  Work  was  performed  at  131  Montrose  Road,  Berkeley,  California,  and 
at  the  University  of  California  Computer  Center.  One  progress  report  was 


3op.  cit.,  p.  1. 

^Chapter  4,  Piledriver  Report  POR  4015. 
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sent  on  January  12,  describing  joint  influence  analysis  in  DL  drift.  This 
document  includes  all  the  information  contained  in  that  progress  report. 

This  report  is  organized  as  follows.  Section  2  describes  the  method 
of  mathematical  modeling  for  tunnels  in  jointed  rock  including  modifications 
made  to  finite  element  analysis  methods  used  herein.  Methods  of  developing 
geologic  cross  sections  from  construction  logs  are  described  and  there  is 
a  discussion  of  the  properties  to  be  assigned  to  the  rock  blocks,  including 
deformability  and  strength.  A  new  criterion  of  failure  for  Piledriver  rock 
was  incorporated  in  the  study  and  is  described  in  this  section.  The  pro¬ 
perties  of  joint  surfaces  are  discussed  and  the  methods  of  sampling  and 
testing  joints  to  determine  their  deformability  and  strength  properties 
are  reviewed.  Typical  shear-displacement  curves  are  presented  for  the 
various  kinds  of  joint  surfaces  that  are  known  to  occur.  Extensive  data 
on  joint  strength  and  deformability  properties  are  presented  in  a  series 
of  tables  compiled  from  many  sources  around  the  world.  Finally,  at  the 
end  of  this  section,  the  application  of  the  modeling  technique  to  the  DL 
drift  is  discussed. 

Section  3  describes  the  use  of  the  joint  influence  diagram  method  to 
DL  and  CR  drifts.  The  method  of  constructing  and  using  joint  influence  dia¬ 
grams  is  reviewed,  and  the  state  of  stress  around  each  of  the  tunnels  adopted 
for  the  analysis  is  presented.  Analyses  were  made  every  10  feet  throughout 
CR  lateral  and  CR  north  drift  and  throughout  DL  1  and  DL  2.  At  every  cross 
section  considered,  the  results  of  analysis  are  presented  in  a  series  of 
figures  and  compared  with  the  actual  damage  patterns  revealed  after  the  blast. 
Computer  programs  used  in  this  work  are  presented  in  the  Appendix. 
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2.  MATHEMATICAL  MODELING  OF  TUNNELS  IN  JOINTED  ROCK 


2.1  GENERAL  -  FINITE  ELEMENT  METHOD 

The  finite  element  method  was  used  to  construct  mathematical  models 
of  selected  cross  sections  in  C  and  D  drifts.  The  methods  used  were  des¬ 
cribed  in  the  two  previous  reports**.  Several  modifications  and  additions 
were  introduced  in  the  course  of  this  investigation,  to  allow  for  incremental 
loading  of  large  systems  of  rock  blocks  and  joints  under  plane  strain  condi¬ 
tions  with  large  displacements.  The  computer  program  used  in  the  analyses 
is  included  in  Appendix  1  together  with  explanatory  notes  on  its  operation. 

The  results  of  computations  can  be  no  better  than  the  quality  of  the 
input  data:  ’These  are:  the  structural  geology  of  the  cross  sections;  the 
properties  of  the  rocks;  and  the  properties  of  the  joints. 


2.2  GEOLOGICAL  DATA 

The  general  geology  of  the  Piledriver  drifts  is  well  known  from  geo¬ 
logical  logs  of  the  walls  and  roof  of  the  drifts  prepared  by  Corps  of  Engineers 
geologists,  J.  Zeltinger  and  Harold  Jack.  Fig.  2.1  presents,  for  example, 
the  logs  for  DL  1  and  2.  These  logs  provide  the  information  for  preparation 
of  detailed  cross  sections  showing  the  outlines  of  joint  bounded  blocks 
around  the  tunnel.  Fig.  2.2  shows  how  these  detailed  sections  are  prepared 
by  projection  and  extrapolation.  Thus,  cross  sections  of  the  drifts  can 
be  constructed  as  in  Fig.  2,3. 

The  detailed  geologic  sections  present  data  on  the  distribution  of 
joints  and  defects  and  the  shape  of  rock  blocks  behind  the  tunnel  walls. 

From  these  data,  a  finite  element  computing  mesh  can  be  compiled.  In  the 
subsequent  analysis,  obviously  the  shapes  of  the  blocks  delimited  by  the  net¬ 
work  of  discontinuities  is  very  important.  The  confidence  with  which  this 
is  known  is  greatest  closest  to  the  wall  of  the  tunnel  where  the  extrapolation 
of  geological  data  required  is  minimum,  fortunately,  this  is  the  most  important 
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region  for  the  analysis.  Unfortunately,  joints  which  parallel  the  tunnel 
axis  behind  the  wall  may  not  be  seen  in  the  tunnel  and  the  only  way  they 
can  be  placed  on  the  cross  sections  Is  by  regularly  repeating  joints  of 
the  same  orientation  that  do  intersect  the  tunnel.  However,  a  slight  in¬ 
accuracy  in  the  placement  of  a  single  joint  will  lead  to  fictitious  slivers 
and  wedges  of  rock  which  may  tend  to  fail  by  bending  and  crushing  as  the 
joint  blocks  move;  thus,  it  is  important  to  define  the  correct  cross  section 
as  accurately  as  possible.  Use  of  a  bore  hole  stratascope,  camera,  or  tele¬ 
vision  device  would  add  a  new  dimension  of  precision  to  the  development  of 
geologic  cross  sections  and,  therefore,  is  very  strongly  recommended  for 
studies  of*  this  type. 

The  accuracy  of  the  geologic  projections  can  be  checked  in  transition 
sections  where  tunnels  go  from  one  size  to  another,  in  end  sections  where 
a  drift  terminates  and  in  sharp  bends  wher^  drifts  are  passing  through  joints 
previously  seen  around  the  bend.  The  quality  of  projection  in  the  example 
gLven  for  DL  1  and  2  is  presented  in  Fig.  2.4,  where  the  projected  end 
section  can  be  compared  with  the  actually  logged  section.  The  general 
treads  are  correct  but  there  is  much  inaccuracy  in  the  configurations  of 
individual  blocks.  The  primary  difference  is  the  offsetting  of  the  D  joints 
by  the  A  shears,  a  factor  not  considered  in  making  these  projections. 

Greater  precision  is  achieved  close  to  the  walls  than  in  the  center  of  the 
end  plate  or  far  behind  the  walls. 


2.3  ROCK  PROPERTIES 

The  properties  of  the  Piledriver  rock,  previously  described6,  were 
determined  by  the  Missouri  River  Division  Laboratory  by  direct  pull  tension, 
and  unconfined  and  triaxial  compression  tests.  The  intact  rock  specimens 
acted  as  brittle,  linear  elastic  solids  at  confining  pressures  as  high  as 
3000  psi.  The  tensile  strength  was  1450  psi  and  the  average  unconfined 
compressive  strength  was  30,500  psi.  The  modulus  of  elasticity  was  11.0  to 
12.0  >'•  10^  psi  in  compression  and  8.0  x  106  psi  in  tension. 


^Technical  Report  Number  2  (op.  cit.,  p.  1) 
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The  criterion  of  failure  adopted  for  previous  work  based  on  the  test 
data  was  adapted  from  the  Griffith  and  Modified  Griffith  theories  of  failure. 
This  is  a  two  dimensional  approach  in  which  the  intermediate  principal  stress 
is  ignored.  A  new  approach  is  introduced  herein  where  the  three  principal 
stresses  are  all  considered  by  using  a  power  law  in  terms  of  octahedral 
stresses. 

(1)  Power  Law  Failure  Criterion 

The  criterion  of  failure  of  rocks  may  be  expressed  in  terms  of  three 
principal  stresses,  or  invariants  formed  from  them.  In  a  three  axis  prin¬ 
cipal  stress  space,  the  criterion  of  failure  would  delimit  a  closed  solid 
surface  within  which  the  coordinates  of  any  point  define  a  stable  combina¬ 
tion  of  principal  stresses.  In  detail,  it  is  not  known  what  the  shape  of 
the  failure  surface  is  for  the  various  classes  of  rocks.  However,  all  rocks 
have  a  relatively  small  tensile  strength  and  show  an  increase  in  compressive 
strength  with  confining  pressure  at  a  rate  which  decreases  with  higher 
confinements.  Jaeger^  has  suggested  that  the  failure  surface  for  rocks 
may  be  a  solid  of  revolution  about  the  line  making  equal  angles  with  the 
stress  axes  (the  hydrostatic  axis) .  He  suggested  a  paraboloid  of  revolution 
obtained  by  rotating  the  Griffith  Parabola  about  the  hydrostatic  axis. 

Any  solid  of  revolution  about  the  hydrostatic  axis  can  be  expressed  as  a 
curve  in  a  plane  defined  by  octahedral  shearing  stress  (toct)  as  ordinate 
and  octahedral  normal  stress  (0oct)  as  abscissa;  these  quantities  are 
expressible  in  terms  of  principal  stresses  as  follows: 

T0ct  *  3  Aoi  -  02)  2  +  (o 2  -  a3)  2  +  (a3  -  0^)2  (2-1) 

and 

°oct  "  f  (°1  +  <?2  +  o3)  (2-2) 

Rather  than  constraining  the  failure  law  to  conform  to  a  specific 
curve,  e.g.  parabaloid  or  cone,  a  more  general  power  law  relation  between 
Toct  and  0oct  has  been  adopted,  i.e. 


^Proceedings  8th  Symposium  on  Rock  Mechanics,  AIME,  1966,  p.  3-57/ 
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Toct  "  N  +  D  <JoctB  (2-3) 

This  is  simply  an  empirical  criterion  of  failure  conforming  to  the  constraint 
that  it  represent  a  solid  of  revolution  about  the  hydrostatic  axis  in  prin¬ 
cipal  stress  space.  To  apply  this  criterion,  Eq.  2-3  shall  also  be  written 
in  logarithmic  form,  as  follows: 

l°8l0  (Toct  “  =  A  +  B  log1Q  Goct  (2-4) 

(tfoct  >  0) 

where  A  has  been  written  in  place  of  log^Q  D.  D,  B,  and  N  are  strength 
constants  of  the  material.  The  dimensionless  constant  B  describes  the  cur¬ 
vature  of  the  failure  surface.  D  and  N,  whose  values  depend  on  the  system 
of  units  employed,  describe  something  like  the  friction  and  cohesion.  N 
is  the  value  of  the  octahedral  shear  stress  when  the  octahedral  normal  stress 
is  0.  It  could  be  determined  directly  from  a  pure  shear  test  on  a  plate  of 
rock  (o^  *  -Oy  =  0) ,  or  from  a  torsion  test  of  a  solid  cylinder  of  rock. 

N  can  be  estimated  as  the  intercept  on  the  Toct  axis  of  the  failure  curve, 
plotted  from  test  data  in  the  QQCt  ~  Toct  plane.  Another  method  of  estimating 
N  would  be  to  calculate  N  from  a  formula  derived  for  the  special  case  of 
Mohr-Coulomb  theory. 

A  three-dimensional  generalization  of  the  Mohr-Coulomb  theory  is  repre¬ 
sented  by  Eq.  2-3  if  B  is  equal  to  1.  From  triaxial  test  data,  one  may 
calculate  a  Mohr-envelope  and  determine  the  constants  <j>  (the  angle  of  internal 
friction),  and  c  (the  shear  strength  intercept  or  cohesion) .  Toct  and  °oct 
may  also  be  calculated  from  triaxial  compression  test  results,  as  follows 
(compare  with  Eqs.  2-1,  2-2) 

Toct  *  — J  (°lf  ”  ?)  (2-5) 

^oct  =  °lf  +  2P  (2-6) 

3 

where  p  equals  the  confining  pressure  and  *  the  axial  pressure  at  the 
time  of  failure.  Substituting  these  expressions  in  Eq.  2-3  yields  the 
expression 
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N  + 


+ 
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For  the  case  where  B  »  1,  Eq.  2-7  simplifies  to 


3N 

>7  -  D 


+ 


2D  + 
vT  -  D 


P 


(2-7) 


(2-8) 


The  analogous  formula  for  the  Mohr-Coulomb  theory  is: 

*  2c  tana  +  tan2a  p  (2-9) 

where  a  *  (45  +  <f>/ 2) . 

By  comparing  Eqs.  2-8  and  2-9,  one  may  write  the  two  identities: 


2  c  tana 


(2-10) 


2D  +  /2 
/2  -  D 


tan2a 


(2-11) 


which  may  be  solved  for  N-  to  yield 


=  2/2  c  tan(45  +  4>/2) 

2  +  tan2(45  +  <j>/2) 


(2-12) 


Eq.  2-12  may  be  used  to  gain  an  estimate  of  N  in  order  to  determine  the 
constants  A  and  B  from  actual  failure  data.  One  may  also  solve  Eqs.  2-10 
and  2-11  fr'r‘  D  to  yield: 


/2  [ tan2(45  +  <j>/2)  -  1] 

2  +  tan2(45  +  (j»/2) 


(2-13) 


Eq.  2-13,  together  with  Eq.  2-12,  provides  values  for  the  strength  constants 
D  and  N,  with  B  »  1,  representing  an  extension  of  the  Mohr-Coulomb  theory  to 
three  dimensions  by  rotation  about  the  hydrostatic  axis  to  form  a  cone. 
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However,  it  would  be  more  sensible  to  determine  the  best  values  of  the 
strength  constants  by  empirical  means. 

The  procedures  for  determining  the  values  of  D,  B,  and  N  to  fit  actual 
test  data  are  as  follows: 

1.  Convert  the  breaking  stresses,  in  laboratory  or  field  strength 
tests,  to  octahedral  normal  and  shear  stresses  using  Eqs.  2-1 
and  2-2.  In  the  triaxial  compression  test,  for  example,  with 
confining  pressure  p  and  axial  stress  at  failure  Oj_  ,  (a^  = 

and  ^2  =  °3  =  P) 

(°lf  -  p) 

aoct  "  - 3 -  +  P  (2-14) 


4  <«if  -  p) 


(2-15) 


while  in  the  uniaxial  tension  test,  with  tensile  strength  =  -ot, 


(2-16) 


(2-17) 


2.  Determine  N.  It  may  be  computed  from  Eq.  2-12  (<{>  and  c  must  be 
determined  first),  or,  if  data  exist  in  the  region  of  low  octa¬ 
hedral  normal  stresses,  e.g.  tension  or  torsion  or  simple  shear 
test,  it  may  be  determined  graphically  by  plotting  the  values  of 
aoct  anc*  Toct  at  failure  on  arithmetic  graph  paper  and  sketching 
a  smooth  curve  to  find  the  value  of  Toct  when  <?oct  =  0.  If  data 
are  available  only  from  confined  and  unconfined  compression  tests, 
the  latter  method  will  not  be  sufficiently  accurate. 

3.  determine  D  and  B  by  plotting  the  data  to  log-log  scale.  The 
ordinate  is  log-^Q  (Toct  “  N)  and  the  abscissa  is  log^Q  (aoct) • 

The  best  fit  straight  line  yields  the  constants  A  and  B  (Eq.  2-4); 
then,  D  =  log-^A. 
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In  employing  the  graphical  procedure  above,  one  must  be  careful  to 
avoid  prejudicing  the  data  because  of  scatter  in  repeated  unconfined  com¬ 
pression  or  tension  tests.  The  octahedral  stresses  calculated  from  any  set 
of  unconfined  compression  tests  must  lie  along  the  straight  line  Toct  * 

Jl  ooct.  Similarly,  scattered  data  from  repeated  uniaxial  tension  tests  be 
along  the  line  TQCt  *  -/l  Coct.  To  avoid  prejudicing  the  data,  for  results 
of  unconfined  compression  or  tension  tests,  only  one  point,  the  one  corres¬ 
ponding  to  the  average  strength,  should  be  converted  to  octahedral  stresses 
and  graphed. 

Figs.  2.5  and  2.6-  show  strength  data  plotted  in  this  way  for  Piledrlver 
rock.  In  Fig.  2.5,  the  arithmetic  plot  of  xoct  versus  aoct  based  on  com¬ 
pression  tests  is  curved  slightly  downward.  The  uniaxial  tension  strength 
of  -1450  psi  yields  the  point  Toct  *  +685  psi,  0Qct  *  -484  p«i.  The  inter¬ 
polated  value  of  N  is  1300  psi.  Eq.  2-12,  with  <j>  *  56°  and  c  *  3600  psi, 
gives  N  =  2,540  psi.  In  Fig.  2-6,  log-^Q  (tocC  -  N)  is  plotted  against 
log1Q  (ooct)  for  both  values  of  N  and  for  N  *  0.  The  set  of  points  corres¬ 
ponding  to  N  =  1300  psi  provide  the  least  scatter  about  a  straight  line  fit. 
The  rock  strength  constants  for  the  three  cases  are  as  follows: 

(1)  (2)  (3) 

N  -  0  N  -  1300  psi  N  *  2450  psi 

D  =  6.60  D  -  2.95  D  -  1.14 

B  =  0.88  B  *  0.91  B  -  1.00 

The  second  case  was  adopted  as  the  failure  law  for  the  Piledriver  rock,  i.e. 

Toct  =  1300  +  2*95  °oct0*91  (2-3.8) 

For  a  parabaloid  (Griffith  theory  in  three  dimensions),  B  *  0.5,  while 
B  *  1.0  for  a  cone  (a  generalization  of  Modified  Griffith  and  Mohr-Coulowb 

O 

theories).  Franklin0  reviewed  data  from  many  triaxial  tests  and  found 
that  most  data  fall  between  the  modified  Griffith  and  Griffith  theories, 
meaning  that  the  failure  surface  curves  downward  slightly  as  confining  pres¬ 
sure  is  increased,  but  in  general,  not  so  much  as  predicted  by  the  Griffith 

^Franklin,  J.  A.,  A  Strength  Criterion  for  Rock,  Imperial  College 
Rock  Mechanics  Research  Project  Number  6,  1968. 
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theory.  The  Piledriver  quartz  monzonite  strength  data  fall  within  the 
observed  range, 

(2)  Failure  Law  When  0oct  Is  Negative 

When  the  octahedral  normal  stress  becomes  tensile,  mathematical  dif¬ 
ficulty  is  experienced  in  working  with  Eq.  2-3.  Hoek^  suggested  that  the 
controlling  strength  parameter  is  simply  the  uniaxial  tension  test  in  this 
region,  i.e.  the  rock,  breaks  when  the  minor  principal  stress  is  less  than 
the  uniaxial  tensile  strength.  This  criterion  has  been  adopted  herein. 
Another  possibility  would  be  to  adopt  the  maximum  extension  strain  criterion 
advanced  by  Trollope^  in  which  the  greatest  extension  strain,  in  any  com¬ 
bined  stress  state,  is  unsafe  if  it  exceeds  the  strain  at  failure  in  a  uni¬ 
axial  tension  test.  It  would  be  a  simple  matter  to  introduce  this  criterion 
in  the  computations. 

The  properties  adopted  for  the  rock  are  summarized  in  Table  2.1. 


2. A  JOINT  PROPERTIES 

(1)  Importance  of  Joints 

Much  has  been  written  in  recent  years  about  the  significance  of  planes 
of  weakness  on  the  stability  of  structures  and  excavations  in  hard  rock. 

It  is  quite  possible  that  a  person  reviewing  the  recent  publications  and 
reports  in  the  field  of  rock  mechanics  might  become  wearied  by  all  of  these 
references  to  fractures,  discontinuities,  faults  and  their  impact,  he  might 
feel  that  the  problem  has  been  overstated.  Nothing  could  be  farther  from 
the  truth.  In  an  excavation  in  hard  rock,  such  as  the  granitic  rock  of  the 
Piledriver  Drift,  if  joints  and  defects  of  the  rock  mass  are  not  considered, 
it  is  very  hard  to  demonstrate  why  a  drift  should  fail,  even  if  it  is  un¬ 
lined.  Post-shock  evaluation  of  the  test  sections  at  intermediate  range 


^Hoek,  E.,  Rock  Fracture  Under  Static  Stress  Conditions,  CSIR  Report 
MEG  383,  October  1965. 

■^Duncan,  J,  M.  and  Goodman,  R.  E.,  Methods  of  Analysis  for  Rock  Slopes, 
U.  S.  Arajr  Engineers,  Waterways  Experiment  Station,  1968,  p.  45. 
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reveals  conclusively  that  the  relative  movement  of  joint-bounded  blocks  waa 
the  most  important  single  criterion  in  causing  damage  to  lined  and  unlined 
sections. 

To  demonstrate  the  importance  of  joints  in  a  tunnel  in  hard  rock,  the 
computer  methods  used  herein  were  applied  to  the  study  of  a  circular  tunnel 
section  in  homogeneous  hard  rock  with  the  properties  of  the  Piledriver 
quartz  raonzonite.  (These  properties  are  listed  in  Table  2.1.)  Fig.  2.7 
shows  the  finite  element  mesh  (PD-5)  used  for  this  analysis.  A  static  pres¬ 
sure  was  applied  horizontally  in  increments  of  4000  psi,  reaching  a  peak 
pressure  of  32,000  psi.  Simultaneously,  vertical  compression  of  0.43  *  the 
blast  pressure  was  applied  to  represent  the  confinement  in  the  wave  front. 

With  these  boundary  conditions,  there  is  no  tension  developed,  theoretically, 
at  any  point  around  the  opening,  but  very  significant  compressive  stresses 
developed  in 'Che  tangential  direction  at  the  roof  and  floor.  No  element 
in  the  finite  element  mesh  experienced  failure  at  any  stage  in  the  incre¬ 
mental  loading  process,  even  when  32,000  psi  had  been  reached.  Unquestionably, 
the  wall  rock  at  the  top  and  bottom  of  the  opening  would  have  broken.  How¬ 
ever,  the  surface  layer  of  finite  elements  cannot  represent  the  stresses 
on  the  surface  but  rather  at  the  element  centers,  which  for  this  mesh  were 
approximately  1/32  of  a  radius  behind  the  wall.  Even  though  this  is  a  short 
distance  behind  the  wall,  the  stress  gradient  is  so  steep  that  the  stresses 
had  already  reached  tolerable  values  and  no  failure  was  predicted  by  the 
failure  criterion  adopted. 

In  order  to  examine  a  situation  with  more  severe  loading  conditions  a 
second  case  was  investigated  where  the  ratio  of  free  field  pressures  was 
0.25,  i.e.  the  vertical  pressure  was  25%  of  the  blast  pressure.  This  would 
correspond  to  a  value  of  Poisson's  ratio  of  0.2  and  is  more  drastic  than 
the  actual  wave  motion  believed  to  have  been  experienced  in  the  Piledriver 
Drifts.  At  this  confining  ratio,  some  tension  develops  in  the  0°  and  180° 
walls.  The  progress  of  failure  with  increasing  increments  of  applied  pres¬ 
sure  is  plotted  in  Fig.  2  .8.  With  the  first  increment,  4,000  psi  blast 
pressure,  no  elements  had  failed;  at  8,000  psi,  two  elements  had  failed  in 
the  tensile  region  (Fig,  2,8a).  The  volume  of  rock  which  was  actually  dis¬ 
turbed  by  this  very  significant  stress  is  not  great  and  it  is  quite  possible 
that  had  there  been  a  rock  bolt  reinforcement  or  steel  lining,  nothing  very 


severe  in  the  way  of  change  would  have  occurred  in  the  tunnel.  By  the 
time  pressure  of  32,000  psi  had  been  reached,  damage  had  extended  to  only 
five  elements,  encompassing  only  those  elements  near  the  horizontal  diameter 
adjacent  to  or  j.ust  behind  the  wall  (Fig.  2,8c).  No  elements  are  shown 
as  failed  at  90°  (roof)  as  the  stresses  at  element  centers  were  tolerable, 
even  at  p^  =  32,000  psi. 

While  the  type  of  failure  pattern  described  above  jnay  be  relevant  to 
the  case  of  an  unreinforced  hole  of  small  diameter  for  the  installation  of 
telephone  conduit,  it  has  very  little  to  do  with  the  case  of  lined  or  bolted 
16  foot  diameter  opening  in  rock.  The  localized  crushing  and  slabbing  of 
rock  which  Was  observed  underground  in  some  of  the  drifts  must  represent  the 
stress  concentration  effect  of  geological  weakness  planes  on  the  inter¬ 
vening  blocks.  In  order  to  make  an  analysis  of  this  factor,  it  is  necessary 
to  gain  samples  of  the  weakness  planes,  to  t$st  them,  to  determine  their 
strength  and  stiffness,  and  to  construct  a  finite  element  mesh.  The  subject 
of  sampling  and  testing  natural  joint  surfaces  is  not  well  known  in  the 
United  States  but  has  been  advanced  _n  the  last  two  years  to  a  significant 
level  of  achievement  in  Western  Europe.  The  following  discussion  relates 
to  some  of  the  findings  of  these  investigations. 

(2)  Methods  of  Sampling  Joints 

In  order  to  determine  theTmechanical  properties  of  natural  joint  surfaces, 
-it  is  necessary  to  produce  samples  of  significant  size  in  the  laboratory  for 
shear  testing.  In  some  cases,  the  natural  joint  surfaces  are  so  regular* 
smooth  and  planar,  that  the  creation  of  artificial  joints  by  sawing  of  other¬ 
wise  intact  cores  provides  realistic  specimens  for  analysis.  Another  method 
is  by  breaking  beams  irj,  the  field  or  in  the  laboratory  and  testing  the  arti¬ 
ficial  shear  or  extension  features  as  developed.  One  may  also  sample  joints 
by  selecting  blocks  in  an  outcrop  which  are  mated,  carefully  sculpting  and 
removing  them  and  fitting  them  into  a  laboratory  shear  machine.  It  is  also 
possible  to  drill  across  joints  in  the  field,  either  normal  to  the  joints, 
parallel  to  the  joints,  or  at  a  pre-set  angle  to  the  joints.  Finally,  it  is 
possible  to  remove  blocks  containing  joints  by  drilling,  trepanning  or  even 
wire  sawing  in  the  field.  All  of  these  methods  will  be  discussed. 
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Artificial  Joints.  A  large  number  of  friction  tests  have  been  Con¬ 
ducted  in  laboratories  by  sawing  intact  cores  of  rock.  The  two  pieces  may 
then  be  mounted  in  a  direct  shear  box  by  setting  in  Portland  cement  or 
epoxy  resin,  or  fitted  inside  a  triaxial  device.  In  the  latter  case,  the 
diamond  saw  cut  must  be  at  about  45-60°  to  the  ends  of  the  core.  Tests  of 
this  type  were  conducted  by  the  Corps  of  Engineers  in  the  Missouri  River  Div. 
Laboratory.  They  will  be  discussed  in  the  next  section.  It  is  doubtful 
that  an  unmodified  diamond  saw  cut  is  a  very  good  model  of  a  natural  joint 
surface.  Joints  owe  their  origin  to  breakage  of  rock  under  extension  or 
shear.  Extensional  joints  are  characteristically  rough.  While  shear  features 
are  smooth,  they  invariably  contain  a  filling  of  gouge,  clay,  or  crushed 
rock.  A  diamond  saw  specimen  can  model  a  natural  joint  if  a  gouge-like 
layer  is  developed  on  the  diamond  cut  surface.  Research  in  this  area  is 
being  pursued  at  the  University  of  Illinois  by  Mr.  J.  Coulson  under  the 
direction  of  Professor  Donald  Deere.  Diamond  saw  specimens  are  polished  to 
specified  smoothness  and  then  sheared  at  high  normal  pressure.  This  pro¬ 
duces  slickensides  and  crushed  material  on  the  joint  surface  which  resembles 
slickensided  surfaces  sometimes  seen  on  naturally  occurring  shear  joints. 

The  repeated  cycles  of  loading  of  such  specimens  are  then  reasonable  models 
of  the  actual  in  situ  mechanical  characteristics  of  the  joint  surface. 

Fig.  2*9  shows  a  typical  specimen  of  an  artificial  joint  in  the  University 
of  Illinois  test  series. 

Dr.  M.  DeFreitas  and  Dr.  John  Knill  in  the  Department  of  Geology  at 
Imperial  College,  London,  have  been  studying  the  shear  behavior  of  arti¬ 
ficially  induced  extension  joints  in  granite.  Rock  blocks  containing  joints 
are  brought  to  the  laboratory  from  the  quarry  where  they  have  been  produced 
in  hard  rock  by  impact  of  simple  beams.  The  matching  rock  pieces  are  laid 
on  the  laboratory  table  and  a  mold  for  the  shear  box  is  accurately  positioned 
(Fig.  2.10a).  The  joint  specimen  is  cemented  into  the  mold  with  epoxy  cement 
and  then  the  mold  is  removed  and  the  specimen  is  placed  in  the  direct  shear 
machine  for  testing.  The  joints  produced  and  tested  in  this  manner  are  ty¬ 
pically  rough  extension  joints,  with  a  waviness  of  several  centimeters  (Fig. 
2.10b).  Dr.  DeFreitas  has  been  studying  the  effect  of  joint  roughness  on 
the  shear  characteristics  and  in  fact  has  been  measuring  the  shape  of  the 
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surface  with  a  wheel  and  stylus.  Joints  produced  in  this  way  may  be  very 
realistic  models  of  actually  occurring  joints  of  extension  origin. 

Artificial  joints  may  also  be  produced  in  the  testing  machine  by  inci¬ 
pient  shear  failure  of  intact  specimens.  Both  triaxial  and  direct  shear 
specimens  can  be  used.  Dr.  Jaeger,  at  Australian  National  University,  has 
experimented  along  these  lines. 

Field  Sampling  of  Open  Joints.  Very  often,  joint  surfaces  of  interest 
in  the  field  are  loose  or  open  so  that  one  block  rests  on  another  with  per¬ 
fect  fit  but  without  any  cohesion.  In  this  case,  it  is  possible,  with  a 
geology  pick  or  crowbar,  to  remove  blocks  containing  joints.  It  may  not  be 
necessary  to  hold  the  upper  and  lower  blocks  together  for  shipment,  but  only 
to  mark  their  relative  positions  so  that  they  may  be  replaced  in  proper  orien¬ 
tation  in  the  laboratory.  Professor  Krsmanovich  in  Sarajevo,  Yugoslavia, 
has  sampled  a. large  number  of  joints  in  limestone  by  this  method.  The  samples 
are,  of  course,  disturbed  in  the  sense  that  any  original  closing  of  aperture 
tangentially  or  normally  which  may  occur  during  the  actual  loading  in  the 
field  cannot  be  reproduced  in  the  laboratory  specimen. 

Drilling  Joint  Samples.  It  is  possible  to  obtain  a  relatively  undis¬ 
turbed  sample  of  a  joint  by  drilling  through  it  with  a  core  barrel.  There 
are  three  methods  of  doing  this  —  normal  drilling,  parallel  drilling,  and 
inclined  drilling.  The  firm  Coyne  and  Bellier,  in  Paris,  France,  has  ob¬ 
tained  samples  at  several  dam  sites  by  drilling  perpendicular  to  joint 
planes,  so  that  the  joint  occupies  a  position  parallel  to  the  ends  of  the 
core;  after  trimming  in  the  laboratory  and  plotting  in  Portland  cement,  the 
circular  cross-sectioned  joint  surface  is  tested  in  a  direct  shear  machine. 
Samples  from  two  inches  diameter  to  as  much  as  nine  inches  in  diameter 
have  been  obtained  in  this  manner  by  Coyne  and  Bellier  (Fig.  2.11).  It  is 
possible  that  drilling  across  the  joint  by  a  rotary  drill  will  seriously 
disturb  the  joint  surface  as  one  block  turns  on  the  other.  In  fact,  it 
might  be  impossible  to  sample  very  weak  joints  or  open  joints  by  this  tech¬ 
nique. 

A  modification  of  the  perpendicular  drilling  technique  is  being  inves¬ 
tigated  by  Goodman  and  Mahtab  at  present  in  a  project  at  the  University  of 
California,  Berkeley.  In  this  method,  a  small  hole  is  first  drilled  across 
the  joint  and  a  rock  anchor  is  installed  and  tightened;  then  the  small  hole 


is  over-cored  so  that  the  sample  is  obtained  without  relative  movement  be¬ 
tween  the  two  blocks.  This  has  worked  satisfactorily  in  the  laboratory  on 
cores  of  sandstone  but  has  not  yet  been  tried  in  the  field. 

Coyne  and  Bellier,  the  Portugese  National  Civil  Engineering  Laboratory 
and  the  mining  group  at  Imperial  College  have  drilled  samples  parallel  to 
the  plane  of  joint  features  so  that  the  joint  specimen  occupies  a  diametral 
plane  of  a  core  sample.  Fig.  2.12  shows  such  specimens  after  a  test;  these 
NX  size  (2-1/8")  specimens  were  encased  in  Portland  cemer.t.  Coyne  and  Bel¬ 
lier  have  also  drilled  larger  specimens  (Fig.  2.13a).  Dr.  David  Pentz, 
at  Imperial  College,  England,  has  also  obtained  samples  by  this  technique. 

In  this  work,  he  used  a  nine  horsepower  diamond  drill  to  obtain  cores  nine 
inches  in  diameter.  The  cores  were  drilled  without  a  core  spring;  when  the 
full  depth  of  drilling  had  been  reached,  the  diamond  bit  was  removed  and  a 
core  spring  and  dummy  bit  were  installed.  Then  the  sample  was  pulled  by 
the  core  spring,  jacking  the  drill  against  the  collar,  to  produce  a  tension 
break  at  the  end.  The  core  containing  the  joint  specimen  could  then  be  re¬ 
moved  from  the  hole. 

After  removal,  the  joint  is  bound  so  that  it  will  not  be  disturbed  in 
transit  and  ptrafflnsA  (Fig.  2.13b)  or  otherwise  protected  so  that  the  natural 
moisture  content  will  be  preserved.  On  return  to  the  laboratory,  the  core 
is  fitted  in  the  shear  box  by  imbedment  in  either  concrete,  as  Coyne  and 
Bellier  have  done,  or  in  epoxy  resin,  as  Imperial  College  has  done.  Samples 
obtained  this  way  can  be  very  large,  up  to  1600  square  centimeters.  A  sur¬ 
prisingly  large  number  of  samples  have  been  obtained  by  Coyne  and  Bellier 
and 'Imperial  College  (Fig.  2,14). 

It  is  also  possible  to  drill  specimens  for  triaxial  joint  tests  by  ori¬ 
enting  the  drill  at  45-60°  to  a  naturally  occurring  joint.  A  triaxial  test 
with  a  naturally  occurring  joint  at  this  controlled  orientation  is  called  a 
multi-stage  triaxial.  Tests  of  this  type  were  performed  by  the  Missouri 
River  Division  of  the  Corps  of  Engineers  on  Piledriver  rock.  Such  tests 
were  also  performed  on  limestone  by  Heuze  and  Goodman^.  Professor  Jaeger 
in  Australia  has  conducted  many  tests  of  this  type. 

H-Heuze,  F.  and  Goodman,  R.  E.,  Mechanical  Properties  and  In  Situ  Beha¬ 
vior  of  the  Chino  Limestone,  Crestmore  Mine,  Proceedings  9th  Symposium  on 
Rock  Mechanics,  AIME,  1968. 
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Removal  of  Blocks  Containing  Joints.  It  is  possible  to  remove  large 
chunk  samples  of  rock  relatively  undisturbed  and  thus  obtain  large  numbers 
of  samples  in  the  laboratory.  One  method  of  obtaining  such  samples  is  by 
drilling  overlapping  holes  with  a  drilling  template  to  create  a  series  of 
slots.  The  U.  S.  Bureau  of  Reclamation  took  several  large  chunk  samples 
from  the  Grand  Coulee  Third  Power  House  by  this  technique.  The  firm  Coyne 
and  Bellier  has  introduced  a  wire  sawing  technique  to  speed  up  and  reduce 
the  cost  of  sampling  large  blocks.  Their  technique  consists  of  drilling  the 
corners  of  the  eventual  block  by  conventional  drill  holes;  the  wire  sawing 
apparatus  consists  of  'two  rods  with  pulleys  on  the  end  which  are  fitted  into 
the  drill  holes  (Fig.  2»15).  The  wire  runs  between  the  pulleys  which  are 
fixed  into  two  corner  holes  of  the  eventual  block  and  fed  out  of  the  drill 
holes  to  the  continuation  of  the  wire  loop.  In  principle,  four  slots  could 
be  cut  in  this  fashion  to  outline  the  four  walls  of  the  block;  but,  in  prac¬ 
tice,  several  of  the  sides  were  trepanned.  Then  the  bottom  of  the  block  is 
freed  by  sawing  with  an  ingeniously  contrived  loop  in  the  wire.  Eurenius  and 
Fagerstrom,  of  the  Swedish  firm  Vattenbyggnadsbyran,  have  taken  a  large  num¬ 
ber  of  block  samples  at  a  Syrian  dam  site  in  a  soft  rock  by  cutting  with  a 
power  saw.  Direct  shear  specimens  were  performed  in  the  laboratory  by  nor¬ 
mal  drilling  of  the  seams  contained  within  the  block  samples.  Other  methods 
of  cutting  blocks  which  may  prove  eventually  to  have  application  for  taking 
block  samples  are  electron  beam  cutting,  laser  cutting,  and  possibly  hydrogen 
arc  cutting. 

(3)  Methods  of  Testing  Joints 

Triaxial  Test  with  Oriented  Joint  —  Multi-stage  Triaxial.  As  stated 
above,  it  is  possible  to  conduct  a  triaxial  test  with  an  oriented  joint  so 
that  the  eventual  failure  is  by  the  mode  of  slipping  along  the  pre-selected 
surface.  The  specimen  is  first  failed  under  some  initial  confining  pressure, 
followed  by  step-wise  increase  of  confinement  and  sliding  by  increasing  axial 
pressure.  This  yields  a  value  for  the  joint  angle  of  friction.  This  tr?“ 
may  be  reasonable  in  stiff  joints  where  there  is  no  filling  or  compressible 
rock  zone  between  the  walls  of  the  joint.  A  finite  element  analysis  of  tri¬ 
axial  specimens  containing  inclined  joints  more  compressil  le  than  the  sur¬ 
rounding  rock  was  performed  as  a  collateral  experiment  in  the  course  of  the 
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previous  contracts.  The  results  indicated  very  uneven  stress  distributions 
along  the  joint  and  totally  different  concentrations  of  stress  in  the  upper 
and  lower  walls  of  the  joint.  Thus,  the  multi-stage  triaxial  is  probably 
not  a  good  test  for  seams  or  joints  of  shear  origin. 

Laboratory  Direct  Shear  Tests.  The  direct  shear  test  is  a  pertinent 
method  of  test  for  evaluating  the  relative  sliding  tendency  of  joint  blocks 
in  an  excavation  or  foundation.  However,  it  is  not  a  simple  matter  to  pro- 
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duce  a  reasonably  uniform  stress  distribution  inside  a  direct  shear  box.  If 
joints  were  perfectly  smooth,  it  would  be  possible  to  align  the  joint  inside 
a  direct  shear  chamber  and  displace  the  top  half  relative  to  the  bottom  with 
forces  almost  colinear.  Ideally,  these  forces  would  be  shears  along  the  top 
and  bottom  blocks  and  the  blocks  would  be  very  thin.  However,  since  the  top 
and  bottom  blocks  must  have  some  thickness  in  order  for  the  specimen  of  rock 
to  be  sufficiently  stiff,  it  is  not  practical  to  force  the  top  over  the  bottom 
by  shears  along  the  top  and  bottom  surfaces.  Therefore,  it  is  customary  to 
push  or  pull  the  top  relative  to  the  bottom  by  applied  forces  and  reactions 
transmitted  from  the  sides  of  the  block.  These  forces  introduce  overturning 
moments  into  the  shear  box  which  must  produce  contrary  moments  in  reaction 
from  the  top  and  bottom  of  the  shear  box,  causing  non-uniform  distribution 
of  normal  pressure  along  the  sliding  surface.  One  possible  result  is  indi¬ 
cated  in  Fig.  2.17*  If  the  joint  surface  is  rough,  furthermore,  it  is  neces¬ 
sary  to  provide  a  sufficient  gap  between  the  top  and  bottom  supporting  blocks 
to  allow  for  the  full  amplitude  of  the  irregularity  of  the  joint  surface. 

This  introduces  a  further  overturning  moment  and  a  tendency  for  non-uniform 
stress  distribution  in  the  shear  box.  There  are  further  moments  in  the  ver¬ 
tical  plane,  introduced  by  hinging  around  high  points  of  rock  along  the  sliding 
surfaces. 

Another  problem  of  direct  shear  tests  concerns  the  lateral  boundary 
conditions.  In  evaluating  the  strength  of  rock  joints  for  purposes  of  slope 
stability  analysis,  it  is  reasonable  to  conclude  that  the  prototype  blocks 
will  be  in  a  condition  of  plane  strain,  i.e.  that  no  side  displacements  will 
be  permitted.  In  order  to  achieve  this  boundary  condition  in  a  direct  shear 
machine,  only  one  direction  of  displacement  is  permitted.  This  is  the  Casa- 
grande  boundary  condition  designed  into  the  Casagrande  direct  shear  machines, 
distributed  in  soil  mechanics  laboratories  around  the  world.  With  irregular 


joint  surfaces  there  is  the  possibility  of  rotational  tendencies  in  a  hori¬ 
zontal  plane  because  of  rotary  moments  on  high  points  along  the  joint  surface. 
Thus,  the  upper  or  lower  block  tends  to  swing  into  the  sides.  In  the  yilane 
strain  type  of  shear  box,  this  is  impossible  because  of  the  lateral  confine¬ 
ment  offered  by  the  sides.  However,  the  resulting  distribution  of  side 
pressures  is  very  irregular  and  is  not  known;  therefore,  the  test  results 
may  be  quite  unreal.  Another  approach  to  this  problem  is  to  omit  the  sides 
of  the  shear  box  altogether.  The  problem  of  achieving  a  uniform  stress 
distribution  in  a  direct  shear  box  is  difficult.  A  through  study  by  modern 
analytical  methods  is  needed.  Shear  strength  characteristics  determined  in 
direct  shear  tests  depend  on  the  boundary  conditions  achieved.  Unfortunately, 
there  is  no  uniformity  in  the  construction  of  machines  in  use  around  the  world 
and  it  is  possible  that  the  results  of  one  worker  will  not  agree  with  the 
results  of  another  even  though  identical  specimens  are  tested. 

The  most  impressive  and  elaborate  direct  shear  machine  in  the  world 
was  developed  by  Dr.  David  Pentz  and  colleagues  at  the  Royal  School  of  Mines 
at  Imperial  College,  London,  England  (Fig.  2«  16) .  In  this  machine,  the 
horizontal  force  required  to  maintain  the  displacement  in  a  single  direction 
is  controllable  and  can  be  measured.  The  shear  box  takes  samples  up  to  nine 
inches  in  width.  Elaborate  controls  of  pressure  and  rate  of  displacement  have 
been  incorporated  in  its  design.  Unfortunately,  it  is  so  recent  a  development 
that  relatively  few  samples  have  been  tested  as  of  this  date. 

Casagrande  type  shear  boxes,  designed  for  soil  mechanics  use,  have  been 
adapted  for  joint  samples  in  many  laboratories  around  the  world.  In  addition, 
plane  strain  shear  boxes  of  much  larger  size  and  load  capacity  have  been 
constructed  for  work  on  rock  specimens.  The  most  important  of  these  machines 
is  the  one  used  by  Coyne  and  Bellier  in  Paris  and  developed  by  the  S.E.I.L.* 
(Fig.  2«lS) •  It  takes  samples  of  dimensions  0.3  meters  by  0.5  meters  or  less. 
It  is  desirable,  in  a  direct  shear  machine,  that  the  normal  load  follow  the 
position  of  the  center  of  the  test  surface  which  is  in  contact  at  any  time. 

In  the  Coyne  and  Bellier  machine,  this  is  achieved  by  providing  hinges  that 
allow  the  normal  pressure  ram  to  displace  with  the  relative  displacement  across 

*S.E.I.L.  *  Societe  d'Equipment  Industriel  et  de  Laboratoir 


21 

the  joint.  Movement  of  the  axis  of  normal  load  with  the  moving  block  can 
be  allowed  by  employing  dead  weights  or  hangars.  The  large  direct  shear 
machine  employed  at  Imperial  College  by  DeFreitas  and  Knill  (Fig.  2.19), 
and  the  fine  direct  shear  machine  employed  at  the  University  of  Illinois 
by  Deere  and  Coulson  (Fig.  2.20),  are  of  the  plane  strain  type  and  use  dead 
weight  loading. 

Professor  Krsmancv'.c  in  Yugoslavia  constructed  a  large  machine  with¬ 
out  sides  (Fig.  2.21a).  Specimens  for  this  machine  have  the  dimensions 
40  x  40  x  20  centimeters;  giving  a  shear  surface  of  about  1600  square  cen¬ 
timeters;  the  maximum  normal  stress  is  about  40  kilograms  per  square  centi¬ 
meter  and  the  maximum  shear  stress  is  80  kilograms  per  square  centimeter 
on  these  large  areas.  The  Krsmanovic  device  is  designed  in  such  a  way  that 
the  applied  shear  force  is  inclined  about  4°  relative  to  the  surface  of 
shearing.  Inclination  of  the  applied  force  is  one  way  of  overcoming  some  of 
the  difficulties  with  rotary  moments  described  above.  However,  a  large 
inclination  of  the  applied  force  causes  normal  stresses  to  vary  appreciably 
with  the  shear  stress  applied. This  makes  it  very  difficult  to  test  specimens 
under  low  normal  stress.  The  Krsmanovic  device  uses  jacks  on  both  sides 
and  ball  bearings  both  above  and  below  the  shear  blocks  so  that  the  center 
of  the  surface  which  is  in  contact  at  any  time  does  not  move  during  the 
test  (Fig.  2.21b).  Therefore,  it  is  possible  to  use  a  fixed  loading  frame 
to  apply  the  normal  pressures. 

In  Situ  Direct  Shear  Tests.  It  is  possible  to  perform  direct  shear 
tests  in  the  field.  These  tests  are  accomplished  by  selective  excavation  to 
form  a  block  on  a  natural  surface  of  weakness.  The  surface  does  not  need  to 
be  horizontal. 

There  is  no  standard  method  of  performing  this  test  .  The  block  is 
invariably  confined  in  a  reinforced  concrete  or  a  composite  steel  concrete 
frame.  The  shear  load  is  usually  slightly  inclined.  The  normal  load  is 
sometimes  provided  only  by  the  block's  own  weight. 

Blocks  on  inclined  surfaces  are  under  shear  stress  before  the  test 
begins.  To  extend  the  elastic  information  of  the  test,  the  block  may  be 
pushed  up  the  hill,  through  zero. 


^Standards  for  in  situ  shear  tests  are  being  considered  at  present  by 
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Few  tests  of  this  type  have  been  performed,  in  hard  rocks,  in  the 
United  States,  owing  to  the  expense  in  labor  required.  Many  such  block 
shear  tests  have  been  conducted  by  European  organizations,  but  relatively 
few  on  natural  joint  surfaces.  Several  examples  are  shown  in  Figs.  2.22, 

2.23  and  2.24. 

(4)  Results  of  Shear  Tests 

Results  of  a  large  pumber  of  laboratory  direct  shear  tests  and  several 
in  situ  block  shear  tests  have  been  summarized  in  Tables  2.2  and  2.3.  These 
data  were  gathered  mostly  from  unpublished  reports  generously  made  available 
by  the  orgaiizations  concerned. 

The  load-deformation  curves  for  these  tests  have  been  represented  by 
values  for  shear  stress  and  displacement  at  the  peak  and  residual  points. 

In  some  cases,  a  yield  point  was  reported  as  Well  in  the  early  part  of  the 
curves.  Though  the  possible  combinations  of  values  for  these  quantities 
are  numerous,  several  generalizations  can  be  made  about  the  characteristics 
of  the  results  for  different  types  of  weakness  surfaces.  The  load  defor¬ 
mation  curves  were  classified  into  four  types  as  drawn  in  Fig.  2.25  and  sum¬ 
marized  in  Table  2  4. 

Type  1  stress  deformation  curves  rise  very  steeply  to  peak  stress  at 
very  low  deformations  and  then  quickly  fall  to  a  residual  value  which  may 
be  one  third  or  less  of  the  peak.  Healed  joints,  and  incipient  fractures 
belong  to  this  class. 

Type  2  curves  also  develop  the  peak  stress  at  very  low  deformations  but 
do  not  fall  so  sharply  to  the  residual,  further,  the  peak  strength  is  only 
slightly  above  the  residual.  Artificial  joints  formed  by  diamond  saw  catting 
rock  specimens  belong  to  this  class.  Polished  surfaces,  type  2a,  show  less 
rate  of  decay  of  strength  after  the  peak  than  do  unpolished  saw  cuts  (2b). 

The  shear  stiffness  is  higher  than  natural  joints,  but  is  still  low  enough 
in  number  (ca.  10^  kg/cm^)  to  contribute  significant  deformation  beyond 
that  of  the  intact  rock  if  the  joints  are  repeated  at  intervals  of  several 
feet  in  a  hard  rock  such  as  quartz  monzonite. 

Type  3  curves  show  lower  stiffness  and  numerous  secondary  peaks.  Second 
peaks  may  be  as  large  as  the  primary.  These  curves  result  when  testing  clean, 
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rough  joints,  such  as  extension  joints  in  granite.  The  irregularities  in 
the  load-deformation  curve  result  from  overriding  of  successive  asperities. 
Undisturbed  samples  of  such  joints  may  show  an  initial  concavity  as  tangen¬ 
tial  aperture  is  closed  (type  3b). 

In  classes  1-3,  it  mattered  little  whether  the  surfaces  were  dry  or 
soaked  with  water.  In  type  4  curves,  on  the  other  hand,  the  behavior  is 
modified  grossly  by  change  in  the  water  content.  These  curves  are  typical 
of  sheared  zones,  clay  filled  joints,  and  shale  partings.  At  low  moisture 
contents  (4a),  the  curves  show  high  stiffness  and  great  differences  between 
peak  and  residual  displacements,  although  generally  not  as  severe  as  with  the 
healed  joints  (type  1).  When  wet,  the  stiffness  is  greatly  reduced.  Thin 
zones  (4b),  when  wet,  show  "strain  hardening"  behavior  with  ultimate  strength 
considerably  greater  than  the  point  of  maximum  curvature  (denoted  "peak") . 
Thick  seams  have  elastic-plastic  stress  deformation  curves  when  wet  (4c). 

All  of  the  in  situ  tests  (Table  2.3)  yielded  type  4  curves,  probably  because 
these  expensive  tests  are  reserved  for  the  worst  discontinuities. 

It  is  difficult  to  increase  the  water  content  of  sheared  zones  and  clay 
seams  to  desired  values.  The  field  moisture  content  can  be  retained  by 
careful  sealing  in  the  field.  To  increase  the  moisture  content  beyond  this 
value,  soaking  is  necessary.  Soaking  the  specimens  without  surcharge  may 
lead  to  excessive  moisture  contents  —  well  beyond  those  expectable  in  the 
field.  Since  the  permeability  becomes  low  under  normal  pressure,  soaking 
under  surcharge  tends  to  wet  the  edges  more  than  the  center.  Coyne  and 
Bellier  have  tried  introducing  water  into  drill  holes  parallel  and  also  per¬ 
pendicular  to  the  seam,  with  moderate  success. 

Typical  values  for  the  peak  displacements,  tangential  stiffness,  and 
ratio  of  peak  to  residual  stresses  are  given  for  all  types  of  surfaces  in 
Table  2.4.  It  is  emphasized  that  these  are  only  typical  values;  consider¬ 
able  deviation  exists  in  all  these  quantities.  Both  the  peak  displacement 
and  shear  stiffness  tend  to  increase  with  increasing  normal  stress.  When 
variation  in  normal  pressure  is  ignored,  the  deviation  of  peak  displacements 
is,  however,  far  less  than  the  deviation  of  stiffness.  In  fact,  if  no  test 
were  possible^ as  a  rough  estimate  of  initial  shear  stiffness,  one  could  cal¬ 
culate  the  strength  from  the  friction  angle  (which  is  usually  between  30° 
and  40°)  and  divide  it  by  the  peak  displacement  value  selected  from  Table  2.3. 
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One  of  the  most  interesting  conclusions  from  the  test  results  reported 
in  Table  2.2  is  the  large  displacement  required  to  reach  residual  strength 
values  —  frequently  several  centimeters. 

Data  on  displacements  normal  to  the  shear  plane  were  not  presented 
in  Tables  2.2  and  2.3.  In  all  laboratory  tests  where  normal  displacements 
were  given,  shearing  displacement  was  accompanied  by  dilations  (thickening 
of  the  joint  zone).  The  peak  dilation  occurs  well  after  the  peak  stress 
at  low  normal  pressures,  but  occurs  right  at  the  peak  shear  stress  in  tests 
run  under  high  normal  pressure.  Several  of  the  in  situ  block  shear  tests 
showed  slight  closing  of  the  joint  in  the  early  part  of  the  shear-displace¬ 
ment  curve;  however,  dilation  always  occurred  eventually. 


2.5  APPLICATION  TO  PILEDRIVER  -  DL  DRIFT 

The  methods  described  in  this  section  were  used  to  model  the  cross 
section  of  DL  Drift  at  station  0  +  70.  This  section  is  at  the  intersection 
of  DL  1  and  DL  2.  The  drift  was  16  feet  in  diameter,  containing  pre-ten- 
sioned  rock  bolts  on  a  regular  pattern,  with  one  bolt  per  3.25  square  feet 
of  the  wall.  There  are  two  layers  of  chain  link  fabric  held  by  the  bolts. 

Originally,  it  was  planned  to  make  a  number  of  models  at  many  different 
cross  sections.  However,  the  model  became  so  complex  that  it  demanded  a 
great  deal  of  time.  An  effort  was  made  to  reduce  the  complexity  of  each 
geologic  cross  section;  however,  results  of  study  with  the  DL  drift  section 
showed  that  more,  rather  than  less,  complexity  would  be  required  to  achieve 
meaningful  results.  Further  refinement  would,  unfortunately,  call  for  know¬ 
ledge  of  the  detailed  geologic  and  material  properties  exceeding  that  which 
was  known  before  the  test.  The  purpose  of  the  analysis  presented  here  is 
to  judge  our  capacity  to  perform  a  meaningful  quantitative  analysis  of  an 
actual  geological  cross  section. 

The  finite  element  mesh  prepared  for  DL  drift  is  presented,  without 
nodal  point  and  element  numbers,  in  Fig.  2«26a.  This  is  the  most  complex 
finite  element  mesh  known  to  the  author.  Its  900  nodal  points,  667  total 
elements,  and  398  joint  elements  completely  filled  the  core  storage  on  the 
IBM  7094  computer  available  at  the  time  it  was  prepared.  However,  the  new 
CDC  6400  computer  with  extended  core  storage  recently  installed  at  the 
University  of  California  allows  a  five-fold  increase  in  the  problem  size 
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and  a  significant  enlargement  of  the  allowable  bandwidth.  Debugging  the 
mesh  was  made  possible  by  using  a  plotter  to  draw  the  element  boundaries 
as  punched  on  the  input  cards.  (Fig.  2*26  is,  in  fact,  the.  machine  plotted  . 
mesh . ) 

Results  of  two  cases,  the  parameters  of  which  are  summarized  in  Table 
2.5,  will  be  described.  In  case  1,  the  joints  were  assigned  very  stiff, 
strong  properties,  as  if  they  are  only  incipiently  developed,  or  are  healed. 
In  case  2,  the  joint  strength  was  reduced  and  the  deformability  assigned  to 
joint  and  rock  elements  after  failure  was  very  great.  In  a  third  case,  the 
different  sets  of  weakness  surfaces  were  treated  differentially;  the  shears 
(A  joints)  were  assigned  the  weakest  and  most  deformable  properties  and  the 
mineralized  B  joints  were  assigned  the  stiff,  strong  properties  of  incipient 
joints.  The  results  from  study  of  case  3  did  not  contrast  strongly  with 
case  2,  and  therefore  are  not  presented  here. 

The  failure  criterion  adopted  for  the  rock  in  case  2  was  more  severe 
than  for  prototype  Piledriver  rock  as  N  was  taken  equal  to  zero.  The  actual 
rock  being  slightly  stronger  than  this,  the  extent  of  rock  breakage  indicated 
from  computations  is  somewhat  too  great. 

The  results  for  cases  1  and  2  are  presented  in  Figs.  2*27  and  2.28. 

The  displacements  are  plotted,  to  scale,  around  the  walls  of  the  tunnel  and 
at  selected  points  behind  the  wall.  The  blast  side  is  to  the  left  in  these 
figures.  Elements  that  have  failed  are  indicated  by  "R"  if  they  are  rock 
elements  and  by  "j"  if  they  are  solid  elements.  Double  lines  indicate 
opening  of  joints. 

(1)  Case  1  —  Results 

In  Fig.  2.27a,  displacements  are  plotted  to  the  scale  of  the  drawing 
after  the  second  5000  psi  load  increment  and  the  failed  elements  are  shown 
after  the  total  pressure  reached  5000  and  10,000  psi.  In  Fig.  2.27b,  dis¬ 
placements  and  failed  rock  elements  are  shown  after  15,000  psi  (third  incre¬ 
ment),  the  failed  joints  were  not  shown  owing  to  a  mistake  which  disturbed 
the  printing  of  joint  stresses  after  the  third  increment. 

The  results  plotted  in  Fig.  2.27  show  that  the  elements  tended  to 
fail  along  the  steep  A  joints  above  and  below  the  tunnel.  This  joint  sys¬ 
tem  is  roughly  perpendicular  to  the  blast  front.  A  large  failed  zone 
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was  localized  in  the  tunnel  roof  above  the  lee  side.  Rock  blocks  in  this 
zone  were  in  tension  resulting  from  eccentric  loading  as  the  blocks  shifted 
along  joints.  Little  additional  damage  was  done  after  the  first  Increment 
of  the  blast  (9  additional  failed  elements  in  the  second  increment  as  com¬ 
pared  to  32  in  the  first).  The  opening  of  joints  and  the  relative  displace¬ 
ment  of  adjacent  nodal  points  is  shown  after  15,000  psi  on  Fig.  2.27b.  The 
whole  blast  side  moved  into  the  opening  almost  uniformly. 

(2)  Case  2  —  Results 

In  Fig.  2.28,  a  much  more  drastic  but  generally  similar  behavior  pat¬ 
tern  is  shown  for  case  2.  Here  the  properties  of  the  rock  and  joints  were 
weaker  both  before  and  after  failure.  The  displacements  after  10,000  psi 
were  so  large  that  the  mesh  could  not  be  used  for  additional  increments. 
Failures  were  plotted  after  the  first  increment  of  pressure  (2500  psi)  and 
the  displacements  were  plotted  after  the  second  increment  with  cumulative 
pressure  equal  to  5000  psi.  On  the  blast  side  of  the  tunnel,  the  wall  moved 
in  1.5  feet,  while  on  the  lee  wall,  the  movements  were  negligible.  Very 
small  displacements  resulted  from  the  first  2500  psi  increment  but  there 
were  numerous  failed  elements,  the  great  deformability  assigned  to  failed 
elements  for  the  second  increment  led  to  the  very  large  displacements  plotted. 
Joint  openings  occurred  in  the  lower  part  of  the  blast  side  while  blocks  slid 
along  flat  joints  in  the  roof. 

(3)  Comparison  with  Actual  Damage 

The  deformed  shape  of  the  tunnel  and  location  of  broken  rock  zones 
indicated  by  these  analyses  can  be  compared  with  the  actually  observed  failure 
pattern  shown  in  Fig.  2.29.  The  analysis  predicts  large  inward  movements 
of  the  whole  wall,  with  little, rock  breakage,  on  the  blast  side.  The  ana¬ 
lysis  also  indicates  extensive  localized  crushing  of  rock  above  the  roof 
on  the  lee  side.  This,  in  general,  is  what  occurred. 

(4)  Conclusion 

The  computing  mesh  was  far  too  small,  and  the  geological  information 
was  far  too  imprecise  to  allow  more  refined  studies  of  these  sections.  The 
indication  from  the  study  is,  however,  that  given  sufficiently  accurate 


geological  and  materials  property  data,  analysis  can  now  predict  behavior 
of  tunnels  in  jointed  rock.  Further  comments  are  made  on  the  ramifications 
of  this  conclusion  in  Section  4  of  the  report. 
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Table  2.1 


SUMMARY  OF  ROCK  PROPERTIES 
(Compression  Positive) 


DEFORMABILITY 

11.0  x  10^  psi 
8.0  x  10^  psi 
0.3 

STRENGTH 

Triaxial  Data 
4>  =  56° 
c  =  3600  psi 

Tensile  Strength 

-Gt  =  -1450  psi 

Adopted  Equation  of  Failure  Surface  (psi  units*) 

1*  ^oct  >  0 

fails  when  Toct  =  1300  +  2.95(aoct)0'91 

N  =  1300  psi 

D  =  2.95 
B  =  0.91 

2.  aoct  <  0 

fails  when  0m^n  <  -1450  psi 


Modulus  of  elasticity  in  compression 
Modulus  of  elasticity  in  compression 
Poisson's  ratio 


*To  convert  to  psf  units,  multiply  N  by  144  and  solve: 
Dpsf  *  lo8_1  Uog  Dpsi  +  (1  -  B)  log  144] 

In  psf  units: 

N  =  1.87  x  105 

D  =  5.04 

B  =  0.91 


TABLE  2.2  RESULTS  OF  LABORATORY  DIRECT  SHEAR  TESTS  ON  ROCK  SPECIMENS  WITH  JOINTS  SHEET  2  OF  3 
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Filled  joints,  sheared  zones,  shale  partings, 
and  smooth  bedding 

4a  dry  or  slightly  moist  0.12 
4b  wet;  thin  seam  0.26 
4c  wet:  thick  seam  0.30 


Table  2.5 


ROCK  AND  JOINT  PROPERTIES  USED 
IN  STUDY  OF  SECTION  AT  0  +  70  IN  DL  DRIBT 


Case  1  Case  2 


Before 

After 

Before 

After 

Failure 

Failure 

Failure 

Failure 

Rock  Deformability 

E  in  compression  (psi) 

11.0  x  106 

11.0  x  105 

11.0  x  106 

0.01 

E  in  tension  (psi) 

8.0  x  106 

8.0  x  103 

8.0  x  106 

0.01 

v  (Poisson's  ratio) 

0.3 

0.4 

0.3 

0.4 

Rock  Strength  Parameters* ** 

N  psi 

1300 

0 

D  (correspondiag  to  psi) 

2.95 

6.60 

B 

0.91 

0.88 

Tensile  Strength  (psi) 

1450 

Joint  Deformability 

** 

** 

Normal  stiffness  psi/in 

2.9  x  104 

2.9  x  104 

2.9  x  104 

2.9  x  104 

Shear  stiffness  psi/in 

2.9  x  105 

2.9  x  104 

2.9  x  104 

0.6  x  10"3 

Maximum  closure  (inches) 

0.12 

0.12 

Joint  Strength 

Cohesion  psi 

2,7B0 

100 

Friction 

56° 

31° 

*See  Section  2.4 

**Stif fnesses  both  set  to  zero  if  fails  in  tension  (opens) 


FIGURE  Z.l 

PRE-SHOT  GEOLOGIC  LOG 
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EXPLANATION  OF  FIGURE  2.2 

Method  of  Projecting  Geological  Data  frcm  Geological  Log 
of  Tunnel  to  Draw  Cross  Sections 

1*  The  strike  in  relation  to  the  tunnel  can  he  found  by  three  dif¬ 
ferent  methods: 

a.  Revolve  spring  line  in  geological  log  to  tunnel  diameter  and 
measure  strike  a’V; 

b.  Measure  tangent  to  curve  at  the  crown; 

c.  Use  measured  strike  in  the  field. 

2.  Project  the  intersection  of  the  strike  at  the  crown(s)  to  sec¬ 
tions  desired.  For  example,  Joint  is  at  crown  level  at  t  and  r 
respectively  in  sections  A  and  B. 

3.  Project  the  strike  of  the  joint  at  spring  line  level  to  the 
sections  desired;  in  the  example,  Joint  is  at  spring  line  ele¬ 
vation  at  points  n  and  .£  respectively  in  sections  A  and  B. 

k.  Line  tn  is  the  trace  of  the  plane  on  section  A.  line  ri  is 
the  trace  of  the  plane  on  section  B. 

Note:  Points  x  and  y  show  where  the  joint  plane  intersects  the 
tunnel  boundary  in  section  B.  It  does  not  intersect  the  boundary 
in  section  A. 

5.  Find  dip  (6  )  from  distance  sm.  6  =  tan”1  -unpe-^r^^ —  ,  or 
use  dip  measured  in  the  field. 


METHOD  OF  PROJECTING  GEOLOGICAL  DATA 


(PSIxIOOO) 


failure  curve 
corresponding 
to  B  =  1 1 


^  failure  curve  corresponding  to 
tension  and  compression  tests 


failure  curve  corresponding  to  N  =  0 


calculated  value  of  N  corresponding  to  B  =  1 
_ _  tension  test 


4  6  8  10  12  14  16  18  20 

Cr0CJ  (PSIxIOOO) 


FIGURE  2.5  OCTAHEDRAL  STRESSES  AT  FAILURE 
FOR  PILEDRIVER  ROCK 


FIGURE  2.6  LOG -LOG  PLOT  OF  OCTAHEDRAL  STRESSES 
AT  FAILURE 


Figure  2.9  Artificial  joint  specimen  produced  by  longitudinally  diamond  saw 
cutting  an  NX  core  sample.  (Courtesy  Univ.  Illinois,  Dept.  C.E. ) 
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Figure  2.10  Preparation  of  direct  shear  specimen  containing  am  artificial  extension 
joint,  a)  Mold  placed  around  fitted  blocks.  Plane  of  joint  will  be 
protected  with  putty  while  mold  is  filled  with  epoxy,  b)  Direct  shear 
specimens  after  test.  (Courtesy  of  Imperial  College  of  Science  and 
Technology,  Dept,  of  Geology) 


Figure  2.1k  Large  direct  cheer  specimens  obtained  by  drilling  nine  inch  cores 
parallel  to  selected  Joint  surfaces.  (Courtesy  Imperial  College, 
Royal  School  of  Mines) 
a)  specimens  as  received  from  field 
\  b)  after  removing  protective  Jacket 

j  c)  Mounted  in  shear  box,  ready  for  testing  (confining  bands  x 

>  be  cut  before  test 

d)  specimens  after  test 

I  e)  croea  section  through  a  speciaen  shoring  rock  r> 

I  filler  in  epoxy 
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figure  2.15 


Wire  sawing  technique  used  "by  Coyne  and  Bellier  to  extract  block 
s Maples  containing  Joints.  The  vires  run  around  pulleys  (not 
visible)  at  the  ends  of  the  rods.  (Courtesy  of  Coyne  and  Bellier) 
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figure  2.16 


Large  direct  shear  machine  at  Imperial  College,  Royal  School  of 
Mines.  (Dr.  Pentz  to  left.)  (Courtesy  Imperial  College) 


ACTUAL  SITUATION 


DESIRED  STRESS 
DISTRIBUTION 


POSSIBLE  STRESS 
DISTRIBUTION  IF 
ROTATION  IS  IMPOSSIBLE 


FIGURE  2.17  STRESS  DISTRIBUTION  IN  A  DIRECT  SHEAR  BOX 
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Figure  2.18  Casagrande  type  (plane 

strain)  direct  shear  machine  used 
hy  Coyne  and  Belli er.  Takes  sam¬ 
ples  up  to  30  hy  50  centimeters. 
(Courtesy  of  Coyne  and  Bellier) 


Figure  2.19  Plane  strain  direct  sheer 
machine  employed  at  Dept,  of  Geology, 
Imperial  College.  Uses  dead  weight 
loading  on  hangers.  (Courtesy  of 
Dr.  De Freitas) 
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Figure  2.20  Plane  strain  direct  shear  machine  at  University  of  Illinois,  Dept, 
of  Civil  Engineering.  (Courtesy  of  Prof.  Deere) 
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CROSS  SECTION  30 


Figure  2.21  Krsmanovic  shear  machine,  Sarajevo,  Yugoslavia,  a)  Prof.  Krsmanovic 
(right)  and  colleagues  at  shear  machine;  b)  schematic  diagram  (from 
Oeotechnlque,  June  1967,  p.  ll>6) 
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Figure  2.22  In  situ  shear  test  conducted  at  Vouglnns  dam.  (Courtesy  of  Coyne 
and  Bellier) 

a)  excavation  of  a  pillar  containing  the  Joint  to  he  tested 

b)  test  block  before  encasing  in  concrete  frame 

c)  the  whole  test  apparatus 

d)  Jacks  used  to  apply  tangential  force.  Vibration  wire  gauges 
used  to  measure  displacements 
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Shear  test  an  a  plane  of  mechanical  discontinuity  tcith  shear 

movement  within  the  rock  mass. 

1)  Concrete  block  with  measured  points  A  and  B; 

2)  Roller  bearing; 

3)  Strain  gauge  dynamometers; 

d)  Hydraulic  presses  for  pressures. of  SO  or  100  Mp; 

5)  Blaster  levelling  wedge:  a)  Broken-off  rock  round  the 
block;  b)  Cardboard  lubricated  tcith  vaseline;  c)  points 
at  which  the  movement  on  a  plane  of  mechanical  dis¬ 
continuity  is  measured;  l  and  II)  planes  of  mechanical 
discontinuity. 


Figure  2.23  In  situ  shear  tests  conducted  in  Czechoslovakia,  a)  Conventional 
test,  on  a  horizontal  plane;  h)  test  on  a  steeply  inclined  plane. 
(Courtesy  Dr.  Karel  Drozd.  Reprinted  from  Proc.  of  Geotechnical 
Conf.,  Oslo  196 7»  P»  266.) 


Figure  2.24  In  situ  shear  tests  on  a  clay  parting  along  bedding  in  shale  on  a 
natural  slope  at  Khajuri  Kach  dan,  Pakistan.  (Courtesy  Dr. 
Kujundzic,  Jar os lav  Cerni  Institute,  Belgrade,  Yugoslavia) 
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DISPLACEMENTS  AFTER  10,000  PSI, 
FAILED  ELEMENTS  AFTER 
5000  AND  10,000  PSI. 

R  *  ROCK  FAILURE 
J  *  JOINT  FAILURE 
R  AND  J  SIGNIFY  FAILURE 
ONLY  AFTER  10,000  PSI 

0  2  4  FEET 

i _ I _ i 

SCALE (OF 

DRAWING  AND  DISPLACEMENT) 

- ORIGINAL  POSITION 

-  DEFORMED  POSITION 

- TUNNEL  WALL 


FIGURE  2.27o  DL  DRIFT  0+70 
DISPLACEMENTS  AND  FAILED 
ELEMENTS  FOR  CASE  I 


DISPLACEMENTS  AND 
FAILED  SOLID  ELEMENTS  (R) 
AFTER  15,000  PSI. 

FAILED  JOINTS  NOT  SHOWN 
DOUBLE  LINE  *  OPEN  JOINT 

024  FEET 

i _ i _ i 

SCALE  (OF  DRAWING  AND 
DISPLACEMENT) 


- ORIGINAL  POSITION 

- DEFORMED  POSITION 

- TUNNEL  WALL 


FIGURE  2.27  b  DL  DRIFT  0+70 
DISPLACEMENTS  AND  FAILED 
ELEMENTS  FOR  CASE  I 


DISPLACEMENTS  AFTER  5000  PSI. 
FAILED  ELEMENTS  AFTER  2500  PSI 

024  FEET 

i _ i _ i 

SCALE  (OF  DRAWING  AND 
DISPLACEMENT) 

R  ROCK  FAILURE 
J  JOINT  FAILURE 

(BY  OPENING  OR  SHEARING) 

- ORIGINAL  POSITIONS 

-  DEFORMED  POSITION 


FIGURE  2.28  DL  DRIFT  0+70 
DISPLACEMENTS  AND  FAILED 
ELEMENTS  FOR  CASE  2 


BLAST  SIDE  TO  LEFT 

8  =  BROKEN  ROCK  IN  MESH 

M  =  BROKEN  ROCK  COVERING  INVERT 

- ORIGINAL  POSITION 

- DEFORMED  POSITION 

FIGURE  2.29  POST  SHOCK  CROSS  SECTION 

STATION  0+70  -  DL  DRIFT 
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3.  ANALYSIS  OF  CR  AND  DL  DRIFTS  BY  JOINT  INFLUENCE  DIAGRAM  METHOD 


3.1  DESCRIPTION  OF  PROCEDURES 

The  inclusion  of  joint  elements  with  prescribed  stiffness  and  strength 
values  in  the  powerful  finite  element  method  allows  testing  of  sophisticated 
mathematical  models.  The  stress  distributions  employed  in  evaluating  the 
strength  of  a  tunnel  section  by  this  method  take  into  account  the  effect  of 
discontinuities,  of  heterogeneities,  and  of  local  failures.  Tunnel  liners 
and  rock  bolts  can  be  incorporated.  However,  the  discontinuities  are  en¬ 
tirely  two-dimensional  in  their  representation.  Further,  the  method  is  rela¬ 
tively  tedious  if  large  numbers  of  joints  are  involved. 

Another  approach,  termed  the  ubiquitous  joint  method  in  previous  reports, 
permits  an  examination  of  the  influence  of  varying  distributions  of  joints 
at  a  large  number  of  sites  with  relatively  little  effort.  The  method  consists 
of  computing  regions  of  influence  for  prescribed  joint  orientations  in  a 
given  stress  field  and  comparing  the  influence  diagrams  with  site  geology. 

The  steps  in  the  analysis  are  as  follows: 

1.  Construct  cross  sections  of  the  drifts  showing  the  location  and 
orientation  of  joints  and  faults.  Divide  the  joints  into  sets 
on  the  basis  of  preferred  orientation. 

2.  Calculate  the  region  of  influence  of  each  joint  set  assuming 
some  state  of  stress  acting  as  a  result  of  the  blast.  A  simple 
Kirsch  solution^  stress  state  was  assumed,  with  the  blast  re¬ 
placed  by  a)  a  static  load  equal  to  P  and  b)  confinement  in  the 
perpendicular  directions  equal  to  (v/1  -v) (P) ,  where  P  is  the 
peak  pressure  of  the  blast  and  v  is  Poisson's  ratio. 

3.  Overlay  the  geologic  sections  on  the  influence  diagrams  and  note 
the  locations  where  joints  of  each  set  can  slip. 

4. '  Using  the  results  of  (3),  and  considering  the  whole  network  of 

joints  at  each  section,  sketch  the  probable  zones  of  rock  fall- 


12obert,  L.  and  DuVall,  W. ,  Rock  Mechanics  and  the  Design  of  Structures 
in  Rock,  Wiley,  1967,  p.  101. 
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out  and  failure.  In  rock  bolted  sections,  the  capacity  of  the 

13 

rock  bolts  to  prevent  rock  movement  is  considered  at  this  stage  * 

Fig.  3.1a*  above  the  sequence  of  test  sections  and  stationing  referred 

to  in  this  text.  Fig.  3* lb*  shows  the  main  joints  and  faults,  and  labels 

those  sections  where  major  failures  occurred.  All  joints  and  sheared  tones 

occurring  in  D  drift,  CR  Fast  West,  and  CR  Worth  were  listed  and  assigned 

to  one  of  nine  joint  sets  (A-X)  on  the  basis  of  orientation.  Geologic  cross 

sections  were  constructed,  as  described  in  Section  2.1,  every  10  feet  in  CR 

drift  and  in  DL  1  and  2.  At  each  section,  the  mean  orientations  of  each  of 

the  joint  sets  were  calculated  and  joint  influence  areas  were  determined, 

assuming  4  equals  31°,  and  c  equals  either  100  or  1000  pal.  She  stress 

field  used  to  calculate  the  joint  Influence  cones  correspond  to  the  Kirsch 
ik 

solution  for  stresses  around  a  circular  tunnel,  in  a  biaxial  stress  field, 
superimposed  on  the  component  of  the  blast  pressure  that  is  longitudinal  to 
the  tunnel,  as  listed  in  Table  3-1* 

'Ite  computer  programs  used  in  this  study  are  presented  in  Appendix  2. 
The  computer  was  used  to  generate  the  stress  fields  around  the  openix^s,  by 
the  Kirsch  solution,  and  to  calculate  and  plot  the  joint  influence  regions 
for  given  joint  parameters. 

In  this  analysis,  it  has  been  assumed  that  the  joints  do  not  modify 
the  stress  distribution,  l.e.  that  they  are  stiff.  Further,  it  is  assumed 
that  rock  bolted  sections  and  lined  sections  can  be  studied  without  taking 
into  account  their  strengthening  effect  until  the  last  step;  in  other  words, 
that  they  similarly  do  not  modify  the  stress  distribution.  Both  of  these 
assumptions  are  Indefensible  on  detailed  examination  but  they  simplify  the 
problem  to  such  a  degree  that  it  becomes  practical  to  apply  this  procedure. 
The  question  to  be  explored  here  is  to  what  extent  this  simplified,  step 
by  step  joint  influence  analysis  yields  reasonable  tunnel  behavior  estimates. 
Ihls  will  be  discussed  section  by  section. 


♦Reproduced  from  FOR  4015,  Chapter  4. 

13U.  S.  Army  Engineer  District  Cmaha,  Technical  Reports  2  and  3 
^op.  cit.,  p.  58. 
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3.2  CR  EAST-WEST  DFUT  —  RESUETS  OP  ANALfSIS 

Analyses  were  made  every  10  feet  in  CR  drift  beginning  at  station 
0  +  20  (CR  2A).  The  excavated  diaaeters  and  lining  types  of  the  sections 
in  this  drift  are  listed  in  Table  3.2. 

Results  of  analyses  are  reported  for  the  rock  bolted  or  unlined  sec¬ 
tion  from  0  +  60  to  1  +  70.  The  strike  and  dip  of  joints  and  sheared  zones 
deterai^ed  from  Figs.  3-lc,d,  pre-construction  logs,  are  listed  in  Table  3.3. 
In  Table  3*^  the  strike  relative  to  the  tunnel  axis  (0),  the  dip  of  the 
joint  down  from  horizontal  (D),  and  the  dip  (  6)  relative  to  0°  along  the 
radius  from  the  weapon  point  are  listed  for  each  joint  set  at  each  station, 
From  these  data,  the  angle  (0)  between  the  normal  to  the  joint  plane 
and  the  axis  of  the  tunnel  has  been  determined  as  described  in  Technical 
Report  No.  3^*  Jbint  influence  diagrams  have  been  constructed  for  all 
the  sets  of  joints  for  the  stress  state  corresponding  to  a  blast  pressure 
of  10,000  psi.  Joint  sets  C  and  6  were  found  to  have  no  influence  region 
around  the  tunnel  at  a  number  of  stations,  as  described  by  the  note  "no 
Joint  failure"  in  Table  3«^«  Two  cases  were  analyzed,  corresponding  to 
assumptions  of  Joint  cohesion  of  100  and  1000  psi.  Influence  diagrams 
for  these  cases  are  shown  in  Figs.  3.2a,b  and  3*3a,b  respectively. 

Four  types  of  figures  (a,  b,  c,  d)  have  been  prepared  for  the  step 
by  step  analyels.  At  each  station  considered,  first  a  geological  cross 
section  is  shown  (a).  (These  are  to  a  standard  size  but  varying  scale). 

The  second  figure  in  each  set  (b)  shows  only  the  joints  and  shears  which 
occur  within  their  zone  of  influence.  It  is  derived  by  superposition  of  the 
Influence  diagrams,  for  each  joint  set  in  turn,  on  the  geologic  cross  section 
(a).  The  third  and  fourth  figures  in  each  set  (c  and  d)  show  the  expectable 
poet-shock  profile  for  the  tunnel;  these  figures  are  obtained  by  superim¬ 
posing  figure  b  on  the  geologic  cross  section  to  determine  which  blocks 
have  kinematic  freedom  to  move  into  the  opening.  Figures  c  and  d  corres¬ 
pond  to  cohesion  of  100  and  1000  pei  respectively.  No  expectable  post¬ 
shock  profiles  were  prepared  for  concrete  or  composite  lined  tunnels.  For 
rock  bolted  sections,  it  was  assumed  that  the  bolts  and  wire  mesh  could 
restrain  only  up  to  several  feet  of  broken  rock  from  being  accelerated  into 
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the  tunnel.  Figs.  3*4  through  3*14  present  the  results  of  the  analysis 
for  the  CR  (3f)  drift. 

(i)  station  0  +  60  CR3  -  Fig.  3«4(a-d) 

The  analysis  weights  the  A  joints  heavily  in  the  roof  and  on  the  blast 
side  Just  below  spring  line.  D  and  E  joints  at  240°  and  270°  are  also  sig¬ 
nificant.  The  closely  fractured  zone  in  the  roof  is  formed  of  uniafluential 
C  joints  but  their  intersections  with  the  A  joints  define  blocks  which  would 
tend  to  fall  into  the  tunDel.  Thus,  the  analysis  predicts  rock  movement 
inward  at  90-100°,  21*0-270°  and  340-350°. 

Actual  damage  can  be  seen  on  the  post-shock  cross  section  at  0  +  59 
(Fig*  3*l£)*  The  closely  spaced  fractures  in  the  roof  did  not,  in  fact-, 
cause  roof  falls  in  this  section  except  where  they  intersected  the  north 
42  west  shear  and  this  led  to  a  major  rock  movement  at  about  110°.  Also, 

k 

towards  the  invert  at  the  base  of  the  lee  wall,  there  was  a  minor  rock 
movement  inward.  Rxere  was  very  substantial  inward  movement  on  the  blast 
wall  from  slightly  above  0°  down  into  the  invert.  In  total,  the  damage 
to  the  section  is  not  substantial.  Tae  comparison  of  the  actual  damage 
pattern  and  its  relationship  to  the  geology  with  the  damage  pattern 
predicted  by  the  joint  influence  diagram  analysis  shows  conformable  areas 
of  agreement.  It  is  interesting  that  rock  holt  deformeter  data  at 
station  0  +  59  were  lost  at  100°  due  to  the  extensive  rock  movement 
in  this  region.  These  deformeter  data  also  documented  large  rock  movement 
on  the  blast  wall  side  —  4.7  inches  increase  in  length  over  the  16  foot 
long  deformeter  length. 


(2)  Station  0  +  70  CR3  -  Fig*  3«5(a-d) 

The  geologic  section  at  0  +  70  shews  flat  lying  joints  closely  spaced 
above  the  roof,  a  shear  zone,  called  1-2,  in  the  blast  wall  with  apparent 
dip  towards  the  opening,  and  a  series  of  steep  3  Joints  behind  the  lee  wall. 
The  cross-section  is  very  complicated.  The  joint  influence  analysis  predicts 
very  strong  influence  of  2  A  Joints  which  intersect  the  roof  and  are  inclined 
with  apparent  dip  towards  the  weapon  point.  The  results  of  the  analysis  are 
shown  in  Fig.  3.5  c  and  d.  A  major  rock  fall  from  the  roof  was  predicted 
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together  with  a  rock  fall  from  the  lee  wall  down  to  about  spring  line  de¬ 
limited  by  the  steeply  dipping  B  Joints.  Inward  movement  is  also  predicted 
from  the  blast  wall  due  to  the  interaction  of  D  and  F  joints  and  the  1-2 
shear  zone.  The  actual  damage  pattern  at  this  location  can  be  seen  on  Fig. 
3.16  for  station  0  +  71.  Major  rock  falls  occurred  from  the  roof.  Some 
rock  was  broken  but  was  restrained  by  the  mesh  at  about  80°.  From  90-120°, 
major  rock  movement  occurred  and  this  movement  continued  with  decreasing 
thickness  down  into  the  lee  wall  to  about  spring  line,  apparently  along  the 
steeply  inclined  joints.  Also  there  was  inward  movement  on  the  blast  wall. 
The  rock  bolt  deformeter  data  at  0  +  71  were  lost  in  the  roof  because  of  the 
rock  falls.  Also,  very  large  inward  movements  were  observed  on  the  blast 
wall  side.  The  agreement  with  the  joint  influence  analysis  appears  to 
be  good. 

(3)  Station  0  +  80  CR4  -  Fig.  3«6(a-d) 

The  geological  features  of  the  cross  section  at  0  +  80  are  a  series  of 
steeply  inclined  B  joints  far  into  the  lee  wall,  a  closely  fractured  zone 
of  G  joints  below  the  invert,  and  a  series  of  shears  and  joints  of  the  A 
orientation  traversing  the  tunnel  and  inclining  toward  the  weapon  point. 
There  is  also  a  regular  spacing  of  C  joints  traversing  the  tunnel  and  in¬ 
clined  away  from  the  weapon  point.  Absent  from  the  section  are  B  joints 
traversing  the  tunnel  and  A  joints  above  the  roof  of  the  tunnel.  As  will  be 
seen  in  the  discussion  of  CR2  and  the  unlined  unsupported  drill  slot,  the 
intersection  of  A  and  B  joints  in  the  roof  provides  the  situation  for  large 
rock  movements  into  the  tunnel.  The  joint  influence  analysis  predicts  rock 
movement  at  100-120°,  mainly  on  A  joints.  The  actual  damage  to  be  seen  in 
Fig.  3.16  and  the  cross  section  for  station  0+79  shows  large  rock  falls 
restrained  by  the  wire  mesh  from  about  20-90°,  These  falls  were  defined  by 
the  intersection  of  vertical  B  joints  and  flat  lying  A  joints  in  the  roof. 
The  analysis  predicted  fall- in  of  rock  from  the  roof  at  100°  but  totally 
failed  in  predicting  the  major  rock  breakage  between  the  blast  wall  and  the 
roof  at  20-90°.  This  is  because  of  the  failure  of  the  geological  mapping 
to  define  the  intersecting  features  in  this  region  of  the  roof.  It  again 
calls  need  for  the  use  of  a  stratascope  or  bore  hole  camera  in  a  mapping 
program. 
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(4)  Station  0  +  90  to  1+20  CR2  and  Brill  Slot 

Fig,  3*15  is  ®  longitudinal  profile  along  the  center  line  of  CF  drift 
from  CR3  to  the  end  of  the  CR  Tee  intersection.  It  is  cleer  from  examina¬ 
tion  of  this  figure  that  a  major  rock  fall  occurred  above  CR2  and  above 
the  drill  slot.  The  drill  slot  was  totally  unsupported  and  it  is  probable 
that  this  area  progressed  towards  CR2  and  into  the  end  of  CR4.  Post  shot 
examination  revealed  flat-lying  joints  in  the  roof  above  the  drill  slot  and 
above  CR2  which  had  gone  undetected  in  the  mapping  program  because  they  do 
not  intesect  the  walls.  Geologic  cross  sections  at  station  0  +  90,  1  +  CO, 
1+10  and  1+20  are  presented  in  Pigs.  3.7>  3*8,  3»9  and  3»10>  together 
with  the  joint  influence  analysis  for  these  stations.  At  0  +  90,  the  weak¬ 
ening  effect  of  an  intersection  of  A  joints  and  B  joints  in  the  roof  is 
clearly  demonstrated.  Large  fall-out  of  rock  forming  a  cathedral  over¬ 
break  from.90-130o  is  shown  in  the  predicted  post  shot  appearance.  At 
station  1  +  00,  A  Joints  in  the  rocf  similarly  are  predicted  to  cause  fall¬ 
out  from  the  roof  between  90-110°.  The  predicted  fall-out  was,  in  fact, 
much  smaller  than  that  actually  occurring,  as  evidenced  by  the  post  shot 
cross  sections  in  Fig.  3»l6.  The  analysis  at  station  1+10  does  not  at 
all  correlate  with  the  observed  damage.  In  this  station,  as  can  be  seen 
in  the  cross  section  in  Pig.  3.16,  very  major  roof  falls  occurred.  The 
geologic  cross  section  at  station  1+20  shows  a  profusion  of  B  joints 
traversing  the  tunnel  and  a  series  of  C  joints  intersecting  the  B  joints  in 
the  roof.  The  Jointing  Is  closely  spaced  above  the  roof.  The  joint  orien¬ 
tations  are  favorable  to  progressive  stoping  upwards.  Thus,  the  predicted 
cross  section  for  1  +  20  shews  a  very  large  roof  failure  between  60-120° 
in  the  roof  as  well  as  rock  movements  inward  from  the  walls  above  the  in¬ 
vert.  This  pattern  of  fall-out  agrees  very  closely  with  that  actually 
observed  at  station  1  +  20,  Pig.  3*16. 
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(5) ‘  Station  1  +  30  and  Station  1  -f  40  CR  Tee  -  Figs.  ^Tj;i  and  3,12 

The  geology  in  section  1  +  30  is  characterized  by  a  series  of  I  shears 
on  the  blast  wall  and  a  series  of  A  joints  inclined  towards  the  blast  and 
traversing  the  tunnel.  This  pattern  of  joints  leads  to  a  predicted  rock 
breakage  diagram  which  is  not  very  consequential,  consisting  principally  of 
a  rock  fall  cathedralling  upward  between  two  shears  at  100-120°.  However, 
at  station  1  +  40,  a  series  of  B  joints  occurred  inside  the  blast  wall. 
Existence- of  the  joints  is  very  well-known  because  of  the  log  of  the  per¬ 
pendicular  leg  of  the  Tee  in  CR  north  drift.  These  B  joints  can  be  projected 
from  the  north  leg  of  the  Tee  into  the  wall  to  their  precise  positions.  They 
do  not  continue  to  1  +  30  because  they  are  offset  by  another  shear.  The 
effect  of  the  B  joints  in  the  wall  is  to  create  a  large  mass  of  rock  bounded 
by  sheared  joint  planes  which  could  be  accelerated  into  the  opening.  An 
inspection  of  the  actual  damage  at  station  1  +  40  in  the  cross  sections  of 
Fig.  3.16  reveals  that  a  large  kidney  of  broken  rock  supported  only  by  the 
wire  mesh  occurred  between  0°  and  45° .  The  cross  section  also  showed  ex¬ 
tensive  falls  from  the  lee  roof.  These  are  also  inferred  from  the  predicted 
blast  damage  of  the  joint  analysis.  In  summary,  it  appears  that  the  rock 
damage  at  1  +  40  is  very  well  predicted  by  the  joint  influence  analysis. 

(6)  Station  1  +  50  CR  Tee  Intersection 

The  stress  distribution  in  the  Tee  intersection  is  unknown  and  the 
joint  influence  diagram  method  cannot  be  applied  to  this  station. 

(7)  Station  1  +  60  -  Fig.  3.13 

The  geology  of  this  section  shows  a  series  of  A  shears  on  the  lee  wall, 
a  series  of  closely  spaced  H  shears  above  the  region  of  the  lee  wall  at  about 
120-140°,  and  a  series  of  B  joints  inside  the  blast  wall.  The  B  joints  behind 
the  blast  wall  are  a  major  weakening  factor  and  cause  inward  movement  of  rock 
material  above  and  below  0°.  In  addition,  rock  blocks  are  expected  to  drop 
out  of  the  roof  at  about  140°  because  of  the  effect  of  the  H  joints.  The 
actual  damage  at  this  station  can  be  seen  in  Fig.3«l6  .  There  was  extensive 
inward  movement  of  rock  and  failure  of  rock  bolts  on  the  blast  wall.  In 
actual  fact,  the  damage  was  much  more  extensive  than  that  predicted  by  the 
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joint  influence  analysis.  The  effect  of  the  Tee  intersection  on  the  stress 
distribution  is  not  known;  the  use  of  the  Kirsch  solution  for  the  stresses 
around  the  walls  at  this  station  is  certainly  not  correct. 

(8)  station  1  +  70  -  Fig.  3-1^ 

A,  3,  and  G  joints  compose  the  geological  section  at  1  +  70. 

A  joints  are  inclined  toward  the  weapon  point  and  primarily  intersect 
only  the  lee  wall,  lhere  are  only  two  A  joints  intersecting  the  blast 
wall.  Vertically  inclined  B  joints  occur  at  the  crown  and  invert.  !0iese 
joints  intersect  A  and  G  joints.  The  G  joints  traverse  the  tunnel  and 
are  inclined  away  from  the  weapon  point.  Joint  influence  analysis  of 
I+70  weighs  heavily  the  intersection  of  A  and  B  joints  above  the 
spring  line  on  the  lee  wall  (120°).  Intersection  A  and  G  joints  on 
the  lee  wall  below  the  spring  line  (210°)  also  produce  a  weakened  area. 
There  are  two  other  minor  failure  areas.  These  areas  occur  at  the 
intersection  of  A  and  G  joints  on  the  blast  wall  near  spring  line 
(350°)  and  B  and  G  joints  above  the  spring  line  (20°).  Rock  falls 
are  predicted  in  the  roof  at  120°  and  near  the  invert  at  210°.  Two 
minor  fallout  areas  are  also  predicted  in  the  blast  wall  at  20°  and 
350°  near  the  spring  line.  Fig.  3*l6a  has  no  section  at  1  +  70  since 
the  last  surveyed  cross  section  is  at  1  +  67.  Section  1  +  67  shows 
rock  falls  in  the  roof.  These  falls  begin  just  before  the  crown  and 
extend  to  the  lee  wall  near  the  spring  line.  Rubble  so  filled  the 
rest  of  the  tunnel  that  no  further  comparisons  could  be  made.  The 
visible  portion  of  the  tunnel  generally  conforms  with  what  was 
predicted. 


3.3  CR  NORTH  DRIFT  -  RESULTS  OP  ANALYSIS 


CR  North  drift  begins  at  the  CR  Tee  intersection  (station  i  +  52  of 
CR  Drift)  and  heads  towards  the  weapon  point.  Analyses  were  made  every 
10  feet  beginning  at  station  1  +  60  in  the  north  leg  of  CR  Tee.  The 
excavated  diameters  and  lining  types  are  shewn  in  Table  3* 5* 

The  absolute  and  relative  orientations  of  joint  sets  were  determined 
from  Pig.  3*17#  the  pre-construction  logs,  as  previously  described  for 
CR  lateral  drift.  Table  3*6  presents  the  strike  and  dip  data  for  the  dif¬ 
ferent  stations.  Table  3*7  lists  the  joint  orientation  parameters  for  con¬ 
struction  of  joint  influence  diagrams.  Joint  set  H  did  not  occur. 

In  this  drift,  the  loading  is  considerably  different  than  in  CR  lateral, 
as  the  applied  pressures  are  less  than  half  of  the  blast  pressure  and  are 
hydrostatic  (see  Table  3«l)« 

Since  this  drift  is  oriented  at  right  angles  to  CR  laterial,  the  da¬ 
maging  B  joints  of  CR  lateral  were  of  no  significance  in  CR  north.  The  I 
joints  were  also  found  to  have  no  influence  in  the  north  drift.  However, 
other  joint  sets  were  found  to  have  even  larger  influence  regions  than  in 
CR  lateral.  The  joint  influence  diagrams  for  c  *  100  psi  and  c  =  1000  psi 
are  plotted  in  Fig.  3.18a  and  b  respectively. 

(l)  Station  1  +  60  North  Leg  of  CR  Tee  -  Fig.  3.19 

The  geology  of  this  section  shows  a  series  of  B  Joints  on  the  right  wall 
and  several  joints  of  sets  C  and  G  orientations  traversing  the  tunnel.  ThL  s 
section  was  supported  by  rock  bolts.  The  Joint  influence  analysis  predicted 
breakage  of  the  rock  and  falls  from  the  right  wall.  The  actual  damage  can 
be  seen  in  Fig.  3*28,  a  long  profile  along  the  north  drift,  and  Fig.  3*29/ 
the  set  of  post  shock  cross  sections  at  varying  stations.  The  cross  section 
at  station  1  +  6l  shows  very  extensive  damage  in  the  right  wall  exactly  where 
predicted  by  the  joint  influence  analysis.  (The  right  wall,  looking  towards 
the  weapon  point,  is  0°  by  convention  in  CR  north  drift).  The  actual  extent 
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of  fall- in  of  rock  at  this  station  is  much  more  extensive  than  predicted  by 
the  analysis.  However,  the  agreement  on  location  of  damage  is  excellent. 

(2)  Station  1+70  Transition  to  Drill  Slot  -  Fig.  3«20 

This  section  is  unlined.  A  series  of  E  joints  occur  on  the  right  wall  at 
about  0°,  inclined  steeply  from  right  to  left.  Another  set  of  E  joints  cross 
the  tunnel  at  the  invert.  It  is  possible  that  these  aTe  the  same  joints  but 
have  been  offset  along. the  shears  of  set  I  that  occur  above  the  roof.  There 
are  also  a  series  of  A  joints  crossing  the  tunnel  at  wide  spacing.  The  joint 
influence  analysis  shows  that  the  intersection  of  the  A  and  E  joints  above 
the  right  wall  defines  an  unstable  block,  and  as  this  section  is  unlined,  it 
was  predicted  that  a  major  rock  fall  would  occur  from  this  position.  There 
were  no  post  shot  cross  section  at  exactly  1  +  70;  however,  in  Fig.  3»29 
there  are  cross  sections  at  1  +  67  and  1  +  66,  both  of  which  show  extensive 
movement  of  rock  from  the  roof,  more  towards  the  right  than  towards  the  left 
side. 

(3)  Station  1  +  80  CR  North  Drill  Slot  -  Fig.  3.21 

In  this  large  diameter  unlined  section,  a  series  of  G  joints  cross  the 
tunnel  in  vertical  diameter  and  a  series  of  A  joints  cross  roughly  horizont 
tally.  There  are  probably  more  A  joints  above  the  roof,  but  these  were  not 
seen.  There  is  also  an  F  joint  and  an  E  joint  in  the  right  wall.  The  ana¬ 
lysis  delimited  an  unstable  block  below  the  right  spring  line,  along  the 
intersection  of  the  two  latter  joints.  The  actual  damage  pattern  in  the 
tunnel  was  very  much  like  that  at  the  previous  station,  1  +  70,  namely,  ex¬ 
tensive  rock  falls  from  the  roof.  It  is  not  possible  to  determine  if  heave 
below  the  right  spring  line  actually  occurred  as  predicted  by  the  analysis 
as  this  was  covered  with  rubble. 

(4)  Station  1  +  90  CR  North  Drill  Slot  -  Fig.  3.22 

This  cross  section  shows  extensive  jointing  with  Q  joints  along  the 
left  wall,  a  series  of  A  joints  traversing  the  tunnel,  and  a  shear  zone  of 
orientation  I  above  the  roof.  The  intersection  of  these  features  leads  to 
an  extreme  rock  breakage  prediction,  with  very  extensive  rock  movement-  from 
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the  roof  above  the  left  wall  and  a  high  vault  above  the  right  wall.  If  the 
joints  are  considered  to  possess  great  cohesion,  then  the  major  rock  movements 
remain  only  in  the  left  wall.  The  actual  damage  at  this  location  can  be 
studied  on  the  post  shock  cross  section  at  station  1  +  92,  Fig.  3«29*  There 
is  remarkable  agreement  between  the  prediction  from  the  joint  analysis  and 
the  actually  observed  post  shock  cross  section.  It  is  quite  probable  that 
this  section  would  have  survived,  had  it  been  rock  bolted. 

(5)  Station  2  +  00  to  2  +  38  CR5  and  CR6  Drift 

The  joint  influence  analysis  for  station  2  +  00,  2  +  10,  2  +  20,  2+30 
and  2  +  38  are  presented  in  Figs.  3*23  to  3*27*  No  extensive  damage  is  pre¬ 
dicted  at  any  one  of  these  locations.  Though  small  blocks  were  delimited  by 
the  intersections  of  unstable  joints,  it  is  considered  that  the  direction  of 
the  blast  would  not  be  such  as  to  cause  a  small. volume  of  rock  to  penetrate 
the  wire  mesh  or  pull  out  the  rock  bolts.  No  very  large  volumes  of  rock  slip¬ 
ping  on  joints  were  delimited  in  the  roof  or  walls  and  it  was  believed,  in 
making  the  analysis,  that  the  wire  mesh  and  rock  bolts  would  be  sufficient 
to  handle  them.  The  post  shock  cross  sections  for  these  stations  can  be  seen 
in  Fig.  3,29*  There  was  virtually  no  distortion  or  rock  falling  in  these 
drifts,  but  the  entire  drifts  were  displaced  upwards.  There  is  no  way  that 
the  joint  influence  analysis  can  foresee  absolute  movement  of  an  entire  drift. 


3.4  DL  1  AND  2  -  RESULTS  OF  ANALYSIS 

DL  drift  is  a  lateral  drift  at  closer  range  than  CR.  Analyses  were 
made  every  10  feet  in  DL  1  and  2,  16  foot  diameter  rock  bolted  sections, 
from  station  0+50  to  0  +  80. 

Fig.  3.30  is  a  geological  log  of  the  walls  of  DL  1  and  DL  2  made  before 
the  blast  by  Corps  of  Engineers  geologists.  The  data  shown  on  the  log  were 
extrapolated  beyond  the  walls,  allowing  geologic  cross  sections  to  be  drawn 
at  stations  0  +  50,  0  +  60,  0  +  70  and  0  +  80  (Figs.  3*32  to  3»35)»  The 
joints  were  divided  into  five  sets,  A,  B,  C,  D,  E,  on  the  basis  of  orientation 
as  summarized  in  Table  3.8.  Zone  F  is  a  group  of  closely  spaced  A  joints; 
it  probably  represents  a  fault.  (The  nomenclature  of  joints  is  entirely 
different  from  that  adopted  in  the  previous  discussion  of  CR  drift.) 
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Unfortunately,  B  Joints  parallel  the  tunnel  axis  and  cannot  he  seen  in 
the  roof.  Post  shock  study  of  DL  revealed  the:  occurrence  and  significant 
role  played  by  B  Joints  in  the  roof,  particularly  above  the  unreinforced 
drill  slot  (not  analyzed  here). 

Joint  influence  regions  for  this  drift  are  plotted  in  Fig*  3«31*  Ap¬ 
plication  of  this  figure  to  the  geologic  sections  in  the  manner  previously 
described  resulted  in  the  post  shock  predictions  of  Figs.  3*32  to  3*35* 

Major  roof  falls  are  predicted  at  stations  0  +  60  and  0  +  70,  and  a  major 
rock  bolt  failure  is  predicted  on  the  lee  side  from  0  +  70  to  0  +  80. 

The  actual  damage  can  be  studied  in  Fig.  3*36,  the  po6t  shot  longitu¬ 
dinal  profile,  and  Fig.  3*37,  the  post  shot  cross  sections  in  DL  1  and  2. 

In  these  sections,  the  weapon  point  is  to  the  right;  the  post  shot  predic¬ 
tions,  unfortunately,  were  prepared  with  the  weapon  point  to  the  left.  So, 
an  inversion  of  either  figure  set  will  be  needed  to  make  a  comparison. 

Rock  movement  in  the  rock  bolted  sections  was  manifested  by:  1)  whole¬ 
sale  inward  displacement  of  seemingly  intact  walls;  2)  rock  breakage  and 
fall-out  at  the  surface,  but  sustained  by  the  chain  link  fabric;  3)  rock 
breakage  and  fall-out  not  sustained  by  the  chain  link  fabric;  and  4)  rock 
bolt  failure  with  inward  bulging  of  the  wall  and  pull-away  of  the  rock  bolt 
washer  plates  from  the  fabric. 

Failure  type  (3)  occurred  extensively  in  the  crown  from  the  drill  slot 
to  station  0  +  66,  the  vaulted  roof  reaching  an  estimated  maximum  of  15  feet 
above  the  pre-shot  crown  in  the  drill  slot.  This  may  have  been  a  progressive 
failure  propagating  into  DL  1  from  the  less  reinforced  and  larger  drill 
slot  opening.  A  and  B  joints  controlled  the  rock  breakage. 

Failure  type  (4)  occurred  on  the  lee  wall  mainly  between  0+70  and 
0  +  85  in  DL  2.  Floor  heave  occurred  extensively  throughout  the  drifts. 

The  multiple  depth  rock  deformeters  at  station  0  +  70  were  destroyed 
on  the  lee  wall  (l80°)  and  in  the  invert  (2J0°).  The  response  of  the  blast 
wall  (0°)  and  the  crown  (90°)  indicates  that  most  of  the  permanent  rock 
movement  at  these  points  occurred  within  the  first  3  feet  of  depth. 

The  predicted  post  shock  cross  sections  have  several  points  of  agree¬ 
ment  with  the  observed  deformations.  Crown  failure  in  DL  1  was  much  more 
extensive  than  the  predicted  fall-out.  This  is  particularly  true  at  station 
0  +  50,  where  the  crown  should  have  held  except  for  pyramiding  upward  along 
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the  A  and  C  joints  above  the  lee  shoulder.  The  roof  failure  nay  have  been 
a  progression  from  the  drill  slot  or  may  reflect  occurrence  of  several  B 
joints  in  the  roof;  as  previously  noted,  these  joints  parallel  the  tunnel, 
therefore  cannot  be  predicted  above  the  tunnel  by  observations  on  the 
vails.  In  cross  sections  at  0  +  70  and  0  +  80,  the  observed  movement 
of  the  lee  vail  and  bolt  failures  therein  is  grossly  matched  by  the 
predicted  cross  sections  based  on  joint  influence  analysis. 

3.5  SUMMARY  OF  RESULTS 

The  predicted  and  actual  damage  patterns  in  CR  East  West  drift,  CR 
XOrth  drift  and  DL  drift  are  compared  in  Tables  3*9,  3*i°>  and  3*H*  The 
joint  influence  analysis  predictions  of  tunnel  damage  will  not  stand  up 
under  very  detailed  comparisons  vlth  the  actual  damage.  But  in  gross  fea¬ 
tures,  the  comparison  is  good.  It  must  be  noted  that  the  analysis  is  very 
sensitive  to  the  quality  of  geological  information.  The  geological  data 
presented  are  very  complete  in  comparison  vith  usual  geological  studies 
performed  for  other  types  of  vorks.  However,  considerably  more  detail  about 
the  locations  and  continuity  of  individual  joints  and  shears  is  required 
to  refine  the  quality  of  the  results. 
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Table  3.1 

APPLIED  PRESSURES  TO  DEVELOP 
STRESS  DISTRIBUTIONS  USED  IN  JOINT  ANALYSIS 


Applied  Pressures 

Sections 

Radial  -  0° 

Radial  -  90° 

Longitudinal 

DL  and  CR  (EW) 

P 

(v/l-v)P 

(V/1-V)P 

CR  (North) 

(v/l-v)P 

(v/l-v)P 
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P  ■  peak  pressure  of  blast 


v 
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CR  (EW)  DRIFT  -  LINING  TYPE  AND  TUNNEL  DIMENSIONS 
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Table  3-3 
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Drift  CR  -  Lateral  Joint  Set  Strikes  and  D’ps  at  Each  Station 


Strike 


"A” 

"B" 

"C" 

"D" 

"E" 

"p" 

"G” 

"H" 

»tjl» 

0‘ 20 

H53W(12) 

N62W(13)  N32E(»0 

H72E(2) 

- 

- 

- 

- 

H!i017(l)i 

0-30 

H50W(5) 

n68w(7) 

M31E(3) 

i/i9E(*0 

Ni9E(i) 

N72E(3) 

N' 5E(l) 

- 

- 

- 

Ktl,7(<i)i 

0i-'i0 

W50W(9) 

M67W(5) 

N32E(1) 

H';0E(6) 

- 

- 

- 

N!i1V/(3) 

0+50 

Hl-'fW(3) 

N7vW(3) 

- 

N1i6e(6) 

N38v/(l) 

- 

- 

N29'2(l) 

0f60 

N*“0V/(8) 

K78w(»i) 

N23E(3) 

HC5E(3) 

Hli6E(6) 

N  38'.7(i) 

- 

- 

KGV7  (1) 

0^70 

Nil'*/ (7) 

N79W(‘i) 

N13E(5) 

n8oe(i) 

N49E(9) 

173817  (1) 

- 

- 

- 

0r80 

nWwC  ) 

H73W(7) 

N13E(3) 

- 

nV>e(*;) 

- 

N35E(!«) 

-• 

” 

0+90 

N5017  ( 5 ) 

K73W(iO) 

- 

- 

iMe(8) 

- 

- 

- 

- 

1.00 
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;i7hw(2) 
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- 
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• 

M  3E(2) 

- 

- 

]  >-10 

N+9-'(8) 
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- 

n;.ve(2) 

- 

H'f8w(:;) 

- 

f 

“  ! 

1+90 

N  9W(*  ) 

m6»iw(8) 

*K’2E(3) 
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Table  3-J* 


CR  -  Lateral  -  Joint  Orientation  Parameters 


Joint 

Set 

<P 

D° 

0? 

6° 

(13)  B 

-3 

83 

88.0 

-87 

(7) 

-9 

83 

82.8 

-8l.l 

(5) 

-3 

90 

90.0 

-87.0 

(3) 

-11 

-85 

-8*1.8 

79.0 

00 

-1*4 

79 

73.7 

-76.3 

00 

-15 

78 

77.6 

-75.4 

(7) 

-9 

8h 

83.8 

-8l.l 

do) 

-9 

88 

88.0 

-81.5 

(2) 

-10 

90 

90.0 

-80.0 

(2) 

-16 

90 

90.0 

-7*f.O 

(16) 

-13 

89 

89.O 

-77.0 

(1) 

-15 

90 

90.0 

-75.1 

-11 

90 

90.0 

-79.0 

0) 

-it 

-88 

-39.0 

86.0 

(7) 

-3 

-87 

-87.0 

82.0 

(10) 

-2 

-8*4 

-8*4.0 

88.1 

(9) 

+3 

90 

90.0 

87.0 

(3) 

+4 

+88 

8*4.0 

36.0 

(3) 

•7 

90 

90.0 

83.0 

(6) 

-3 

79 

79.0 

-87.0 

00 

-5 

88 

88.0 

-85.0 

(-25) 

+86.*: 

81.3 

-55.2 

83.7 

\  Joint 

!  set 

0° 

D° 

<P 

6° 

CO  C 

-89 

-77 

-*i.2 

13 

(3) 

90 

-62 

•0 

-28 

00 

83 

-66 

-15.2 

-25 

(1) 

83 

-66 

-15.2 

-25 

(3) 

85 

-76 

-19.3 

-15' 

CO 

77 

-68 

-29.0 

-25.5 

(3) 

77 

-66 

-26.7 

-27.3 

* 

00 

7*i 

-7*4 

-*41.8 

-22.5 

(12) 

76 

-75 

-*>1.9 

-20.5 

(9) 

78 

-67 

-25.9 

-25.7 

Cd 

3l 

-59 

-1*4.*' 

-32.2 

(3) 

80 

-*i8 

-10.8 

-*43.0 

(2) 

77 

-50 

-1*1.9 

-**1.7 

(11) 

89 

-*i6 

-1.0 

-*t4.0 

(56) 

-10.2 

-40.2 

00 

-.0  71  0 
.0  71.0 


3.0  76.’. 


Table  3-4  Continued 


CR-Lateral 


6°  $ 


0*20  G 

- 

- 

- 

- 

0-30 

- 

- 

- 

- 

0-40 

- 

- 

- 

- 

m 

0+50 

- 

- 

- 

- 

o-6o 

- 

- 

- 

- 

(8) 

0+70 

-8i4 

65 

*12.6 

-25.9 

(1*0  o<8o 

-8i 

67 

+20*3 

-24  6 

(2) 

0+90 

- 

- 

- 

- 

3-00 

-75 

81* 

68.1 

-16.!* 

1+10 

- 

- 

(2) 

1»20 

-71 

75 

50.5 

-2*1 . 1 

1+30 

- 

- 

- 

- 

(0 

1-5)0 

-79 

6^> 

22.2 

-27.5 

1  +50 

- 

- 

- 

~ 

0) 

J  .60 

-90 

61' 

0 

-25.0 

(?) 

1+70 

-6o 

68 

51.2 

-26.7 

(•) 

l"8o 

-64 

55 

32.0 

-'2.8 

(■) 

1+90 

-66 

51 

26.5 

-1*5.0 

(.) 

P+00 

-70 

60 

30.6 

-35.7 

0) 

2J10 

- 

- 

- 

- 

2-19 

-79 

59 

17.6 

-33-1 

(  o) 

Mean  or 

Range : 

30.0 

-l>0.0 

Joint 

Set 

0° 

D° 

0° 

6° 

H 

- 

- 
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- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

(1) 

15 

-65 

-6l!.2 

-76.5 

- 

• 

• 

- 

- 

- 

- 

- 

- 

- 

- 

- 

(1) 

-7 

-65 

-64.8 

83.7 

(1) 

-28 

-59 

-65-8 

23.0 

(3) 

('*) 

-28 

-59 

-55-8 

-65 

23.O 

80 

Note:  Joints  cease  to  have  any  weakening  effect  if  their  angle  with  the 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  problems  for  analysis  posed  by  a  jointed  hard  rock  mass  are  not 
insurmountable.  Provided  that  sufficient  geological  and  engineering  data 
about  the  rock  can  be  garnered,  it  is  possible  to  estimate  the  behavior  of 
an  excavation  under  given  loading  conditions. 

The  Piledriver  test  has  demonstrated  conclusively  that  the  pattern  of 
rock  breakage  and  failures  of  rock  bolted  and  unlined  underground  galleries 
at  intermediate  range  depends  primarily  on  geological  conditions.  Since 
these  conditions  exist  before  the  blast,  there  is  no  a  priori  reason  why 
they  cannot  be  discerned  with  sufficient  detail  to  define  the  geological 
conditions  sufficiently  well  for  analysis.  The  only  real  limitation  is 
economic.  (It  is  another  question  to  define  the  loading  conditions  with 
certainty. ) 

With  the  development  of  bore  hole  television  and  photography  devices, 
the  expense  of  very  detailed  geological  mapping  should  not  be  excessive 
for  a  gallery  serving  a  valuable  purpose.  It  is  proposed  that  in  a  future 
project,  an  effort  be  made  to  make  use  of  these  devices  so  that  geological 
cross  sections  can  be  prepared  definitively  showing  the  location  of  joints 
and  fractures  to  a  distance  of  at  least  one  diameter  around  the  tunnel. 

Methods  are  becoming  available  for  determining  the  strength  and  de- 
formability  properties  of  weakness  planes  in  rock.  Some  of  the  available 
techniques  of  sampling  and  testing  discontinuities  in  rock  are  reviewed  in 
Section  2.4.  In  a  future  project,  after  the  geological  study  has  offered 
a  classification  of  the  weakness  surfaces,  grouping  discontinuities  toge¬ 
ther  that  are  alike  mechanically,  it  is  proposed  that  a  program  of  sampling 
and  testing  be  conducted  to  determine  the  characteristic  strength  and  stiff¬ 
ness  of  each  group. 

It  is  believed  based  upon  the  work  of  this  project,  that  an  analysis 
premised  upon  information  of  the  quality  proposed  above  would  be  trust¬ 
worthy.  Therefore,  it  should  be  possible  to  design  reinforcement  for  an 
opening  making  optimum  use  of  the  reinforcing  materials.  In  particular, 
rock  bolt  patterns  can  be  unbalanced  to  strengthen  the  section  in  accor¬ 
dance  with  localized  needs. 
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APPENDIX  1 


FINITE  ELEMENT  ANALYSIS 


COMPUTER  PROGRAM 


University  of  California,  Berkeley 
College  of  Engineering 


170 


COMPUTER  PROGRAM  FOR  ANALYSIS  OF  JOINTED  ROCK  SYSTEMS* 

WITH  RESIDUAL  STRESS  AND  ACCELERATION  OPTIONS 
-  PLANE  STRAIN  -  NONLINEAR  ALONG  TANGENTS  -  FINITE  STRAIN  - 


PURPOSE 


The  purpose  of  this  computer  program  is  to  determine  displacements 
throughout  a  jointed  rock  body,  stresses  within  reck  blocks,  and  normal 
and  tangential  stresses  on  rock  joints.  Tne  j..ints  have  no  tensile  strong! 
and  a  finite  shear  strength.  The  blocks  are  orthotropic  or  isotropic, 
linear,  elastic  solids.  A  failure  criterion  is  included,  for  the  rock: 


oct 


=  N  +  D  a 


oct 


B 


and  for  the  joints: 

s  =  c  +  cr  tan(|> 

For  tensile  values  of  sigocta'nedral,  failure  criterion  in  the  rock  is 
simply: 


°jnin  CTt 


=  uniaxial  tensile  strength) 


INPUT  DATA 

A.  IDENTIFICATION  CARD  -  (72H) 


Columns  1  -  72  of  this  card  contain  information  to  be  printed  with 
results . 


B.  CONTROL  CARD  -  (415,  2F10.2,  6l5) 


Columns  1  -  5 

6-10 
11  -  15 
16  -  20 
21  -  30 
31  -  40 
4i  -  45 

r 

46  -  50 
51  -  55 


Number  of  nods.],  points  (900  maximum) 

Number  of  elements  (670  maximum) 

Number  of  different  materials  (12  maximum) 

Number  of  boundary  pressure  cards  (200  maximum) 
Acceleration  in  X  direction 
Acceleration  in  Y  direction 

Number  of  approximations  (never  less  than  2  for 
jointed  body) 

Code  for  running  multiple  problems  (0  if  second 

problem  follows,  1  if  no  second  problem  follows) 
Residual  stress  code  (reads  residual  stress  if  not 
equal  to  0) 


¥  Adapted  from  plane  stress  program  by  E.  L.  Wilson,  July  1965,  and  joint 
stiffness  and  joint  stress  subroutines  by  R.  L.  Taylor  and  R.  £.  Goodman, 

July  1967.  2/68 

7/68  Revision 
9/68  Revision 
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5 6  -  60  Joint  cut-off  number:  ell  materials  vith  higher 
numbers  are  joints 

6l  -  65  Code  to  read  punched  deck  output  (effective  strain 
deck)  to  restart  problem  at  point  of  termination 
of  last  run  (0  =  no  deck  read;  1  «  deck  read)* 

66  -  70  Number  of  joint  elements  (hOO  maximum) 

C.  MATERIAL  PROPERTY  CARES 

The  following  group  of  cards  must  be  supplied  for  each  different  solid 
material  and  joint  material.  Hie  last  numbered  solid  material  defines  the 
solid's  properties  after  tension  or  shear  failure. 

Solid  Materials 

First  Card  -  (15,  F10.0) 

Columns  1  -  5  Materials  identification  (any  number  from  1  to  12) 

6  -  15  Mass  density  of  materials 


Second  Card  -  (7F10.0) 


Columns  1-10 

Tensile  strength  (input  as  positive  quantity) 

11-20 

E  compression 

21  -  30 

V 

31  -  hO 

E  tension 

hi  -  50 

N 

51-60 

61  -  70 

?}  Toct,-  n  *  N  +  D  a  B 

B  (failure)  oct 

Joint  Materials 

First  Card  -  (15) 

Columns  1  -  5  Number  of  Joint  materials 

Second  Card  -  (7F10.0) 

Columns  1-10  Normal  stiffness,  pef/ft 

11-20  Tangential  stiffness,  pef/ft 
21  -  30  Joint  cohesion,  psf 


*If  number  of  joint  materials  was  0,  the  previous  run  will  not  have  yielded 
punched  output. 
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31  -  4o  Joint,  friction  angle ,  degrees 
4i  -  50  Tangential  stiffness  after  slip  occurs 
51  -  60  Joint  closure  allowable  (ft.)j  input  as  negative 
quantity 

MOTE:  If  joint  is  non-stiff,  the  value  of  stiffness  is  probably  l/lO  -  l/lOO 
of  the  stiffness  of  the  intervening  block,  e.  g.  ^block  -  1*0  x  lOo  psf, 
block  thickness  =  5<0  ft,  ttjen  block  stiffness  is  2.0  x  107  psf/ft  and  a 
low  joint  stiffness  {t  10°  psf/ft. 

■D.  NODAL  POINT  CARDS  -  (15,  F5-0,  4F10.0) 


One  card  for  each  nodal  point  with  the  following  information.  Joint 
elements  are  obtained  by  double  rows  of  nodal  points  at  the  same  coordinates. 


Columns 


1  -  5 

Nodal  point  number 

6  -10 

Number  which  indicates  if  displacements  or  forces 
are  to  be  specified 

11  -20 

X  -  ordinate 

21  -30 

Y  -  ordinate 

31  -4o 

XR 

hi  -50 

XZ 

NOTE:  If  the  number  in  column  10  is 

0  XR  is  the  specified  X-load  and  XZ  is  the  specified  Y-load. 

1  XR  is  the  specified  X-aisplacement  and  XZ  is  the  specified  Y-load. 

2  XR  is  the  specified  X-load  and  XZ  is  the  specified  Y-displacement. 

3  XR  is  the  specified  X-displacement  and  XZ  is  the  specified  Y-displacement. 


All  loads  are  considered  to  be  total  forces  acting  on  an  element  of  unit 
thickness.  Nodal  point  cards  must  be  in  numerical  sequence.  If  cards  are 
omitted,  the  emitted  nodal  points  are  generated  at  equal  intervals  along  a 
straight  line  between  the  defined  nodal  points.  The  necessary  temperatures 
are  determined  by  linear  interpolation.  The  boundary  code  (column  10),  XR, 
and  XZ  are  set  equal  to  zero. 


E.  ELEMENT  CARDS  -  (615) 


Solid  Elements 


One  card  for  each 


element 


1,2 


Columns  1  -  5 

6-10 
11  -  15 
16  -  20 
21  -  25 
26  -  30 


Element 

Nodal  Point  I 

Nodal  Point  J 

Nodal  Point  K 

Nodal  Point  L 

Material  identification 


For  a  right-hand  coordinate  system,  order  nodal  points  counter-clockwise  around. 
2Maximum  difference  between  nodal  point  I.D.  must  be  less  than  27. 
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Element  cards  must  be  In  element  number  sequence.  If  element  cards  are 
omitted,  the  program  automatically  generates  the  emitted  information  by  incre¬ 
menting  by  one  the  preceding  I,  J,  K,  and  L.  The  material  identification  code 
for  the  generated  cards  is  set  equal  to  the  value  of  the  last  card.  The  lest 
element  card  must  always  be  supplied. 

Triangular  elements  are  also  permissible;  they  are  identified  by  repeating 
the  last  nodal  point  number  (i.e.  I,  J,  K,  K).  One-dimensional  bar  elements 
are  identified  by  a  nodal  point  numbering  sequence  of  the  form  I,  J,  J,  I. 


Joint  Elements 


Nodal  Point  Naming  convention  for  joint  elements: 


Nodal  points  must  be  numbered  I,  J,  K,  L  proceeding  from  bottom  left  corner 
counterclockwise  to  upper  left  corner.  Bottom  and  top  are  defined  by  x’,  y ' 
system  of  coordinates  created  by  rotation  6  (  l8oo)  fro*  x  to  x',  where  x' 
is  along  the  length  of  the  element. 


F,  PRESSURE  CARDS  -  (215,  XFIO.O) 


One  card  for  each  boundary  element  which  is  subjected  to  a  normal  pressure. 


Columns 


As  shewn  above,  the  boundary  element  must  be  on  the  left  as  one  progresses 
from  I  to  J.  Surface  tensile  force  is  input  as  a  negative  pressure.  Joints 
cannot  be  placed  on  boundaries. 


G.  INCREMENTAL  LOADING  (15F5*2) 


1 


m 


One  card  giving  multiples  of  input  pressure  to  be  applied  successively 
and  cumulated. 


Example ; 

For  a  problem  with  3  pressure  increments  equal  to  1000  psi,  4000 
psi,  and  2000  psi  (total  pressure:  TOGO  psi),  with  1000  psi  applied 
in  F,  the  following  card  would  be  required: 

Column  5  1 

Column  10  4 

Column  15  2 


NOTE:  The  number  of  entries  on  this  card  must  equal  the  numoer  cf  approx¬ 
imations  indicated  on  the  control  card. 

H.  SKEW  BOUNDARIES 

If  the  number  in  columns  5  -  10  of  the  nodal  poin'  card',  is  /t.xw  than  0, 
1,  2,  or  3,  it  is  interpreted  as  the  magnitude  of  an  ar.;;!e  in  decrees. 
This  angle  is  shown  below. 


The  terms  in  columns  31  -  50  of  the  nodal  point  card  are  then  interpreted 
as  follows: 

XR  is  the  specified  load  in  the  s-direciion 

XZ  is  the  specified  displacement  in  the  n-direction 
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The  eagle  6  oust  always  be  input  ea  e  negative  angle  and  my  range 
frc*  >*001  to  -l80  degrees.  Hence,  +  l7<5  degree  is  -the  earn  as  >179*0 
degrees.  The  displacements  of  these  nodal  points  which  are  printed 
by  the  prograa  are: 

u^  *  the  dlsplacesent  in  the  s -direction 
u^.  *  the  displacement  in  the  n-dlrection 

I.  RESIDUAL  STRESS  CARDS 

If  residual  stress  code  on  the  control  card  is  not  equal  to  0,  one 
card  for  each  eleaent  must  be  supplied  with  the  following  information; 

Columns  L  -  5  Element  mother 
6-7  Blank 
8-22  x-stress 
23  -  37  y-stress 
38  -  52  xy- stress 

If  cards  are  emitted,  the  stresses  for  the  emitted  elements  are  assumed 
equal  to  those  supplied  on  the  last  card.  If  all  elements  have  the  same 
stresses,  only  the  1st  and  last  element  stresses  my  be  supplied. 

If  eleaent  is  a  Joint 

Columns  1  -  5  Element  number 
6-7  Blank 

8-22  Initial  tangential  stress 
23  -  37  Initial  normal  stress 

J.  EFFECTIVE  STRAIN  CARDS 

These  are  punched  by  computer  after  last  cycle  of  a  run.  If  asked 
for  by  1  in  column  65  of  the  control  card,  the  effective  strain 
deck  is  read  if  placed  behind  the  last  data  card  in  the  data  deck. 


OUTPUT  INFORMATION 

Hie  following  information  is  developed  and  printed  by  the  program: 

1.  Reprint  of  input  data 

2.  Nodal  point  displacements*,  for  the  preceding  cycle  and  cumulative 
3*  Stresses  at  the  center  of  each  solid  element*;  principle 

stresses  are  cumulative 

4.  Stresses  and  relative  displacements  (of  center  top  and 
bottom)  of  each  Joint  element*,  cumulative 


*2.,  3. ,  and  4.  are  printed  after  each  cycle. 
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NOTE: 


SIGMAX  (maximum  tension) 


New  element  corner  positions  are  used  to  calculate 
new  stiffnesses  on  each  cycle  (finite  strain) 


17T 


*  P  P  C  N  D  I  X  J 


FORTRAN  PROGRAM  FOR  I HCR£M£NTAL  LOADING  Of  PLANE  STRAIN  STRUCTURES 
PROPERTIES  Of  SOL  10  HATFULS  ARE  NODI  PICO  If  STRESSES  INDICATE  fAUURt 
,»  n  *  *°Wf*  LA*  ,M  Tf*MS  0f  OCTAHEDRAL  STRE5SE5 

UP  TO  400  JOINT  CLEMENTS  CAN  8E  INTRODUCED  BANDWIDTH  •  |0 

WRITTEN  f OR  USE  ON  CDC  *400.  PROGRAM  LENGTH  120000  OCTAL 

ADAPTED  FROM  PLANE  STRESS  ROUT  J  NC  Of  PROf  C  WILSON. 


PROGRAM  MAIN  UNPUT, OUTPUT*  TAOf*.  TAPE9.  PuNCMi 

c  plane  strain  with  joints,  cumulative  loading,  non  linear  along  tangcn 

COMMON  NUMNP.NUMEL«NUHMAT.NUHPC#  ACfLX.ACCLY,  N , VOL  » T £ MP *  Ml Yp£ , 

1  M£OU?l*  Cl*. 12),  ROU?)*XXNNfl?)»Rf900|.2l9PO).URl9OO1.u2l900l. 

2  COOC«90OJ,!6C(?0DJ.jACJ?e0J.PA(?0O».ANGLCU).$|GI)O».LBAD 
COMMON  /A«G/  RR«<S>.Z7f «5».S«10.10).P110»,U(*»,lmU),DDI30). 

I  hh<6*10>.RR<4).22U),C<4*4)»h<6,]0)»0<6*6)*FIS.1C).TP(6I*X|(10) 
2»ECt?l .1  *<670. SI  .CPSi670t.MTAGl4?0»,RSTRSU)  .RES  lp(  670,  31  * 

1  NRES.PREACT,NUHJT,S|GStl3.670» 

COMHON/6ANA«G/M0ANO.NKMei.»:feJ  jfcO».AU60.#0) 

COMMON/ JOl  NT/N$hEll*<COO£ I 400  J  .fNST  UCO)  ,f  TSI  UOL)  1*00)  » 

2  CPSN$T(430).  PPPlil.NNN.NJ 
DIMENSION  BST  < 10071  *  FACTORIISI 

C 

C  READ  ANO  PRINT  Of  rONt«0L  IMEOkmaTION  AND  MATERIAL  PROPERTIES 

load-o 

SO  READ  1000. MED. NUMNf.WJMCi.*NUMMAT.NUMPC»ACELX,ACtLr.  NP.  NEND.NReS» 

1  N  SHELL  .NSTAGE  .N'JNJT 

PRINT  2000.  MfD.NUMMP.NUMEL.NUNMAT.NUMPC.ACCLX.ACCLr.NUMJT.NP 
C 

If  INSHCLL  *E0.  0)  NSMELL  •  13 
SR  00  59  M-l.NUMMAT 

REAO  1001.  MTTPC.  ROIMTVPE) 

If  CMTYPE  »GT •  N SHELL )  GO  TO  560 
PRINT  2011.  mtype.  roimtypei 
READ  1005.  IEIJ»MTYPEI.J-1.7) 

PRINT  2010.  ICU.MTYPf  |,J*1,7) 

GO  TO  39 

5*0  READ  1005.  I E I J. MTYPE ) , J* 1  ,6 ) 

PRINT  2017,  MTVPC 

PRINT  2016. IE  I J*MIYPE ) » J* 1 ,*> 

59  CONTINUE 
C 

. . . . . . . . . . . . 

C  READ  ANO  PRINT  Of  NOOAL  POINT  DATA 


*3  L*l*l 

JflN-ll  40*41.42 

42  DO  46  1*1,3 

46  RCSIDI  l,  !  )*R£SID(L<*1«! ) 

GO  TO  43 

41  PRINT  1007.  It.  IRES1DIX* I ) • I *1 *3  >.  K*NL»K  1 
IflNUMCL-NI  40.45.47 
40  PRINT  1008 *N 
LOAD  «  1 

43  CONTINUE 


C  REAO  OUTPUT  Of  PREVIOUS  RUN  TO  RESTART  PROBLEM 

. . . . . 

DO  312  N*1 iNUMJT 

312  KCDOEIN)  >0 
KNP*2*NUMNP 
00  313  K*1*RNP 

313  BSTIO.O.O 
IflNSTAGE.EO.O)  GO  TO  314 
PRINT  1017 

READ  1018*  (N.MTAG<N)*CPSfN).SIGST<l*NJ.SIGSTf?,Nt.S!6ST(3«N)< 

1  N*1.HUNEL) 

PRINT  10U.  (N.MTaGI  N)  .CPSlN).SlGST(lfN)  .SIGST  t  ?.N)  .SI  GST  1 3  »N)  • 

1  N'l.NUMEL) 

PRINT  1019 

READ  1020»^CODf  !Nl.fPSNSTCN).EPS?ST|N|.fNSTIN|.fTST<N).N*l.NUMjn 
PRINT  1020* I XCOOE  1*0 »CPSNSTlN),fPSTSTINl«fNST<N|.fTSTlNI*N*l .NUMJT I 
PRINT  1021 

READ  1022  *  CN.DST 1 ?»N-1 l»BSTt?»H)»N*l .NUMRP) 

PRINT  102? » I N. BST I 2*N- l )»8ST I 2*N) .N*1*NUMNP| 

. . . . . . . 

C  DETERMINE  BAND  WIOTH 

. . . . . . 

314  J*0 

00  340  N>l*NUM£L 

DO  340  I *1 #4 

DO  323  L*1 .4 

XK*!ABS< IXIN.I l-tXtN.l) I 

If IAK.LE.39)GO  TO  313 

PRINT  20U»  N 

LBAD*LBAD*1 

313  IFIKK-J)  325.323*320 
320  J«XX 
323  CONTINUE 
340  CONTINUE 
*BANO*2»J*2 
PRINT  2012*M6ANO 
PRINT  2019.L6A0 
If  IL6AD.NE.  OlSTOP 

C  * SOLVE  NON-LINEAR  STRUCTURE  BY  SUCCESSIVE  APPROXIMATIONS 


PRINT  2004 
L-0 


ff(KRES.NE.O)  00  TO  34? 
00  341  NM.NUMFL 


60  READ  1002*  N.COOE IN ) *R I N) *Z IN) *UP <  N ) «u£ IN) 

NL*L* 1 
2X-N-L 

0R*IR(K)~R  IlM/ZX 
DZ* I Z I N> -Z I L ) ) /2X 
70  l*LM 

If IN-L )  100*90*00 
00  COOEIL ) *0.0 
RIL)*RIL-I)40R 
zuwil-iwdz 

UR(L)*0.0 
UZIU-0.0 
GO  TO  70 

90  PRINT  2002*  IK.CODCU  )  *RCK»  ,ZU  )  *URU|  ,U2  U)  ,K*NL.NI 
IFINUMNP-NJ  100*110*40 
100  PRINT  2009.  N 
LBAD  •  1 
GO  TO  60 
110  CONTINUE 

. . . . . . . 

C  READ  AND  PRINT  Of  ELEMENT  PROPERTIES 

. . . . . 

PRINT  2001 

N*0 

130  READ  1003*  H*I1X(M*1)*1«1,3) 

•  140  N*N*1 

If  IM-N)  170*170,130 
150  IX<N,l)*lX|N-l.llU 
1XIN»2)-1X(N-1*?U1 
IXIN,3)*IX|N-1*3)*1 
|X(N,4)*tX(N-t.4)M 
1XIN,3>»IX1N-1.M 

170  PRINT  2003.  N*UXIN. 11*1*1.31 

If  IN-N)  110*100*140 
110  IF  (NUMEL-N)  190*190*1)0 
190  CONTINUE 

. . . . . 

C  REAO  ANO  PRINT  Of  PRCSSltfC  BOUNOARy  CONDITIONS 

. . . . . . . 

IF  INUMPC.EO.  0)  GO  TO  310 
290  PRINT  2005 

DO  300  LO.NUNPC 
REAO  1004.  IBCIU.JBCIL)*PRIL) 

I  300  PRINT  2007.  I6CIL I  * J6C I L ) ,PRtL I 

PRINT  2029 

READ  20)0*  If AC TOR  1 11*1*1* NP I 
PRINT  2031,  I  FAC TOR 1 1 )*!•!# NP) 

310  CONTINUE 

,  c  . . •••••••••••• 

C  REAO  AND  PRINT  Of  INITIAL  OATA  FOR  THE  PROBLEM 

C  . . ••••»•#•••••« . . 

!  IFINRES.CO.O)  00  TO  43 

PRINT  1006 
L*0 

t  4T  READ  1007.  N.  (RCSIOIN. I ) *  1 • 1 . 3 » 


DO  341  I«l.) 

RESIDIN* I ) *0*0 
3*1  CONTINUE 
34?  CONTINUE 

IFINSTAGE.NE.O!  GO  10  300 
00  330  N-l.NUMEL 
00  307  1*1,3 
307  SI0STII*N»*0.0 
MTAGINI*! 

350  EPSINI*0*0 

DO  360  N*1,NUNJT 

epsnstini*o. 

epststin>*o* 

FN5T(N)*0. 

FTSTIN) *0* 

360  CONTINUE 
370  XNP*2»NUMNP 
DO  373  XM.XNP 
373  BSTIK) *0*0 
C 

300  DO  300  NNN-I. NP 
C  FORM  STIFFNESS  MATRIX 
C 

00  40)  XNM.NUMNP 
RIKN)*RlKN)*6Sf I2*N-1 ) 

403  ZUN)*2IRN)*0STI2*RNI 
PRf ACT "FACTOR  I NNN) 

CALL  STIFF 
C 

C  SOLVE  FOR  DlSPLACFMfNTS 

C 

CALL  8ANS0L 
00  400  X*1,RNP 
400  BSTUI*BU)4BSTU) 

C 

PRINT  iOOO*  NNN 
PRINT  2000 

PRINT  2006*  IN.BI2*N-l)*dl2«N),BSTI2»N-n *R$T(?«NI»N*1*NUMNP| 

C 

C  COMPUTE  STRESSES 

c 

CALL  stress 
CALL  JTSTR 
NRES  •  0 
NSTAGE  *•  0 

300  CONTINUE 

IFINUNJT.EO.O)  GO  TO  301 

PUNCH  1010.  IN.N»AGfNI*EPSlNI.SIGSTIl.N|.SlGS7l2.N).SIGSTl3»Nl. 

1  N.l.NUMEL) 

PUNCM1020.UCODElN)*fPSNSIIN).fPSTSTIN).FNST(N»*FTSTIN),N*l,NUHJI) 
PUNCH  1022  *<N»8$tl?*N-l) *8$T l?«N) . N* 1 .NUMNPJ 

301  CONTINUE 
C 

. . . . . . it . . . . 

IF  |  NCND  .CO*  0  )  GO  7C  30 
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. . . 

1000  F0RNATU2A4/4l9.2F10.?,6l9> 

1001  F ORMAT  (  19.  HO. 01 

1002  format  m.  f9.o*  af io.oi 

1001  FORMAT  (619) 

1004  FORMAT  I219.F10.0I 
1009  FORMAT  <  7F10.0J 

1006  FORMAT  (29H1 INITIAL  STRESSES  J H  R 0CX//6H  ELEMENT  4X  9H  X-STRESS  6X 
19H  Y-STRCSS  6X  10K  Xy-S7RfSS//l 

1007  FORMAT  (t9»2X.3E19.4) 

1006  FORMAT  I91H  RES10UAL  STRESS  INPUT  ERROR  N*  151 

1017  FORNATHO  ELEMENT  *TAG  EPS  SIG  1  SIG  2  SIG3 

1  •! 

10U  FORMAT  (2110*4£12.)> 

101?  FORMAT  (•  XCOOE  EPSNST  EPSTST  FNST  FT$T«) 

1020  FORMAT  U9.4E10.2) 

1021  FORMAT (  •  *>  C*JN  X  OISP  CUM  y  01SP  *1 

1022  F0RMATII5.2E20.2J 

2000  FORMAT  UH1  12A6/ 

1  10M0  NUMBER  OF  NOOAL  POINTS - II  / 

2  5 OHO  NUMBER  OF  ElfmFNTS— — - -  19  / 

3  10MO  NUMRER  OF  01FF.  MATERIALS -  13  / 

4  10HO  NUMRER  OF  PRESSURE  CAROS——  IS  / 

9  50M0  X-ACCELCRATION - - - £12.4/ 

4  3 OHO  ^-ACCELERATION--— - £12.4/ 

7  10 HO  NUMBER  OF  JOINT  ELEMENTS  —  11/ 

4  10HO  NUMBER  OF  APPROXIMATIONS -  111 

2001  FORMAT  147M1ELEMENT  NO.  I  J  X  L 

2002  FORMAT  m2.  F  12,2.  ?F12.3.  2E24.7J 
2001  FORMAT  11119.416.111?) 

2004  FORMAT  UORH1NOOAL  POINT  TYPE  X-OROINATE 

1  AD  OR  DISPLACEMENT  Y  LOAD  OR  DISPLACEMENT 

2005  FORMAT  I29M0PRESSURE  BOUNDARY  CONDITIONS/  ?4M 
1URE  ) 

2004  FORMAT  m2.4E20.7l 
2007  FORMAT  (2I6.F12.3J 
2001  FORMA TMONOOAl  POINT  X  DISPLACEMENT  Y  DISPLACEMENT  CUM 

1  X  displacement  CUM  r  oisplacfmcnim 

2009  FORMAT  (26HONOOAL  POINT  CARD  ERROR  N*  111 

2010  FORMAT  U5M0TENS.  STRENGTHS  5HEICI  9X  6HNU  1 IX  4MCITI 

1  10X  6M  N  10X  4H  D  10X  4H  8  /7E15.5J 

2011  FORMAT  (17H0MATERIAL  NUMBER*  13.  15H.  MASS  DENSITY*  C12.4J 

2012  FORMAT  ( /6HM0ANO*  15  /J 

2014  FORMAT  ( *0  XN  XT  C  PMl 

1  XT  FINAL  NAX.  CLOSURE  •/  6E15.4J 

2017  FORMAT  (1H010X.16H  MATERIAL  NUM6ERI51 
2011  FORMAT  <24H  ELEMENT  CARD  ERROR  N*  1 4 J 
2019  FORMAT  (•  LOAD  ••  15J 

2029  FORMAT  (*0  MULTIPLES  OF  INPUT  PRESSURE  SUCCESSIVELY  APPLIED  •  / I 

2030  FORMAT ( 15F3.2  I 

2031  FORMAT  I15F7.2J 

3000  FORMAT  (•!  INCREMENT  N0.«I2J 


Y-ORDINATE  X  LO 
t 

I  J  PRESS 


IF  I VOL.LE.  O.J  LBAD  •  1 
IF  (VOL  .GT.  0,0)  GO  TO  144 
PRINT  2003.  N 

144  IF  ( 1XIN.3J.E0.IXIN.4J J  GO  TO  145 

145  00  150  !!*1«* 

CC*Sfti.lOI/St 10.101 
P(!1I*P1IU-CC*P(10I 
00  150  JJ*1 .9 

150  sm.jj»*sm .jjj-cctsiio.jji 
c 

DO  140  11*1.6 
CC*S(II.9J/SI9»9J 
p( 1 1  J*P( 1 1 J-CC*PI9J 
00  160  JJ* 1 .8 

160  sm»jjj*sm.jjj-cc*s(9.jj> 
c 

C  ADD  ELEMENT  STIFFNESS  TO  TOTAL  STIFFNESS 

c 

165  00  166  1*1.4 

166  LM(!  )>2MX(N»I  1-2 
C 

DO  200  t*1.4 
00  200  K*l»2 
!l*lMt|)«K-XSH!FT 
XX*2«I-2*X 
81 !ll*BI II J4PIKX1 
DO  200  J*1 .4 
00  200  L*l»2 
JJ*LM|  JJ*L-|U1-X5H|FT 
LL*2»J-2*L 
IFIJJJ  200.200*175 
174  IFINO-JJ)  180.195.195 
180  PRINT  2004.  N 
L8A0  •  1 
GO  TO  210 

194  Am.JJI.A(!!.JJ)*SIKX.LU 
200  CONTINUE 
210  CONTINUE 
C 

C  ADO  CONCENTRATED  FORCES  WITHIN  6L0CX 

C 

00  250  H*NL.NM 
X*2»N-XSH|FT 
BIX J*8f X )*U2tN) 

250  B(X-1)*8(X-1J4URINJ 

c 

C  BOUNOARy  CONDITIONS 

c 

C  1.  PRESSURE  6.C. 

C 

IF  (NTA*PC J  2 60.310.260 
260  DO  300  L-l.NUMPC 
I-IBCfO 
J*  JSCILJ 

PP*(PR(Lt/2.J*PRFACT 


SUBROUTINE  STIFF 

c 

COMMON  NUMNP.WWEL.NUMMAT.NUMPC.  acelx.acely.  N.VOt .TEMP.nTYPC. 

1  WE0I12).  El  8  *1 2 J  *  ROI 12  J .XXNNl 1?J. 0(9091 .ZI900J . UR 1 900 J  »UZ 19001 • 

2  C00€ 1 900) . I6C ( 2001 • JNCI 200 J . PR 1 200 J .ANGLE ILJ.SIGl 10 J »LBAD 
COMMON  /ARG/  RRR(5).Z22l3).S(10.10l.PI]0).TT(4J.LM(4J.DDI3.3). 

1  JRK4.10I.RRI4J.2Z(4J.C(4.4).HI6.10).0(6.6J.Fl6.10J.TP(6J.Xl(10l 
2*£E 171 .IXC6T0. 5 J.E*SI670).MTAG( 6701 *RSIRS(4|.R£SID<670. 31* 

J  NRES.PRFACT.NUMJT.SIG5TI 3.4701 
CONMON/BANARG/MBANO.HUM3LK .01 160) »A( 160*80) 

COMMON/ J01NT/NSHEU  .XCOOEUOOJ  .FNST  (400)  .FIST  (400)  .EPSTST  (4001  • 

2  EPSNSTI400).  PPPI81.NNN.NJ 
C 

. . . . . . 

C  INITIALIZATION 

(••••••»»••• . . . . . . . . 

REWIND  9 
N6  •  40 
NO*2«N0 
N02*2«N0 
LB AO  •  0 
NLMBLX*0 
NJ»0 
C 

00  50  N*1 »N02 
6(41*0.0 
00  50  M*l,NO 
90  AIN. M) *0.0 

. . . . •••••••• . ••••••»••••••••••*< 

C  FORM  STIFFNESS  MATRIX  IN  0LOCXS 

. . 

80  NUMBLX*NUM0LX*: 

NH«N8* I NUM0LK  + 1 ) 

nm*nm-nb 

nl*nm-nb*i 

X$MIFT*2»NL-2 

C 

00  210  N*t»NUN€L 
C 

IF  UXCN.9JI  210.210.49 
69  00  40  1*1.4 

IF  I IXIN. I l-NL )  80*70.70 
*  70  IF  UXIN.IJ-NN)  90.90.80 

•0  CONTINUE 
GO  TO  210 
C 

90  IF  I IXIN. 9)  .LE.  N SHELL )  GO  TO  92 
CALL  JTSTIF 

IF  IVOL.GT.  0.01  GO  TO  189 
PRINT  2003. N 
1 8  AO  •  2 
GO  TO  189 
92  CALL  QUAD 

I XIN.9  J*-IX(N#9) 

C 


Dz»izm-zu))*pp 

OR*(RIJ)-Rmi«PP 
264  !I-?*I-XSHIFT 
JJ*2«J-X SHIFT 
IF  (II)  280.280.269 

269  IF  III-NO)  270.270.280 

270  SINA.0.0 
COSA-l.O 

IF  (COOFMM  271  .272.272 

271  SINArSlN ICOOE (1)1 

ccsA.cosicooEin  i 

27?  SIII-ll*BIII-|UIC0SA*D2«SINA«r)R| 

5(111*0(11 )- l SINA»OZ-COS*#OR I 
280  IF  (JJ)  300.300.269 

289  IF  IJJ-ND)  290.290*300 

290  SINA.0.0 
COSA.l.O 

IF  (COOEIJ)l  291.292*292 

291  SINA*SINIC00E IJI) 

COSA-COSICOOflJM 

2 92  B( JJ-ll-BI JJ-! )*ICOSA*DZ*SINA*D» - 
0( JJ)*0( JJ1-ISINA»OZ-COSA«DR» 

300  CONTINUE 
C 

C  2.  DISPLACEMENT  0.C. 

C 

310  00  400*N*NL*NH 

I F  (M-NUMNP)  319.319.400 
319  0*UR(M| 

N*2*M-1-XSH|FT 
IF  (COOE(MJ)  390.400.918 

316  IF  CCOOE«N)-l.)  317,370.117 

317  IF  (COOEIN)-?.)  318.390.318 
316  IF  (CODE CMI»3. )  390.380.390 
370  CALL  MOOIFy ( A .0.NO2.M4AND .N.U) 

GO  TO  400 

910  CALL  MOOIFY(A*0,NO?.mbanO.N,UI 
390  U-UZ(M) 

K*N*1 

CALL  MOOIFYIA.8.NO2.M0ANO.N.U) 

400  CONTINUE 
C 

C  WRITE  BIOCX  OF  E(XjAT IONS  ON  TAPE  ANO  SHIFT  L'P  LOWER  0LOCX 
C 

WRITE  19)  (0IN)»(AtN.H).M*l.M0ANO) »N*I .NO) 

C 

00  420  N*1 ,N0 
X*N*«lD 
0(N)*B(K ) 

01X1*0.0 
00  420  N*1.N0 
AIN.M)*A(X.M) 

420  AIX,N)*0*0 
C 

C  CHECX  FOR  LAST  BLOCK 
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IF  INM-NU*NPl  60«*60.*80 

*•0  continue 


IF  ILflAO  .fO.  01  00  !0  300 

PRINT  ZOOS.  LOAD 

STOP 

300  *£Tl*N 

2003  FOSMt  <26wONCGAIIVC  AREA  ELEMENT  NO.  Ml 
200*  FORMAT  I?9hO^ANO  WIOTM  rxCFfOS  «UOHABtC  |4| 
2003  FO*NAT  «  •  l  A  AO  *  *IM 
fNO 


SUBROUTINE  TRISTFISI , JJ.KCI 
C 

COMMON  NONNP.NUMCt  •N1^»**AT.NU*‘PC«  ACC  L  *  •  ACELY  •  •»  »  VOL  .  Tf  NP.MT  1  *[  « 

1  HC0I17I*  E  *  8  » 17  * »  ROri?l.XXNNtl?).»C900|.ZC9O'n.U»l9r0|.UZC900l. 
?  COOF* 9001 •  I9C(?OO).j8CC20OI.PRC200>.AN0lFMI.<  IGl 10 1.10*0 
COMMON  /A«0/  RSRISl.ZZZC  3)  .SUOflO*  .PI  10»  .1114)  .LWUl.DPC  3.3T* 

1  MM<6»10I.RRI*1.ZZI*»  .CU.Al.Hlfc.lO>  .0« 6 .6 1  *F  *6  *  10 1  •  TP C 6 )  *X»  I  SOI 
?.££(?>  .IXI6T0.S)  .f  P$C  670 1. M7AGI 6705  ♦RSTRSM)  .  Rt 510*670.3) . 

1  NRCS.PRFACT.N'j-iJf.SIOSM  3.6701 
COMMGN/bANARG/MBANO.NC'WBLK.eUSOl  .A<  180*80) 

COMMON/JOlNTFNSMftl.rfOOEUOOl  .FNSTUOOl  #F  * M  U^C I  *f  PMST  1*00 I  ♦ 

2  fPSNST ( *001 •  PPPI81.HNN.NJ 

C 

C  1.  INITIALIZATION 

C 

1*111-11 

LMI2I-JJ 

LM|3l«KK 

C 

RAUHARRUIT 

RR* ?|«RRR( JJT 
RROl-RRRlKKI 
RRHI-RRR* 1 1 1 

*  rzui-zzzun 
ZZ ( 2 1 -ZZZ ( JJ) 

Z2I31-ZZZIKK) 

zzui-zzzim 

c 

83  00  100  1-1*6 
00  90  J-1.10 
FM.JMO.O 
90  HU. J1-J.0 
00  100  J-1.6 
100  DU.J5-0.0 

c 

C  3.  FORM  INTCGRAllGlTUCmOi 

c 

COMM-RRC7l»lZZl3)-ZZll )  I  .RRCI  l  •  C  ZZ  *2 1 -ZZ  1 3 11  MR  ( 3 1  •  *Z?  C 1  l-ZZ  C  2 1  > 


129  FI  I .JI-FI I « J) «Q! I.KMMfK.J) 

130  CONTINUE 
C 

00  HO  1-1.10 

00  HO  K-l  .6 

IF  IHIKtll)  Dl.HC.nS 

131  00  139  J-1.10 

139  SI I.Jl-SII »JI.M*K»!I#F*K.J| 

HO  CONTINUE 
C 

C  FOR*  R£S IWJAl  STRESS  CONTRIBUTIONS  TO  T*F  LOAD* 
C 

ISO  DO  160  1-1.10 

160  PUHPU:-*IC1HHI?,M«I  -RSTRSUl  1 

1  ♦  XH1HHI6.IH*  -RSTRSI2I  1 

?  ♦Xlll»*W3.ll*l  -RSTRSI3I  1 

3  ♦XIUHHIS.II**  -RSTRSI3)  I 

C 

C  ACCELERATION  LOAO* 

C 

CO**-RO*MTYPCl*Xl 1 11/3.0 
00  170  1-1.3 
J-2*L*UI-1 
PlJI-PlJ|-ACfLX*COHM 
170  PlJUT-PI JMT-ACCLY-COMM 

c 

c  FORM  STRAIN  TRANSFORMATION  MATRIX 

C 

*00  00  *10  1-1.6 

00  *10  J-1.10 

*10  MHt  I.JI-MMf  I.JI.HU.J) 

C 

♦IS  RETURN 
C 

END 


SUBROUTINE  QUAD 
C 

COMMON  NUMNP.NUH£L»NUM*AT »NU*PC»  ACCLX.ACELY.  N.VOL .TE*P.MTYP£. 

1  MEDI12).  Eie.171.  R0I12I.XXNNU2)  .RI90C1 .ZI9PC5 .URI90P5 .UZ*900> » 

2  COOEI90OI  ♦  IBCI 200).  J*CI  2001  .PR*  200)  .ANGLE  •*  1  .MCI  101  .LBAD 
COMMON  /ARG/  RRR«S).ZZZr3l.Sll0,10I.PI  105. TTUi.LmU), POO. 31. 

1  HH(6.10I.RRt*5.ZZl*).C(*.*T«HI6«l0l.0(6.65.F*6.10)»TP<6l>X!ll0l 
2.CCI7)  .IX 1670. ST *CPSI670 1.MT AG* 6701  .RSTRSU5.PES  10*670. 35. 

3  NRCS.PRF ACT .NumJT .SlCST *3.6705 

COMMON /0  AN  ARG/M9ANO.M1M0LX.B  1180 1.  A*  160*801 

COMMON/JOINT/NSMELL*XCOOE UOO) .FNST**005 ,FTS7<*001 .rPSTSTI*00| , 

2  EPSNST I *001 •  PPP* 8 1 .NNN.NJ 

c 

l-IXIN.l) 

J-  |X( N.2  5 


Xllll-COMM/2.0 
VOL-VOL*X I C 1 1 

c 

OI2.2l.XMlHCCl.il 
0C2.61-XM1HCU.21 
0C3.31-XM1  HC(*.*I 
0*3.SI-XM1HC(*.*1 
0*9tS)-XMlHC**.*l 
016.61-XM  1 5*<*2»2  5 
C 

108  00  110  1-1.6 
DO  110  J-1.6 

no  o* j.i i-o*  i  •  jt 

c 

C  ♦.  FORM  COEFFICIENT-OISPLACfMfNT  TRANSFORMATION  m*TR|X 

C 

DO<I.l>-«RR*2  5*ZZ*3l-RR*3l-Z/*2>l/COM'- 
001 1.21-*RR*35*ZZI15-RR* 1 I -Z2*  31  |/C(XM 

00* 1 . 3 1 • I RR « 1 1 *Z Z *  2 » -RR « 2 1 -ZZ *1 1 5 /CO^m 
00* 2. 1 1 -  * Zl t 2 1 -ZZ ( 3 1 1 /CO-m 
OOC2.Z)-*ZZ*31-ZZ*  IM/COMM 

DOC  2.31.  C ZZ*  1 1-ZZI2M/C0-M 
OOI3.15-«RR*31-RRI2M/COMM 

00*3*21-1 RR*1 1-RRI 3* I /CO** 

001 3. 3 *• * RR* 2 5-RM Cl  I I/COmm 
C 

00  120  1-1.9 
J-2*L*< 15-1 
Hll.J5-00n.l> 

H(2.J5-OD*2.l I 
HI3.J5-DOI3.tl 

HI*. J.I 5  -  0011*1) 

M*  S* J*1 5  *001 2  *  1 1 
120  HI  6. J* 1 1 -00*  3.1) 

C 

C  ROTATE  UNKNOWNS  |F  RFQUIPCO 
C 

00  12»  2-1.2 
I-LM* Jl 

IF  lANGLEUIl  122U2S.12S 
122  S1NA.SIN* ANGLE (III 
cosa.cosc anglci m 
I J-2*l 

00  12*  K-1.6 
TEM.HIK.IJ-I) 

HtK*tJ-ll«T£M*C0SA«HfK»IJ5-S  NA 
12*  HIK.  1  Jl-  -T£M*$INA*M|K.I  JHCOSA 
12S  CONTINUE 
C 

C  S.  FORM  ELEMENT  STIFFNESS  MATRIX  (HITMOIXm) 

C 

00  130  J-1.10 
00  130  K-1.6 
IF  (HIK.JII  128.130. 128 
128  00  12*  1-1*6 


K-IXIN.3I 

L-tXiN.*) 

NTfPC-lXIN.S) 

C 

c  FORM  STRESS  STRAIN  RELATIONSHIP  FOR  PLANE  STRAIN 
DO  10S  KK-1.7 
105  CEIKKMf  IKK.l  .NTYPEI 
C 
C 

IF  IMTAGIN1-2I  82. 8*. 83 

82  CE135-EEUI 
GO  TO  8* 

83  ECU  5  -Cf  *  31 
C 

8*  €EUI-EEU5/l*l.-EC*2l***2l 
CE  C  31 -CE *  3T/M1.-EC*  2) )  **2 1 
CCI2I-CC*25/U.-EFI2I) 

C 

XX.EFUI/EEI3I 

CCMM-EEU  l/IXX-Ef  *2H-21 

CU»11-C0MM*XX 

C* l»? I -COMM-EE I  2  I 

CU.  31-0.0 

CI2.1HCI1.2T 

CI2.25-COMM 

Cl  2.31-0.0 
CI1.1 T-0.0 

Cl  3. 21-0.0 

C*3.31-.S*Ef* I l/IXX*Ff*25) 

C 

SS-SINIf P$IN* I 
CC-COS*  EPS*  N) 1 
S2»SS*SS 
C2-CC*CC 
SC*SS*CC 
C 

OC  1.1I-S2 
OI1.25-C? 

0(1. 31-0.0 
01 1  .*  5 --SC 
OI2.11-C2 
OI2.25-S2 
0*  2*3 T -0.0 
0I2HI-SC 
013*11-2. »$C 
0*3.21—013.11 
01 3. 31-0.0 
Ol  3.*I—C2*S2 
C 

DO  88  11-1.3 
DO  88  JJ-!.* 

HIM.  JJ  1.0.0 
00  88  K»>1 ,3 

M(  1 1 .  JJI  •  «(  1 1  *  JJ|  «CU  I  .KKHDIKK.JJI 
88  CONTINUE 


l8o 


DO  If  11*1*4 
DO  99  JJ«l*4 
Cl  II  •  JJUO.Q 
DO  19  KK*1*3 

Cl  Il.JJucm  .jJt.Diu.il 
19  CONTINUE 
C 

. . . «MiM...*.»MM«*....**»**..*t...M* 

C  REPLACE  RESIDUAL  STRESSES  FOR  Nt«  CLCmC N1  fit  RSTRSI  I J) 

. . . 

C 

mwtES.CQ.o>  go  to  n? 

DO  ill  SJM.9 

111  RSTRS  HJURESIOU.IJ) 

GO  TO  114 
11?  DO  113  IJ-1.3 

113  RSTRSUJUO. 

114  RSTR$(4)*RSTRS( 51 

c 

C  FORM  OUADRILA7CR4L  STIFFNESS  MATRIX 

C 

ft**())MR< I  )«ft| J|*ft|K|«R(Ll 1/4.0 
ZZZI5>*IZCI>*ZIJ)*ZlK)*ZlLI)/4.0 
DO  94  M* 1 *4 
NM*|X(N.M) 

9)  RRfclMURIMM) 

94  ZZZIMUZIWO 
C 

DO  100  11*1*10 
PI  111*0*0 
DO  95  JJ*1,6 
99  MM( JJ* I ( I *0*0 
DO  100  JJ*1«10 
ioo  sim.jjuo.o 

00  119  11*1*4 
JJ«iXIN*t I 1 

119  AHGL€lin*C00eiJJI/5T.9 
C 

C  FORM  BAR  STIFFNESS 
C 

IF  I |X<N*2>-!X(N*3)1  ?90»?40*?90 
240  TTIll»-eCI4l 
7TI2I— ECU) 

DR-RIJI-RIII 

DZ*ZtJWII) 

XL*SORT(OR**2*DZ**2) 

RRRI3)*f RI 1 )+Rt J) t /2»-2.*EEU)*OZ/*L 
ZZZISUlZIIUZlJU/r.^.'EEUHOR/XL 
CALL  TRISTF 11*2*5) 

GO  TO  190 
290  CONTINUE 

c 

VOL  «  0.0 

IF  (K*NC*  LI  GO  TO  125 


CALL  QUAD 
C 

DO  120  1*1*4 
II*?*! 

JJ*2* I XI Nt  1 1 
P|  II-IUBI JJ-l) 

120  Pilll-blJJl 
C 

00  150  1*1.2 
RRt I »*PI 1*41 
00  150  K*1.6 

150  RRIII*RRII»-SM*9.<»*PIX» 

C 

COKM.St 9.91*51 10. *01-519. 10) *SI »0*9» 

IF  I COMM )  145 *  160*195 

155  PI  9 1* I  SI 10»10)*RRl 1 1-519. 101 *RR( 2 1 1 /CO** 
PI  101*1-51 10<9)*RRm*5l9*9)*RRf  2M/C0MM 
C 

160  00  170  1*1*6 
TPI 1 1 *0.0 
DO  170  K-l.10 

170  TPI  1|*TPI  I  UMMI  I  ,o«P|U 
C 
C 

RRI ll-TPI?) 

RR(2)*TPI6) 

RRl 3S  *0.0 
RRl4)*TP»3l*TPl5» 

C 

176  00  160  1*1*4 
51GI I  URSTRSI 1 1 
RSTRSI 11*0.0 
DO  160  K*l*4 

160  SIGH  l*$!G(l)*C(I.KMRRI<l 
51001*51014) 

00  190  1*1*9 
SIGUUSIGH  )*SIG5Tf  |*mi 
190  5IG5T I l*MI*5IGl I ) 


C  OUTPUT  STRESSES 

. . * . < 

c 

C  CALCULATE  PR!NCIP*L  STRESSES 

C 

CC*ISIGI !l+SlG(?)l/?.0 

BB*(5tGtH -510(2)1/2.0 

CR*SORTCSB**2*5IGO)**2) 

$IG(4)*CC«CR 

S1GI5)*CC-CR 

IFI66.NE.  0.)  GO  TO  195 

EP5IN)*77777. 

GO  TO  196 

199  EPSlN)*A»AN2lSIGl7I.B9l/2. 

196  510161*57. 396*EPSIR) 


120  CALL  TRISTFI 1 *2*3) 

RRRI5l*IRRRl  1  l*RRRI?)*RRRt  3M/3.0 
ZZ2(9)*IZZZIl!*ZZZl2)*ZZZI3))/3.0 
VOL*XtU> 

IF  IVOL.GT*  0.1  RETURN 
PRINT  1000.M 
RETURN 
125  CONTINUE 

CALL  TR! $TF 14*1*9) 

CALL  TRISTFI 1 *2*9) 

CALL  TRISTFI2.3.5) 

CALL  TRISTFI 9*4*5) 

C 

IF IVOL.GT .  9*)  GO  TO  126 
PRINT  1000.N 
124  DO  140  11*1*6 
DO  140  JJ*1.10 
MWI  II.JJUMMI 1 1 » JJ) /4.0 
140  CONTINUE 
C 

150  RETURN 
C 

1000  FORMAT  I  •  NEGATIVE  AREA,  ELEMENT  NO.*  15) 
ENO 


SUBROUTINE  STRESS 
C 

COMMON  NOM*P*NUM£L*NUMMAT*NUMPC*  ACELX.ACELYr  N.VOL *T£MP,MTYPE • 

1  HfDI  12)  »  £16.12).  ROI12I.XXNNU2I  .Rl  900)  .Z 1900)  .URI900 )  .UZ 1 900 1 . 

2  C00C1900 1  *!BC  1 200).  JBC « 200). PR  I  200).  ANGLE  (4 1.SIGMOI.LBAD 
COMMON  /ARG/  PRRI5)»ZZZI5).5(10.10).PI  10).TM4I  .LMU).00(3.)|. 

1  HKI4. 10). RRI 4 ).2Z (4 ).C 1 4 .4) .HI  6. 10) .016.6 ).E 16 *10). TPI 6)*X 1 1 10) 
2*ECI7) .1X1470 .91. EPS! 670) .MTAGI 6701 .RSTRSI 4 | .RCSI0I67O. 9) » 

9  NRES* PREACT. NlPUT.SIGSr*3. 670) 

C v^40N/ BAN ARG/ MBAND* NUM0LK.fi 1 160) .A| 160*60) 

COMMON/ JOInT/NSmClL*KCOO£I 400) .ENST 1400) .ETST UOO) .EPSTST (4001 • 

2  EPSNSTI400)*  PPPlll.NNN.NJ 

. . . . . . * . . * 

C  COMPUTE  CLEMENT  STRESSES 

. . . 

mprint  *0 

C 

DO  900  M*1*«UM£L 
C 

nfail  •  o 

N*M 

I  X I N.5 )• I A8SI IXIN.9I) 

MTVPE*1XIN*51 

C 

IP  (MTTPE.  GT.  NSMELLI  GO  TO  900 

oo  m  x*i  »t 

Ml  ECIX)*EIX*1*MTYP£| 


r 

MTAGI N ) • 1 

IFISIGI4))  261.261.260 

260  MTAGIN)*? 

261  IFISIGI5))  269.269.26/ 

262  MTAGIN)*3 

263  CONTINUE 
C 

C  ••«**.*.*».*l.****«»(«t*»l***t**,M*tC.«***t**l.*...*M.*«..*«. 

c  APPLY  failure  CRITERION  TAUOCT  .Lf.  N  *  D*S!GOCT**b 

C  . . . . * . . . * . 

C 

SIGT  •  fll.MTYPE! 

SlGm*EE(?>*ISIGf4)*S!G(SM 
SIGOCT  «-(S!G(4l«S!G<9  USIGI7)  1/3*0 
IFISIGOCT.GE.  0.0)  GO  TO  290 
IE  1SIGI4).LT.  SIGT)  GO  TO  104 
t  XIN.5 ) »NSM£LL 
NFAIL  *2 
GO  TO  104 

290  TAUOCT  -  SORTIHSIGU>-SIGl3>l**2> 

1  ♦  ClMGl*9l-SIG(7l|**2> 

2  ♦! ISIGl 7)-SIGl4) )**2) )/9.0 
CONN  •  ECU) 

COND  •  EEI5) 

C0N6  •  EEI6) 

STREN  •  CONN  ♦  COnD*IS!GOCT**CONB) 

IEI TAUOCT. LT.STREN)  GO  TO  |0* 

IXIN.9)*NSHELL 
NEAIL  •  1 
C 

104  IE  IMPRINT)  110.109.110 

105  PRINT  2000 
HPR!nT«50 

110  MPRInT*MPR|NT-1 

c 

909  PRINT  2001.  N.RRRl  5 >  »2ZZ I  9 )  •  I  SIGN  ) .  1*1  .6)  *NEA  |L 
300  CONTINUE 
C 

320  RETURN 

c 

2000  FORMAT  I7H1EL.no.  7X  ]HX  7X  JHY  4X  6MX-5tRf SS  49  6mY-S*RESS  3X 
19MX*-STRESS  ?X  l\»MAX-SfRESS  2X  10MMIN-STRESS  7h  ANGLE 

1  2X  1 7M  FAILCO  IF  l  OR  2  1 

2001  FORMAT  I I7t2F6.2*lP9C!2.4.GPlF7.2.llO) 

ENO 


SUBROUTINE  MODIFY! A.fi.Nf O.MfiAND.N.Ul 

c 

DIMENSION  Al  140*60)  .  BlUOt 
C 

DO  259  M-2.M0ANO 


l8l 


trio  235.235. 730 

710  8tK»*8UI-A«*.N|RU 
AfK.MLO.O 
211  K-N.M-1 

IFfNfO-Kl  210.240*240 
2*0  6m*BlKI-AtN.M)*U 
A 1 N.N | «0.0 
250  COKTtHUC 
AINtl 1*1*0 
e«n»-u 
return 

END 


SUBROUTINE  BANSOl 
COH*ON/BANARG/M«,nunNI«:.BI  J40» 

XX-  90 
nl*nn*i 
Nm»NN*NN 
rewind  $ 

RfVlNO  9 
N9*C 

00  TO  150 


REDUCE  EQUATIONS  BY  BLOCKS 


1.  SHIFT  BLOCK  Or  EQUATIONS 

100  NB*NB*1 

00  125  N*I.XN 

XM«hm*X 

6(MI*B(AN) 

BlNM>*0.0 
00  125  m.I.mm 
AfN.*)*A(NM.M) 

125  A(XM#M|«0.0 


2.  Rf ad  xtxT  block  or  equations  ixto  core 

IF  CXUMOLK-XBl  ISO. 700.150 
150  RCAO  «  9 1  (Bfrtl  *CA|N.M»,M*1 .MM) .N*NL.NH) 

IF  (MB)  2C0.100.700 

3.  REDUCE  BLOCK  OF  EOLATIONS 

700  00  ICO  N*l.NN 

IF  (AIX.111  725.300. 725 
225  6lXI*eiX»/AIX.l) 


?  COOC <  *001  •  IBC 1 200>*J*C« 700 1  .PRf  7001  .*NGLf  141.? TGUOI.lBaD 
C0»*<N  /ARC/  RRR<5l.222f5l.5l|O.IOl.P«  10).TTi4i,t«Ul.0Pf  3.31. 

I  *MI4. 10 1  »RR« 4 1 .22 1 4  I. C » *.4 1 .Mff.iOi .Of  6.6 l.F 16.10 ).»P1 61.* 1 1I0> 
2.ff  m.IKliTO.Sl  .FPSC6T0)  .HIAM6T0I. RStR6f4l.Rf5lOI6TO.il. 

1  XRCS.RRFACT.XOHJT.SICSM 5.6TOI 

CO»»n)X' JOIXl/N$HtLL  .KCOOf  ‘  4001. FxSTf  <,091. F  1ST  UOO)»r«>STST,*OOl. 

7  £B$xST(400>»  PPP«8t .NNN.XJ 

DIMENSION  EST If  I 10.101 

COOI VALfXCC  ISI10.10I.CST IFflO.lOll 

dimension  Af4.4i.tRi 7.71 

data  a/?.  .1*  .-t**-?*.  I  *•?•  •■'?•  ••■1..-I* *7**1  •***?♦  **1  **1# *2./ 

RfAl  KS.KX.L 
NJ*MJ»I 

c 

C  REPLACE  RESIDUAL  STRESSES  for  nTk  ELEMENT  BY  REMRSIIJI 

C  . . . . . 

c 

Iff NRCS.E0.G1  00  TO  112 

00  111  u*i.i 

111  RSTRS  (IJMRESIOIX.I Jl 
GO  TO  114 

112  00  113  IJ-1.3 
11)  RS1RSIIJ»*0. 

114  RSTRSf4|.RSTRSl3l 
C 

116  |f*IXfN.ll 
•  JJ*  IXlN.21 

OR*Rf JJI-Rl 1 1 1 
02*2 1  JJ>-2 (III 
L- SCOT  I OR* 0» *02 *02  I 
VOL  •  L 

IFIL.EO.0.1  GO  TO  7CI 
MAT*! XIX. 51 
lX«X.5>*-!XfX.5f 
C**  MATERIAL  PROPERTIES 
KX*E f 1  .MAI  I 

IFtKCOOEIXJ  Gl.O)  60  TO  40 
KS*Ef  2.MAT . 

GO  TO  41 

40  KS*Cf5.MATi 

41  lFlF«sTlNJI.Lr.0.lG0  TO  50 
KX*0. 

KS-O. 

50  !rte»S*SlfNJ‘.GT.  ff6.MATI)G0  TO  704 
KX* 1000.*KX 
704  CONTINUE 

C0MS»KS*L/4. 

C0HX.rM*L/6. 

C  INITIAL I2E 

00  100  1 1  •  1 .6 
Pf  111*0.0 
00  100  JJM.8 
100  ESTIFt If. Jj 1*0.0 


00  775  L*7 »mm 
IF  ufN.UI  730.7T5.7T0 
730  C*Af  N.L  * /A I N. 1 1 
l*N*L-l 
J*0 

00  250  K-L.MM 

J*J*I 

750  Af I.J»*Af I .J)-C*A|N.KI 
Of  I  1*61  I  1  •  A  f  N  *L  I  *R  I H  1 
AfX.LI'C 
775  CONTINUE 
300  CONTINUE 
C 

C  4.  NRITf  BLOCK  OF  REOuCEO  EQUATIONS  OX  T*Pr  7 
C 

IF  fNUMBLK-NBl  375.400.375 
175  VR I Tf  fei  l6fX).fAIX.M).M*7»MM),X.l.XX» 

GO  TO  100 

C  BACK-SUBSTITUTION 

. . . . . 

AOO  00  450  M.| ,xx 
N*NN*1-M 
00  425  K*7»mm 
L*N*K-1 

475  BfNI*BfN>'AfX.KI*0(LI 
NM*N-»NN 
BfXM) *BIN> 

450  A(NM,NDI*glNI 
XB*NB-1 

IF  (RBI  475.500.475 
475  BACKSPACE  8 

READ  115  l8lNI.IAtN.N»,M.2.MMUN*l.NX> 
BACKSPACE  A 
GO  TO  400 

. . ..........I*. . 

C  ORDER  UNKNOWNS  IX  N  ARRAY 

. . . . . . . 

500  K'O 

00  600  XB*I.XUM0lk 
00  600  X*1 .XX 
NH.X.NN 
K*K*  1 

600  Bl K l*A f NM.XBl 

C 

RETURN 

C 

END 


00  160  1*1.4 
J*?*l-  1 

PPPf JI«L*RSTRSI 11/7. 

J*7»I 

160  PPPfJI*L*RST«S«?W?. 

00  161  1*1.4 

161  PPPf I  I .-PPPi I ) 

00  200  11-1.4 
IS-?*II-I 
IN*?*! I 

00  200  J  J* 1.4 
JS  •  2*JJ-1 
JN  •  2*JJ 

EST IF  IIS.JSI  •  C0MS*A1I|.JJ» 

700  ESTIff I X . JN i • CONN* A  f I  I » J  J> 

C  rotate  TO  GLOBAL  COORDINATES 
TRII.l»*OR/L 

Iff  TRf  1  .1  P*t  J.l.  )  RtT'lRN 

TRf 1.71*02/1 

TRf 7.1l*-02/L 

YR*?.?)*DR/L 

00  400  XX* I .4 

00  410  11*1.6 

JJ.  ?*XX-i 

TfMP  •  ESM*  < 1 I.JJ) 

00  410  KKM.7 

f STIFI I  I .JJ)*TEMP*TRf 1 .KKI.f ST |FC I  I .7*NN**TRf7#KK) 

410  JJ*JJ*1 

D<"»  4?0  M*1.8 
J J*2*NN- I 

TfMP  •ESTIFfJJ.il  I 
00  470  KK* I . 7 

ESTIFf  JJ.J  N.  IRt  l.KK»*lfMP.TRf  7*KK  **PS»  If  »  7*XS.  I  i  » 
470  JJ*JJ* 1 
490  CONTINUE 

no  40 i  1*1.4 
J«7*l-I 
I  1*2*1 

Pf  J)  •  PPPI  •  TRi  7.7  i  .PP.  iMi*»w',i. 

401  Pf I  I > *PPPf 1 1 • TRt l ,7 > «» P<  >•  l  ••*  » 

Rf  TURN 

70 J  PRINT  it 4' ,  h 
«f  TURN 

2000  »04M*T«  |*f  « 

7090  fO*MAT,j4^  -  N  .N 

f  xr 


SUBROUTINE  JTS1IF 

COPOX  NUMNP.NUMfL  .XUM-AT  .xu^P^  ,  Alft«.ACfi*. 

I  HCPI17*.  (.8*17*.  RO.  17*.XXNN»  f>. •••'«./  •  4x 


1  HHI4»l0).RRt41.22<4|.C<4.4).H<6.10>«D<&»6>.ri6«10>.TP<6>»Xt<10) 
2«CCI?)*lX<470.SI.CPSC470ltNTAGf4?0>»9STRSUl.RCSI0t47ru3l. 

J  NACS.PRFACT.NUNJT.S1GSTI3.4T0I 
CONMON/WNARG/MtANO.NUNFLK.RIHOl  *AI  160*101 

CONNON/J01NT/NSHCLL.KCOOEU03)  *FNSTI40*H.FTSTI400I  .FPSTSTUOOl . 

2  EPSNSTUOOI.  WllhNIMiNj 

DIMENSION  VUtiUUI 
REAL  L.KN.KS 
M»0 

PAINT  1001 

C  «••*••»£ $TA8l.ISH  OlSPLACF**fNTS  ALONG  ANO  NORMAL  .TO  JOINT 
DO  500  N*l.WJMCL 
MATMXfN.S) 

iriMAT.tC.NSHCLLI  <*0  TO  300 

FTSTtN)*FTSY(M)*RCS10IN. 1 ) 

FNST(N|«rNSTIMt«ftES!0tH»2> 

00  10  X*2 
10  RESIDING  UO. 

(•IXIN.1I 
J«  IXIN.2I 

02«2(Jl-2lll 
ftJ«0.3»tR(J>*Rtl  )l 
2J*0.3M2lJ|42<ln 
L*SORTIOA*OA*D2*D2» 

VOL  •  L 
ON -ON /l 
02*02/1. 

OO  100  21*1*4 
K*|X(N»1 ( I  4 

VC1 I »»-OI2*K-l 1*02*01 2*K)«0R 
100  UUII*8t2«K-ll*OR*ei2*Ki*OZ 
C  *•••••«•  CONFUTE  EFFECTIVE  STRAIN 
XN*0« 

KS*0«0 

IFICPSNSf(M)«G7«0. I  GO  TO  200 
KN*E’1.NAT> 

IFIEFSNSTIMI.GT.EI4.HATI »G0  TO  20} 

KN*1000.«KN 

20}  IFIKCOOCIMI.CT.OI  CO  TO  201 
KS*  C«2»HATI 
GO  TO  200 
201  KS  •  Cl  5. MAT  I 

*00  EP5T*0.3*IUt41-UI 1 ) *U( 3 1 -Of  2 ) ) 

EPSIN1*Q*3*(V(4)~VI 1 1 *713) -VI 2 1 ) 

C  ••••••••••••••••M.MHMIN* 

c  C*^»CM!J*?*i*lIIS*?,rAR  F0^Cr  PCR  0,1,7  teilG™  AW  CA*-CUIATC  strength 

FN*KN«EF$INI*FNST|HI 
205  FT  •  KS*E?ST  *FTSTIM| 

FTST|HI*FT 

FNST|H»*FN 

epsnstini-epsinmepsnstini 


CPSTSTIMI*£PST  ♦  EPSTSTIM* 

CJ  •  CI3.NATJ 

PHI J*  C| 4. HAT  I 

PHIR*FHIJ*}.14139144/1I0. 

SIREN  •  CJ*ABS(FN>«SIN(FHIR>/C0SIFHIR| 

KCOOCINLO 

IF  (ABSIFTI.LT.SMENI  GO  TO  210 
Kcooc<m*i 

210  PAINT  1000.N.RJ.2J.FN.FT.CPSINI.CPST.EPSNSTIHI.EPSTSTIHI.KCOOEIHI 
300  CONTINUE 
RETURN 

1001  FORMAT! *1  CLEM  X  V  NORMAL  STRESS  TANG  STRESS  NOR 

INAL  OISP  TANG  DISP  CUM  NORM  OISP  CUM  TANG  Dt$»  1  SMEARED* I 

1000  FORNAT  113* 2F 4.2. 6E 15.3* 110) 

END 


APPENDIX  2 


JOINT  INFLUENCE  ANALYSIS  CCMHJTER 
PROGRAM 

FORTRAN  PROGRAM  TO  GENERATE  STRESSES  BY  FIRSCH  SOLUTION  FOR  GIVEN  M« 
FIELO  PRESSURES  FORCING  A  CIRCULAR  TUNNEL  AND  TO  CALCULATE  AND  PLOT  THE 

zones  or  influence  around  the  tunnel  or  joints  or  a  own  orientation 

AND  A  GIVEN  COHESION  ANO  FRICTION.  PRODUCES  A  PLOT  ON  ON-LINE  PRINTER 


PROGRAM  FIRSCH  1 INPUT.OUTPUT ) 

C 

COMMON  S1GX,  SIGY.  SIGZ  •  TAUXZ.EACTJT 10) »  C.  PCLUOt*  DCL2T10). 

1  OCNIIO}  .  OOJ2UOI.  OCNIIOI*  HCOI12I*  ITV»CM0U10)  • 

7  ROVCftA!44>]01)»  THETA (44)  .  |PL0T(44.3.10t«  PPLOT  U6»S.lf )  ♦ 

3  FMAXT44.10I.  ITH.  IN,  RADCNO*  NUMJftt.  PKl» 

DIMENSION  CONTORT  SI .ectatioi 
REAL  1PLOT 
C 

c  READ  ANO  PRINT  OP  CONTROL  INfORNAT |ON 

C  READ  ANO  PRINT  OF  CONTROL  INFORNATION 

Iff  READ  lOOO.lKEOTll.l-ltlZl.NP 
IF  TNP.GT.O)  GO  TO  911 

READ  2010.  NUMJ7M.  RADIX.  RADCNO.  THlMC*  PJ.  P2»  GNU.  S.  C*  PHI. 
1  St GYl .SETAT 1 | 

PRINT  1001.  IMEOm.l-l.17l.NUNjTM.RAOlNC.RAOENO.TMlNC.Pl.P2* 

1  OMfS.C.PH! 

PRINT  lOOa.MGYS.SFTATtl 
C 

C  DETERMINE  DIRECTION  COSINES  FOR  BETA  ORIENTATION 
0CNT1 [-COST  BE  TAT ll*3.]4|S9/lt0.0l 
DCN2T  l )-f 1.-0CMT1  I***)/! l.*S**2l 
0CNUI-S0RTJ0CN7T1II 

OCUll*-S«OOl<ll 

DCL2<H*S**2'OCR*MI 

C 

C  GENERATE  SIGF  ANC  S1G  THETA  FROM  XfRSCH  SOLUTION 
IR-IFIX!  TRAOENO-1* I /RAO I  NO 41 
|IH-IF|XT90./TH1NC)41 
P1PP2  •  TPl  ♦  P71  •  0.5 
P1HP2  •  TPl  -  Pj)  •  0.5 
PIR  •  3.1415*27  /  UO. 

PHIR  •  PHI  •  P|R 
00  1  JMtITH 
JPTJ-ll  2.2.3 

2  TlfTAT JI-OnO 
GO  TO  4 

3  TMCTATJ1  •  THETAT J-l I ♦  THlNC 

4  TMCTAR  •  PIR  •  THETAIJI 
THCT2  •  2.0  •  THETAR 

S2  •  SlNTTMCTAR)**? 

C2  •  COSTTMETARi#*? 

SC-  SINT  THETAR I *COST  THETAR I 
00  11  I-l. IR 
IFlt-11  3.5.4 

5  ROVERAT J. I 1-1.0 


GO  TO  7 

4  ROVERAT J. 1 1- ROVERAT  J. I-l HRAOINC 
7  RAD2  •  Tl.O/ROVERATJ.m  ••  2 
SIGRA  •  P1PP2  •  U.O  -  RAD21 

SIGRB  •  P1NP2  •  11.0  -  RAD2*I4.0-3.0vRAD2I I  ®  C0STTHET2I 

SIGN  •  SIGRA  ♦  SIGRB 

SIGTHA  •  PJPP2  •  11.0  ♦  RA021 

SIGTHB  O  P)MP2  •  11.0  ♦  3.0#RAD2*RA02)  •  C05TTHET2) 

SIGTH  •  SIGTHA  -  SIGTHB 

TAURTH  •  -PJHP2  •  (1.0  ♦  RAD2M2.0-3.0+RA0?) I  •  SINTTHET?) 

C  OBTAIN  PRINCIPAL  STRESSES 

R-S0RTII!StGTH-$!GR)/2.)*«?*TAURTH4«2) 

StGM-TSIGTH4S!GRl/2-0 

SIG1-SIGM.R 

SIG3-S1GM-R 

SIG2-0.0 

C  OBTAIN  XYZ  STRESSES 

SIGX-SIGR*C2  -  7.»TAURTH4SC  ♦  SI6TH«S2 
SIG2-  SIGR*S2  ♦  2.*TAURTH»SC  ♦  StGTM*C2 
TAUX2  — TS1GTH-SIGR|«SC  «TAURTH»!C2-S?1 
SIGY-SIG2  ♦  SIGY1 

C  OETCRHINE  FACTOR  OF  SAFETY  FOR  SLIOING  ON  JOINT 
C 

00  101  L*1.NUMJTM 
CALL  JOINTTL) 

C 

C  ASSIGN  A  VALUE  TO  I  TYPE*  1  IF  IN  ZONE  OF  SLIP.  7  IF  OUTSIDE  ZONE 
C 

IF  TFACTJTL)  -  1.0)  BO.  10.  01 
•0  ITYPETI.D  •  1 
GO  TO  101 
•  1  ITVPfll.U  •  2 
101  CONTINUE 
C  END  OF  IR  LOOP 

l l  CONTINUE 
C 

C  ESTABLISH  COORDINATES  FOR  CONTOURS 

C 

00  31  L-l.NUNJTN 
K  •  1 

DO  30  1-2. IR 

IF  TITVPEII.L)  -  ITYPETI-I.U1  21.  30.  23 
21  IPIOTIJ.X.LI  •  1H2 
GO  TO  24 

23  IPLOTT J.K.L1  "  IHf 

24  RPtOTTJ.K.L)  -  ROVERAT J.I-II 
K  •  K»t 

SO  CONTINUE 
31  KMAXI J.L >  •  K-l 
C 

C  CNO  OF  J  LOOP 
1  CONTINUE 
CALL  PLOT 
GO  TO  If* 

999  RETURN 


1B3 


1000  FORMAT  fi?A4.  I7» 

1010  FORMAT  t I10.3F10.4/  4F10.4) 

tOOl  FORMAT I lH1.12A4/2*HGNUNRfR  OF  ORIENTATIONS  IX. 13/ 

1  13H01NCREMENT  OF  RADU  4X.riO.4Z  1 7HD0UTERM0ST  RADIUS  0X.F10.4Z 
213M0INCREMENT  OF  TmETA  4X.flO.4/20NO:)ORIZO*TAL  PRESSURE  4X.F1C. 
32Z13HOVERTJCAL  PRESSURE  7X.F10.?/16K0POlSSON»S  RATIO  3X.F10.4Z 
414H0SL0PE  OF  X2  TRACE  7X.F10.4Z 

3  WOCOMfSION  23X.no.4Z  umoamgle  of  friction 

4  13X.F10.4l 

1007  FORMATUH1.3X.2H  L  4X.  2H  M  IX.  2H  N  | 

1003  FORMAT!  1HOOF1D.4I 

1004  FORMAT! 1M0/  10H  SIGYJ  •  F10.4/11H0BCTAT 1 )  •  3X.F10.il 

1020  FORMAT T4C 10.4 1 

END 

SUBROUTINE  JOINTTL) 

C 

COMMON  SIGZ  SIGT.  SIGZ.  TAUXZ.FACTJ! 10).  C»  DCLT10).  DCL2T10I. 

1  OCN»  •>)) .  DCN2I10I.  OCMtlO).  ME0T121.  ITYPCUOl.IOt. 

2  ROVERAT 44*101 1 .  THETA T46 1 »  IPL0TT44.3.10I •  PPLOTT44.5.10I . 

3  KMAXT44.10).  ITH.  IR.  RAOENO.  NUNJT*.  PhIR 
C 

C  OBTAIN  SHEAR  ANO  NORMAL  STRESSES  ON  JOINT 
C 

S IGV.0CL2 T D *S1GX*0CMT LI ®*?*SI6Y*  DCN2!D*S;gZ*2.«OCN!L l*OCL T L)« 
lTA'JXZ 

TRACT2»TXL«U*SIOX*OCNTLl«TAUXZ)#«24TOCMIC»«SlGri«4?4lOCLTLJ« 
1TAUZ2*0CN!1)*S!GZ)«*2 
YAU  •  SORTCA8STTRACT2-S!GN**2)1 
C 

C  FINO  FACTOR  or  SAFETY  IN  JOINT 
C 

SHSTR-C4S!GN*TANIPm|R| 

IFTTAUI  21.21.72 

21  FA'TJTLI-I.EMO 
GO  TO  23 

22  FACTJTL)*SHSTR/TAU 

23  RETURN 
CNO 

SUBROUTINE  PLOT 
C 

COMMON  SIGX.  SIGY#  SIGZ.  TAUKZ .FACT JT 10) .  C*  DfLMOl.  OCL2UO*. 

1  OCN! 10 I .. 0CN2T 10) .  DCMT10).  H?OT 121 «  ITVPET 10! .101 . 

2  ROVERAT 46* 101 1 .  ThETA!46).  JPL0TU6.5.10) •  PPLOTT44.3.10). 

3  KKAXI 44*10).  ITH.  |R.  RADENO.  NUHJTM.  PHIR 
DIMENSION  AT  101} 

REAL  I PLOT 
C 

IH  •  ITH  ♦  1 

TINES  •  IR  /  TRAOENO  -  1.01 
C 

C  initialize  printing  lihf  FOR  porocr 

00  *00  L*1 .NUNJTM 
*0  00  91  1*2.100 

•  1  A( 1 1  •  IW- 

00  92  1*1.101.3 


92  ATI)  •  IH. 

KOUNT  •  3 
C 

C  PRINT  TITLE  ANO  UPPER  BOROER 
ASSIGN  102  TO  NEXT 

100  PRINT  101.  HCO.  OCLTL).  OCMTL) .  OCNTL) 

101  FORMAT  T1H1Z/Z23X.12A4/M6X.3HL  -F10.4.7H  M  «F10.4*7H  M  -F10 
1.6ZZZI 

GO  TO  NEXT 

102  PRINT  103 

103  FORMAT  T 8X.IM123X.3H1 «S73X.1H??)X.3H7.323X.1H3) 

C 

300  00  320  JX-l.ITM 
J  •  IH  -  JX 

IF  TKNAXIJ.LI)  303.  303.  301 
30!  XX  «  KMAXTJ.L) 

00  302  F-l.KF 

I  -TRPLOTT J.F.LI  -  1.0)  •  TIMfS  ♦  o.5 
102  AT!)  •  IPLOTT J.K.LI 

303  IF  TXOUNT  -  3)  304.  306.  320 

304  PRINT  303.  A 

305  FOCMAT  (IX.IOIAI) 

FOUNT  »  FOUNT  ♦  I 
GO  TO  304 

306  PRINT  307.  THETAT Jl.  A*  THETAIJI 

307  FORMAT  I3X.F2.0»1X*!01A1*F3«0I 
FOUNT  •  | 

30|  IF  TJ  -  2)  320.  310.  313 

310  00  311  ’-2.100 

31)  ATI)  •  1H- 

00  312  1-1.101.3 
312  ATM  -  IH* 

GO  TO  320 

315  00  314  1*2.100 

316  AH)  ■  2H  Z 
320  CONTINUE 

PRINT  103 
C 

C  PRINT  VALUES  OF  PlOTTCO  POINTS 
ASSIGN  400  TO  NEXT 
GO  TO  100 

400  PRINT  401 

401  FORMAT  l73X»',HTHETAnx.'/T3HR/A13XM 
00  *10  JM.ITH 

IF  TFNAXIJ.LI)  4|0*  410.  402 

402  FF  -  FMAXIJ.U 

PRINT  403.  TMETATj).  TRFLOTT J.F.L) .  F*1*RFI 

403  FORMAT  U0X.6F14.7I 
410  CONTINUE 

C 

900  CONTINUE 
RETURN 
ENO 
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^Rock  and  Joint  properties  vere  reviewed  as  a  basis  for  developing  the  analyses. 
The  review  discussed:  failure  criteria  for  the  rock;  and  joint  properties.  Improved 
methods  of  sampling  and  testing  joints  vere  explored  in  detail.  Pile  Driver  test 
sections  vere  used  as  a  basis  for  analysis. 

Two  approaches  of  analysis  of  the  rock  to  predict  performance  of  rock  bolted 
and  unlined  test  sections  were  used.  The  first  method  called  finite  element  was 
used  on  one  section  only.  This  method's  complexity  prevented  its  application  to 
other  sections.  Discontinuities,  heterogeneities,  and  local  failures  can  be  taken 
into  account  by  finite  elements,  but  the  representation  is  only  two  dimensional. 

Joint  influence  analysis,  the  second  method,  was  applied  to  many  sections.  In  com¬ 
parison  to  the  finite  element  method,  Joint  influence  analysis  is  relatively  simple 
to  use.  Joints  and  rock  bolts  are  assumed  not  to  modify  the  stress  distribution  in 
the  joint  influence  analysis;  however,  the  method  represents  a  three  dimensional 
approach . 

Comparisons  were  made  to  post-shot  cross  sections  of  the  Pile  Driver  openings 
to  evaluate  the  methods  of  analysis.  The  finite  element  method  is  more  powerful; 
the  Joint  influence  analysis  ia  easier  to  use.  > 
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